Design and performance of a Collimated Beam Projector for telescope transmission measurement using a broadband light source

Type Ia supernovae are the most direct cosmological probe to study dark energy in the recent Universe, for which the photometric calibration of astronomical instruments remains one major source of systematic uncertainties. To address this, recent advancements introduce Collimated Beam Projectors (CBP), aiming to enhance calibration by precisely measuring a telescope's throughput as a function of wavelength. This work describes the performance of a prototype portable CBP. The experimental setup consists of a broadband Xenon light source replacing a more customary but much more demanding high-power laser source, coupled with a monochromator emitting light inside an integrating sphere monitored with a photodiode and a spectrograph. Light is injected at the focus of the CBP telescope projecting a collimated beam onto a solar cell whose quantum efficiency has been obtained by comparison with a NIST-calibrated photodiode. The throughput and signal-to-noise ratio achieved by comparing the photocurrent signal in the CBP photodiode to the one in the solar cell are computed. We prove that the prototype, in its current state of development, is capable of achieving 1.2 per cent and 2.3 per cent precision on the integrated g and r bands of the ZTF photometric filter system respectively, in a reasonable amount of integration time. Central wavelength determination accuracy is kept below $\sim$ {0.91} nm and $\sim$ {0.58} nm for g and r bands. The expected photometric uncertainty caused by filter throughput measurement is approximately 5 mmag on the zero-point magnitude. Several straightforward improvement paths are discussed to upgrade the current setup.


INTRODUCTION
Type 1a supernovae (SNe Ia) are a class of supernovae that are thought to result from the explosion of a white dwarf (WD) star in a binary system (see Hillebrandt & Niemeyer (2000) and Mazzali et al. (2007) for detailed reviews of SNe Ia physics).Their uniform intrinsic brightness at the 15 per cent level after standardization and their luminosity around 10 10  ⊙ make them particularly useful for cosmology, in particular, to accurately measure the expansion rate  () and to study the properties of dark energy with the determination of the equation of state parameter  (Garnavich et al. 1998;Efstathiou 1999;Copeland et al. 2006).As the number of surveys increases, the statistical uncertainty on cosmological parameters measured with SNe Ia diminishes (Brout et al. 2022).For instance, the Vera Rubin Observatory will observe over 10,000 SNe Ia in its 10-year survey (LSST Science Collaboration 2009).Systematic uncertainties thus become predominant, as demonstrated by Scolnic et al. (2018).These include SNe Ia models, bias selections, galactic extinction, peculiar ★ E-mail: kelian.sommer@umontpellier.fr velocities, local environment, and not the least calibration errors.Indeed, instrumental photometric flux calibration represents the main source of uncertainty for the estimated distances of SNe Ia (Stubbs & Brown 2015) and for cosmological parameters inferred with this probe.We refer the reader to Goobar & Leibundgut (2011) for a review of SNe Ia cosmology and to Betoule et al. (2014); Scolnic et al. (2014) for a review of calibration uncertainties.
This flux calibration issue can be split in two different components: (i) the terrestrial atmospheric transmission, and (ii) the wavelengthdependent instrumental throughput.The first component can be determined by comparing the measured flux of reference stars in the field of view for each photometric band with catalogs of spectrophotometric reference stars.These ultimately rely on 3 pure hydrogen white dwarfs primary CALSPEC standards (Bohlin et al. 2020) for which physical models describe very accurately the spectra (Narayan et al. 2019).Models are internally consistent with Hubble Space Telescope photometry to better than 1 per cent (see Bohlin et al. (2014) for an extended review of absolute calibration methods).To isolate the uncertainty component due to the instrument throughput, multiple efforts with in situ artificial methods have emerged for dif-ferent surveys (Stubbs & Tonry 2006;Stubbs et al. 2010;Doi et al. 2010;Marshall et al. 2013;Regnault et al. 2015;Lombardo et al. 2017).They involve laboratory metrology standards with extremely well-characterized photodetectors.Stubbs & Tonry (2006) proposed a method utilizing tunable lasers to characterize the relative wavelength response of a complete imaging setup relying on a flat-field screen and a calibrated reference photodiode.It was later implemented to determine the relative throughput of the PanSTARRS telescope and the Gigapixel imager with total uncertainty higher than the per cent level, all systematics considered (Stubbs et al. 2010).Doi et al. (2010) developed a monochromatic illumination system -which consisted of a broadband lamp, a monochromator, an integrating sphere and a photodiode -dedicated towards in situ measurements of the response function of the mosaicked CCD imager used in the Sloan Digital Sky Survey (SDSS).They were successfully able to quantify response variations to be lower than 0.01 mag for SDSS  bands over the entire duration of the survey.Their conclusion stated that the main driver behind throughput evolution with time was the seasonal variations described by temperature effect.Therefore, as it is planned for the Vera Rubin Observatory, the throughput measurement operation should be performed regularly across the year.Marshall et al. (2013) developed the DECal instrument which provided both broadband flat fields and narrowband spectrophotometric calibration for the Dark Energy Camera (DECam, see Flaugher et al. 2015) installed on the 4-m Blanco Telescope of the Dark Energy Survey (DES).The flat field system employed LEDs to illuminate individual photometric filters, while the spectrophotometric calibration system involved a standard tunable light source from a monochromator that was directed onto the flat field screen using a specialized fiber bundle and diffuser.Calibrated photodiodes placed along the beam track the telescope's throughput across wavelengths.They successfully deployed a prototype at the Swope 1-m and du Pont 2.5-m telescopes at Las Campanas Observatory in Chile and measured the throughput of the   filters used in the WIRC and RetroCam instruments during the Carnegie Supernova Project with an accuracy of 1%.The DICE (Direct Illumination Calibration Experiment, see Regnault et al. 2015) instrument was composed of multiple narrow-spectrum LEDs ranging from UV to NIR spectral range emulating a point-like source placed at a finite distance from the telescope entrance pupil, yielding a flat field illumination that covered the entire field of view of the imager.They were able to measure the MegaCam and SkyMapper  passbands normalisations terms relative to  band with lower than 0.5% relative uncertainty.Lombardo et al. (2017) built the SNIFS Calibration Apparatus (SCALA) system for the SuperNova integral field spectrograph (SNIFS) which used a monochromatorbased instrument with integrated spheres whose output exit ports were collimated and projects large spots onto the focal plane array of the SNIFS instrument.They assessed the long-term repeatability of the light emitted by SCALA to be better than 1% and evaluated the color-calibration stability of the University of Hawaii 2.2-m (UH 88) telescope, which remained at the per cent level over one year across the wavelength range from 400 nm to 900 nm.
We focus here on the development of a calibration device, called a Collimated Beam Projector (CBP), designed to accurately measure the wavelength-dependent response of the instrument   (), including the optics, the CCD camera, and the filters.A first design was proposed by Coughlin et al. (2016); Coughlin et al. (2018).Their prototype was tested on the CTIO 0.9-m telescope and the DECam wide field imager.The principle and main advantage of the CBP compared to previously cited instruments, is that it relies on illuminating one or multiple spots of the primary mirror with a stable collimated monochromatic light beam.The incident wavefront is parallel upon entering the telescope, resulting in the creation of a collimated beam that covers a fraction of the entire pupil.The full-pupil illumination is synthesized by scanning at multiple angles and radii.Consequently, this design minimizes reflections and stray light within the optical setup, which is a considerable improvement over utilizing a flat-field screen for illumination (Coughlin et al. 2016).It also facilitates the separation of direct illumination beams and ghost spots for analysis and computations.Finally, it measures the output on the detector, which allows the calculation of the throughput.The total transmission function is then rasterized by interpolating between angles and spot positions.
Recently, Mondrik et al. (2023) designed another CBP prototype to measure the StarDICE (Betoule et al. 2023) proof-of-concept 254mm telescope, achieving a ∼ 3 per cent accuracy over the 300 -1000 nm wavelength range.Souverin et al. (2022) built an enhanced version of the CBP for the final StarDICE experiment (Fig. 1) and were able to measure the response of a 400 mm F/4 Newtonian telescope coupled with ugrizy filters and CCD sensor at the 0.1 per cent statistical uncertainty level.This version required a powerful class 4 laser source, unpractical for transportation and exploitation due to its fragility and security concerns.The same kind of design will be used for the Vera C. Rubin Observatory CBP currently in construction, with a laser packed inside a dedicated thermally controlled enclosure (Ingraham et al. 2022).For the Dark Energy Science Collaboration (DESC) 10-year supernovae experiment (The LSST Dark Energy Science Collaboration et al. 2018) of the Vera C. Rubin Observatory, systematic uncertainties on the griz bands zero points and central wavelengths are required to be lower than 1 mmag and 0.1 nm respectively.This is undertaken with the aim of harnessing the capabilities of the instrumentation and survey, leading to substantial enhancements in cosmological constraints.These metrics evaluated on the Zwicky Transient Facility (ZTF) photometric system bands ztf-g and ztf-r (Bellm et al. 2019;Dekany et al. 2020) will serve as benchmarks throughout the remainder of this study.Given that ZTF neither possessed nor has plans to acquire a permanent CBP for bandpasses throughput measurements, there exists a unique opportunity to design an innovative solution that fulfills an actual requirement.
Therefore, this work focuses on the first attempt to develop a "traveling CBP", a transportable version of a CBP for on-site measurements at multiple observatories.We substitute the sensitive tunable laser with a more conventional light source.The primary motivation for switching to broadband sources is their simplified maintenance, cost-effectiveness, and reduced susceptibility to environmental flaws compared to lasers.This versatility renders them more adaptable than lasers, the challenge now being to reach a sufficiently high monochromatic brightness.
The paper is structured as follows.In Section 2, we describe the architecture of the setup and explain the role of each component.In Section 3, we detail the data acquisition procedure, and the noise model and provide a complete analysis with forward modelling techniques.Section 4 presents the CBP throughput curve and other relevant metrics.In Section 5, we address the budget measurement uncertainties including systematics.In Section 6, we discuss the results, potential future improvements, and methods to mitigate some measurement and systematic uncertainties.We draw our conclusions in Section 7.

EXPERIMENTAL SETUP
The experimental setup used in this work is shown in Fig. 2, where the elements of a CBP (Fig. 1) can be recognized.In order to calibrate the CBP telescope and its attached photodiode (PD), we illuminate a solar cell (SC) at the exit of the CBP telescope.This SC has been calibrated using NIST standard calibrators (Larason & Houston 2008;Brownsberger et al. 2022) and has served for the StarDICE telescope calibration.Our goal is to calibrate the throughput of the CBP telescope, using another source than the laser, and relate it to the PD current.The main challenge stems from the fact that the solar cell requires a high illumination intensity, easily reachable with a class 4 laser, but not necessarily with a more standard light source.
In this section, we provide technical details of each element in the instrumental setup.

Traveling CBP light source prototype
We employ a standard broadband 175W Xenon arc lamp (Spectral Products ASB-XE-175) of 5600K color temperature, built with a Xenon parabolic lamp.Its emission spectrum, which lies in the range 200-2000 nm, is shown in Fig. 3.This broadband light is then focused into a Czerny-Turner (Cornerstone 130 1/8 m) monochromator entrance slit with an uncoated lens located at a known distance of the entrance slit.The focal length of the monochromator is 125 mm for an acceptance aperture of F/3.7.The grating available in the monochromator used in this test is a 30x30 mm 2 with 2400 grooves/mm and a blaze wavelength of 275 nm.The manufacturer's documentation recommends setting the entrance and exit slits at approximately 1 mm width for a 2-3 nm spectral bandwidth.Slit heights are adjusted to the maximum of 12 mm to collect as much light as possible.Multiple diffraction orders can coexist in the beam at the exit slit.Therefore, a long pass filter designed to block undesirable higher diffraction orders is inserted between the exit slit of the monochromator and the chopper wheel (introduced below).The transmission curve of this Thorlabs 550 nm Hard-Coated Edgepass Filter, as provided by the manufacturer, is shown in Fig. 3.
Next, a chopper wheel, made of a simple stepper motor driven by an Arduino Nano micro-controller coupled with a ULN2003A unipolar motor driver, rotates a custom 3D printed mask with two apertures in front of a fixed slit (see Fig. 4).Its purpose is to significantly decrease the level of 1/  noise that originates from small fluctuations of the operational amplifier bias voltage which can draw current through the relatively low shunt resistance of the solar cell.The amplitude of the contamination is inversely proportional to the resistance of the solar cell used in this study.The motor rotates at a rate of ∼ 3 Hz.The full analysis that consists of reducing the 1/f noise contribution, accounting for small rotation speed variations and aperture differences is presented in Section 3.

Integrating sphere and spectrograph
We use the General-Purpose Ø50 mm Integrating Sphere (IS) from Thorlabs in order to ensure identical beam fluxes toward the sensors.Its typical reflectance is better than 98 per cent in the range 350-1500 nm (see Fig. 3).A spectrograph connected to the IS through an optic fiber cable monitors spectral shapes of the monochromatic beams exiting the monochromator slit.

Photodiode
We use the SM05PD1B silicon-based photodiode (PD) from Thorlabs, designed to operate in the UV-VIS-NIR range (350-1100 nm), in order to measure the incoming light source intensity.It has a photosensitive area of 3.6 × 3.6 mm 2 .The PD operates in photovoltaic mode to reduce dark current which is important to optimize performance with this type of low-illumination experiments.Electrical charges are sampled by a Keithley 6514 at 50 Hz rate, in conjunction with a 1 power line cycle to get rid of the 50Hz power line influence.
During an integration period, the PD achieves charge levels ranging from 10 to 50 nC.The electrometer's bias current is less than 4 fA, which is negligible in comparison to average photocurrent levels of ∼ 20 pA.A 1 mm hole is placed in front of the photodiode in order to ensure that measurements are within the limits of the electrometer measurement range and do not saturate.(Brownsberger et al. 2022), all the data is obtained from the specifications provided by manufacturers.

CBP telescope and Solar Cell
The CBP telescope was chosen to be compact and with the simplest optical layout in order to avoid internal reflections with field correctors or chromatic aberrations with lenses, which lead to the appearance of ghost patterns (see Murray (1949) for an in-depth explanation).It is an Omegon Ritchey-Chrétien equipped with two hyperbolic BK7 mirrors averaging 92.7 per cent reflectivity.It has a 152 mm aperture and a 1370 mm focal length.The primary mirror is larger than the SC.A mask is mounted at the entrance of the telescope in order to only illuminate the SC photosensitive surface.Multiple tinfoil layers are added to ensure that the assembly is fully light-tight in the optical and infrared range.
The solar cell is a single high-efficiency Sunpower C60 version already studied in Brownsberger et al. (2020Brownsberger et al. ( , 2022)).This SC offers a high quantum efficiency (QE) > 90 and a large collecting area > 10 × 10 cm 2 at a low cost with great portability.Its output charge is measured by a Keysight 2987A multimeter at a sample rate of 500 Hz, based on the acceptable SC time response (at least 200 Hz).Brownsberger et al. (2022) review the absolute calibration procedure relative to the NIST at 0.1 per cent accuracy of the QE between 350 and 1100 nm.One important downside of the SC is the strong dark signal caused by a low shunt resistance  ℎ , which is due to manufacturing defects and causes power losses in solar cells by providing an alternate current path for the photocurrent.Our model of SC was measured with  ℎ = 1.8 kΩ, to be compared with the serial resistance   = 170 mΩ, which is the electrical resistance that hinders the flow of current through the main current path in a solar cell.It is the same SC sample as the one used in Souverin et al. (2022).

DATA ACQUISITION AND ANALYSIS
The monochromator was used to scan wavelengths from 350 nm to 750 nm with a 2 nm step.With the spectrograph, spectra were obtained at five different wavelengths across the considered range of measurements to assess the quality of the monochromatic light beam (shown in Appendix A).We correct the wavelength values read on the monochromator using a third order polynomial fit to the 5 values from the spectrograph.
At each wavelength, the charge accumulation from the SC and PD was read during a 20 s integration time, resulting in 10 4 and 10 3 values, respectively.The monochromator shutter was closed, and the electrometer charges reset between each integration of 20 s.Intertwined to this "signal" sequence, 4 sequences of dark current exposure measurements were performed by leaving the monochromator's shutter closed during the 20 s integration time.We compute the average of the dark current for the PD and subtract it from each time series as it is estimated to behave like stationary noise.The sequence of charges was finally turned into time series of currents by finite difference, using the timestamps provided by the electrometer's clock.As an illustration, Fig. 5 shows in red the 20 s signal in the PD for the  = 550.5 nm setting.The 6 Hz periodicity induced by the chopper is visible.
For the analysis, we adopted a purely modeling-based approach, favoring it over a lock-in detection method for several reasons.Firstly, the prototype bench had a limited existence of less than three weeks in Summer 2022, constrained by material and time limitations.Secondly, the lock-in detection method demands more effort for signals that deviate from being squared or sinusoidal.Lastly, its implementation becomes challenging when dealing with sensors and electrometers with varying bandwidths.The decision to configure the electrometers in charge mode was influenced by replicating methods from previous CBP versions (Souverin et al. 2022;Mondrik et al. 2023).This choice aimed to better account for the laser's pulse-topulse instability, often reaching 10% or more.The chosen hardware integration allowed them to avoid worrying about the fact that time sampling (50 or 500 Hz) was well below the photodiode response time (µs) to laser pulses (ns).Despite experimenting in current mode on the Keysight, where the SC is connected, contributing significantly to the uncertainty budget, no observable improvement was noted.

Photodiode noise model
The PD is affected by readout noise that we estimate experimentally by taking the standard deviation of 4 concatenated dark exposures:    = 0.494 pA.It amounts at worst to 5 per cent of the PD signal in our dataset.This encompasses the combination of both dark noise and electrometer read-noise.
Moreover, despite having low impact on this experiment, the shot noise needs to be accounted for (Chen 2018): with  ℎ () the shot noise in A,  the elementary charge of the electron  = 1.602 • 10 −19 C, Δ  the bandwidth of measurement in Hz (i.e. the frequency range over which the noise is being considered to compute the noise power),  ℎ  () the photocurrent at time  (i.e. the current produced by the photodetector due to incident photons),    the dark current (i.e. the current in the absence of light, often caused by thermal effects or other background sources and taken as the average of 4 dark exposures).Therefore, the total PD photocurrent noise at time  for the scanned wavelength  is estimated as Overall, photocurrent readings are dominated by the readout noise representing on average ∼ 85 per cent of   .Fig. 5 shows the different noise contributions.

Solar cell noise model
The noise structure of the SC is much more complicated, as can be seen in Fourier space on Fig. 6: 1/  and electrical grid contributions largely dominate third and higher orders amplitudes.
We build a noise model for the SC (see equation 3), noted S d (f,   ) with the vector of parameters   = {A, , , } (shown in black on Fig. 6) by fitting in Fourier space the averaged spectrum Ŝ (shown in blue) of 4 SC darks (shown in gray).Hereafter, the overline notation specifies the average whereas the hat indicates an experimental measurement.
where A and  are the 1/  noise parameters, f is the vector of frequencies,  and  account for a small departure from white noise, and   is a Voigt function modeling the 2 spikes at  = 50 Hz and 150 Hz: where   is the peak width and   is the peak amplitude.This model is fitted to the Ŝ data, resulting in the error model subsequently used for the SC signal.
As can be seen on the bottom panel of Fig. 6, the fit is reasonably good, considering the inherent difficulties in modeling all the noise contributions.Table 1 shows the best-fit values with their associated 1- uncertainties.

Photocurrent model
As mentioned in Section 2.1, the chopper wheel shapes the intensity into a periodic signal of frequency , which we model as follows.
First, the time  is folded into a phase  ∈ [0, 2[, where  is the frequency or period of the periodic signal, and  is a phase offset to centre the signal on the [0, 2[ interval.The intensity received through one slit is then modeled as, with( 5) ( ; ,  1 ,  2 ) = ( ; ,  1 ) − ( ; , 1 −  2 ) and ( 6) Equation 7 is the usual sigmoid function.In Equation 6,  2 is defined so that its value is close to  1 .Equation 5 ensures that ℎ evaluates to 1 at  = 0.5, and to 0 at  = 0 and  = 1.Finally, the modeled intensity is where the first function accounts for the first slit of the chopper seen when a run starts ( between 0 and 1) and the second function accounts for the second slit of the chopper ( between 1 and 2) and has the parameters  inverted.The parameter µ accounts for a fluctuating baseline.The amplitude , in pA, is the main parameter of interest to determine.

Signal fitting
For each wavelength step, a  2 -minimization fit of the signal model is performed with lmfit1 on the PD and SC time-series, in time space for the PD and in Fourier space for the SC, using respectively Equations 2 and 3 to define the covariance matrix used in the fit.Though it should be very close to 6 Hz, the signal frequency  is fitted to the data for better accuracy.For the PD, it appears that the baseline parameter  shows a behavior correlated with the intensity of the light, which is not completely understood; as a result, we modify Equation 8 so that the amplitude A is factoring both the two ℎ functions and .Fig. 7 shows the PD data and fitted model for a wavelength corresponding to a SNR in the middle of the dynamic range.One can observe that the frequency is indeed very close to 6 Hz, and that the residuals, though quite small (at most 6 per cent of the signal amplitude ), are structured in a way suggesting some imperfection related to the beam shape being blurred and having tails.Fig. 8 shows the SC data and fit result to observations at a wavelength with medium SNR.The much degraded level of noise in time space, when compared to the PD, is apparent, together with the digitization of the Keysight's readings.

RESULTS
Performing the fits to the PD and SC signals at each wavelength yields the results presented in Fig. 9 and 10, respectively.Some comments are in order.First, the wavelength dependence of the signal amplitudes mostly follows that of the light source.The discontinuity at 570 nm is the result of the addition of the high-pass filter discussed in Section 2.1.Although being without consequences due to the error budget being largely dominated by the SC (see Section 5.2), the fact that the reduced  2 for the PD follows the same trend is indicative of an incomplete error modelization.
We can now proceed with the main results of this work.We first determine the relative throughput curve of our setup.Then, to assess its performance for ZTF photometry, we explore the effect of uncertainties on the zero-point   and the central wavelength  of the ztf-g and ztf-r filters.

Determination of the CBP telescope throughput
The QE of the SC was measured against a NIST-calibrated photodiode, to 0.1 per cent relative uncertainty (Brownsberger et al. 2022).
It therefore measures accurately the number of photons emerging from the CBP.By comparing the ratio of the SC photocurrent to the PD current, we can use the latter as a proxy to retrieve photon flux on other telescopes and determine their throughput 2 .Hereafter, we consider the normalized ratio   /  depicted on Fig. 11.It is assumed that all flux loss between devices is due to the telescope itself (i.e.mirrors with coating, internal baffling).The total statistical uncertainty on the CBP telescope throughput measurement is expressed as, () and    () arise directly from the weighted leastsquares fit presented in Section 3.
For the  filter (see spectral range illustrated as a green area in Fig. 11), we have an average statistical uncertainty close to ∼ 1 per cent for each data point.It increases to ∼ 2 per cent for the  filter (red area in Fig. 11).SC noise overwhelmingly dominates the total statistical uncertainty budget.The PD contribution to the statistical uncertainty is ∼ 0.017 per cent relative uncertainty.The increase seen near 350 nm is due to the decrease in light source brightness.In contrast, the increase seen around 600 nm is due to the low monochromator grating efficiency.Souverin et al. (2022) showed with higher-resolution measurements that the telescope transmission is smooth, as are the quantum efficiencies of the PD and SC.We can use this behavior to interpolate the data with a smooth function, reducing the global statistical error on the response curve.The solid black line is a cubic spline fit.Residuals between the observations and the spline interpolation are shown in Fig. 11 and no structure seems to appear.It would be valuable to conduct repeated measurements to examine the consistency of observations around the fitted smooth spline curve and improve the fit quality.

Errors on magnitude estimation
We can estimate the impact of the uncertainty of the CBP throughput measurement by calculating its effect on photometric measurements.We first compute the effect of the uncertainty on the filter transmission curve on magnitudes.The magnitude   of a flat spectrum measured with a photon counting device through a filter  of transmission   () is with ℎ the Planck constant and  the photon frequency.To evaluate the effect of the CBP throughput error, we draw random samples from a normal distribution N ( =   (  ),  =    (  )) that we inject onto the magnitude equation, and calculate the mean and standard deviation   , of the magnitude distribution we obtain.We propagate uncertainties on the magnitude zero-point for a flat spectrum.For the ztf-g and ztf-r bands respectively, we obtain similar zero-point statistical uncertainties with   , = 0.0049 mag and   , = 0.0052 mag.

Central wavelength uncertainty
The second key metric is the accuracy of the bandpass mean wavelength, i.e. the central wavelength  for a flat spectrum: Here, estimations are done in the same manner as for the zero-point magnitude uncertainty estimation   , by drawing random samples from a normal distribution.Both filters appear determined with subnanometer accuracy.For the ztf-g and ztf-r bands respectively, we obtain    = 0.9 nm and    = 0.6 nm.As the filters have large bandpasses Δ  =  , −  , (135 nm and 160 nm for ztfg and ztf-r filters respectively), relative uncertainties /Δ  of central wavelengths are less than 0.5 per cent for both filters.

EVALUATION OF UNCERTAINTIES
In our efforts to uphold an accurate uncertainty budget for the CBP system, we have made diligent attempts to address the potential effect of the model error.In this section, we establish an assessment of the uncertainty budget and explore three sources of systematic uncertainty that could impact our results.

Overview of potential systematic uncertainties
As our primary focus was to demonstrate the feasibility of the concept itself through the realization of this non-definitive prototype, we did not conduct an exhaustive qualitative and quantitative evaluation of each systematic effect.However, it is essential to recognize that these factors may warrant more in-depth investigation in future iterations of this work to refine the accuracy of our measurements.Below, we detail two main systematic effects affecting this particular experimental setup.
Ambient and stray light contamination : the monochromator manufacturer datasheet3 indicates a 0.03 per cent stray light level.The longpass filter manufacturer datasheet4 gives an optical density (OD) in excess of 4, corresponding to 0.01 per cent transmission in the rejection spectral region.The presence of electronic components in the room (electrometers, PC monitors, and computer status LEDs), while conducting the full wavelength scan, may have impacted the beam characteristics, particularly because the IS was fully open at the beam entrance port.Furthermore, the insulation of the Xenon lamp was not perfect and some light leaked into the room.Careful shielding and light-blocking measures may be necessary to mitigate this effect.
Photodiode temporal and thermal stability: previous work by Results for modeling of the full scan PD data.From top to bottom : frequency  , amplitude , phase , offset , slope , slits widths  1 and  2 , and  2 /  .Bi-modality of  1 and  2 parameters originates from the chopper slits.The linear decreasing trend of the phase offset parameter  is thought to be a constant time lag between consecutive observations at high SNR.The settling time is somewhat constant with a small gap increase caused by the movement of the monochromator's grating turret performed to adjust the grating position that produces the desired wavelength.After reaching the fit lower boundary limit of -0.5, it reaches again a value close to the upper limit of 0.5.RASTI 000, 1-16 ( 2023) Design and performance of a portable CBP 11 The total statistic error is not shown as it overlaps with the SC uncertainty contribution curve.The gradual rises on the sides are due to lower light source power (near UV) and lower grating efficiency (near IR).Eppeldauer (1990) for the optical spectrum range (400-750 nm) showed that the spectral response changes of standard silicon photodiodes were as high as 0.5 per cent / year at constant operating temperature.Based on Larason & Houston (2008) findings, the temperature-dependent change becomes significant only for wavelengths exceeding 1000 nm, where it corresponds to 0.2 per cent/°C.Stubbs et al. (2010) are more conservative with an estimate of 1 per cent/°C above 900 nm and 0.1 per cent/°C for  < 900 nm.Continuous monitoring of the photodiode's temperature may be required to ensure reliable and accurate data acquisition.

Model error
The model is an accurate description of the dataset except for specific features in the residuals  =   −   that we call "model structured error" and is seen at each observation wavelength for the PD signal (see Fig. 12).We suppose these oscillations to be caused by the monochromatic light beam shape being blurred at the entry of the chopper wheel.If the shape has tails, the chopping mechanism may produce such small photocurrent oscillations and does not conform to a perfect sigmoid curve.The goal is to check if these periodic oscillations downgrade the fit quality if they impact the determination of the amplitude parameter   or add any bias.
For SC signals, the residual structure is below the noise level.
Two kinds of simulations are undertaken : (i) Measurement noise only.For the photodiode, we have direct access to a good estimate of the uncertainty   .We simulate a signal with the fit parameters   (  ).We draw random samples from a multivariate normal distribution using the model parameters covariance matrix obtained from the linear regression.Only the amplitude is kept at the original fit value.On these "perfectly" simulated signals, we add a random noise from a normal distribution N (0,   ).The signal is then binned at the multimeter resolution of 0.17 pA.For the solar cell, we simulate the correlated noise by contaminating a perfectly simulated signal with 50 Hz and 150 Hz lines and 1/f  noise directly in the time series.We similarly bin points at the multimeter resolution of 500 pA.
(ii) Measurement noise + modeling errors in the shape of the signal.Residuals between the model and PD observations are interpolated at each wavelength step   .Fig. 12 depicts the shape of these residuals.These are added to the same theoretically perfect simulated signals and standard measurement errors are added afterward.
At each wavelength and for the two simulated cases, we sample and fit afterward   = 10 4 SC and PD signals.The reconstructed parameters provide information on the effect of potential systematic effects introduced or omitted by the model.B1 and B2 give an example of 2- contours obtained at the dataset median wavelength  = 550 nm for both simulation cases.Overall, we find good agreement between model observations and simulations.No significant deviations between simulations and observations appear.We find a small constant bias of ∼0.003 pA on all the simulated PD signals when adding the structured residuals in simulation n°2.This is within the uncertainty of our   estimate and below the electrometer resolution, so the shift can be neglected.The average statistical uncertainty on   stands at ∼ 0.05 pA for the measurement noise simulation and ∼ 0.09 pA for the simulation with the added residuals.The model uncertainty contributions follow the expression represents the measurement uncertainty only (simulation n°1).The non-modeled residuals structure contributes to ∼ 83 per cent of the variance in the   estimate and the remaining 17 per cent originates from the intrinsic measurement noise.
However, the PD accounts for less than 0.05 per cent of relative uncertainty to the total CBP throughput uncertainty budget.The error budget is completely dominated by the SC signal noise.Thus, our model fitting observations without integrating the "model structured error" contribution is sufficiently good as long as the SNR on the SC signal does not drastically increase.Ratios of CBP throughput curves for the considered simulated cases are shown in Fig. 13.No wavelength dependency appears and the ratio of splines stays below the per-mil level.

DISCUSSION
As seen before, the current prototype of the traveling version of the CBP looks very promising.In order to be usable in the context of UV-VIS-NIR photometric survey instruments, some improvements on both the equipment and the acquisition method are required.We discuss their potential implementation hereafter.

Upgrading the light source to increase photon flux
Integration time is one of the key constraints to the design and performance of a CBP, together with the overall stability and precision.Indeed, telescopes operate for whole nights.That makes them available for calibration purposes during the day or more if the weather prevents observations.Moreover, integration times need to be compatible at the millimagnitude level for potential interference from daylight, given that telescope domes do not provide complete light insulation.Therefore, reaching the desired level of uncertainty on the instrument curve throughput determination requires integration times as short as possible.With a 175W Xenon light as we had for this preliminary feasibility test,   = 20 s is insufficient to reach the per-mil level of statistical uncertainty.For that reason, a more powerful light source must be found.Various options exist, such as laser-driven light sources (LDLS) which benefit from laser power without the large volume occupancy or potential hazard.They can also have a broad wavelength range (from UV to IR) with very high spectral radiance (tiny plasma diameters up to 1 mm producing radiance ranging from 1 to 100 mW/mm 2 /sr/nm, much higher than the radiance in our setup of ∼ 0.05 mW/mm 2 /sr/nm at peak), stability and lifetime.Furthermore, collimating light from traditional sources often results in a beam divergence or widening angle, which can be problematic.The emission surface is significantly higher for arclamps with millimeter-sized ones against ∼ 300 µm for the LDLS.The small emitting point of LDLS allows for efficient focusing of the light onto a very small area.Replacing the 175W Xenon light source with the IDIL TLS3-EQ77-UV-VIS-NIR monochromator and LDLS bundle is considered.The coupling between these two parts is achieved with optic fibers.Statistical uncertainty on each wavelength data point can be decreased by one or even two orders of magnitude if the effective radiance entering the IS is proportionally higher.Nevertheless, it is important to stay in tune with the upper flux limit imposed by the telescope to calibrate for the real-use case scenario of the traveling CBP instrument.This issue arises due to the potential saturation of the CCD/CMOS camera, resulting from the injection of an excessive amount of monochromatic flux into the sensor within a brief exposure time.

Speeding up the chopping frequency
One caveat observed in post-processing is that the signal frequency set at 6 Hz is too slow to avoid contamination by the 1/f noise.To quantify the potential gain expected from running at higher chopping frequency, we run simulations using our models.Knowing the fitted parameters of each signal for both SC and PD at each wavelength of the acquisition, we increase the frequency to  2 = 12 Hz and  3  = 18 Hz to evaluate the impact on the estimated amplitudes.Signalto-noise ratio improve by ∼ 25 per cent when going from 6 to 12 Hz and ∼ 50 per cent from 6 to 18 Hz.Above 18 Hz, the overall gain of SNR becomes smaller because the signal peak is already far from the 1/f curve.Thus, rotating the chopper wheel motor at twice or three times the initial frequency or designing a mask with two additional symmetric holes will lead to a noticeable improvement of the SNR on the CBP throughput curve.For the chopping system, off-the-shelf options are plentiful, or even careful redesign to improve our custom DIY device is feasible.

Monitoring sensors temperatures
Accurate and reliable measurements with a CBP require careful monitoring of the temperature of critical components such as the photodiode and solar cell.The responsivity of these components is known to be sensitive to temperature variations, especially in the infrared  > 1000 nm where responsitivity deviations can reach 1 per cent/°C for silicon-based photodiodes.

Double monochromator configuration
Achieving both the desired spectral radiance and bandwidth is difficult with a monochromator and a broadband light source as demonstrated in Stubbs et al. (2010).In the past, a tunable laser was used to overcome this issue for previous CBP versions.In this paper, we have seen that target accuracies are within reach despite using an inappropriate monochromator.Our monochromator model was optimized for particle physics experiments (i.e.Cerenkov effect measurements, explaining the 250 nm blaze wavelength of the unique grating installed).It suffered from poor efficiency in the exploited range of observations averaging at 40 per cent for  = 550 nm (see Fig. 3).
The etendue is also low as the grating width   = 30 mm.As it is squared in the etendue equation, we would benefit tremendously from a grating twice as large (factor 4 of improvement in terms of flux).Upgrading the monochromator with appropriate larger gratings (UV + VIS + IR) could drastically increase the flux entering the IS, the SNR and reduce the integration time all at once.For example, gratings up to 68 mm widths are available for purchase with > 70 per cent efficiency at peak.With this kind of configuration, an increase of around one order of magnitude on the signal-to-noise ratio can be expected at every wavelength sample.Nonetheless, two potential additional systematic errors to address arise from the spectral bandwidth of the beam.While keeping the necessary amount of flux inside the IS to know the system transmission at better than 0.1 per cent level, it is also important to get beam spectra as narrow as possible in order to adequately capture the steep edges of filters (∼ 10 nm for ZTF filters) with a high spectral resolution.The recommendation stands to have optical bandwidths of ∼ 1 nm.The second systematic effect to consider involves the variation in the spectral bandwidth shape with wavelength.In Fig. A1, we observe the wavelength dependency of the FWHM of the monochromatic beam generated by the monochromator.The FWHM decreases from ∼ 3.5 nm at  = 350 nm to ∼ 1.8 nm at  = 750 nm.This change impacts the measurement of the slope of filter edges since the observed CBP telescope response function is the convolution of the monochromatic filter response with the spectral bandwidth of the monochromator's beam.The width evolution leads to systematic and significant differences at the red and blue edges of a filter, causing a tendency to shift the observed filter curve in one direction or another.Adding a second monochromator can tackle multiple sources of systematic uncertainties at once by reducing ambient light contamination (with a fully fibered optical path), enhancing spectral purity (light beam is diffracted one additional time), reducing out-of-band emission, and decreasing the beam spectral bandwidth.Loss of light flux is to be anticipated but may be compensated by carefully injecting light into the entrance slit of the first monochromator, and with twice as large and efficient gratings.

Effect of individual sample uncertainty and wavelength scan step on zero-point and central wavelength precision
In planning the next iteration of the traveling CBP instrument, we present forecasts for statistical uncertainties in zero-point magnitude    and central wavelength   metrics.This assessment serves two purposes: (i) to determine the necessary level of precision for individual observation samples, and (ii) to establish the required wavelength scan step (Δ in nm).Fig. 14 depicts the anticipated uncertainty levels in zero-point and central wavelength for the ztf-g band.We proceed in the same way as in Sections 4.2 and 4.3 but with a single relative uncertainty value for all samples.This relies on the assumption of a flat relative uncertainty distribution across the wavelength range.Then, we evaluate the metrics    and   with Equations 10 and 11 for different wavelength steps Δ.Maintaining a constant wavelength scan step of Δ = 2 nm, as chosen in this study, will allow reaching the mmag threshold in zero-point accuracy and a 0.1 nm threshold in central wavelength accuracy, provided that the uncertainty in individual observation samples for the CBP telescope throughput curve is reduced to approximately 0.4 per cent.This corresponds to a factor 4 to 6 improvement considering ztf-g and ztf-r bands respectively, from the current prototype setup.Consequently, it is important to give careful consideration to the balance between the overall count of observation samples (which directly affects the total integration time for the CBP throughput scan) and the SNR of each observation sample (requiring applicable solutions to mitigate the SC noise).We emphasize that future work should address this trade-off and explore solutions to minimize exposure time (preventing sensor sensitivity variations due to changes in environmental conditions over large periods of time) and maximize SNR in the SC.

CONCLUSIONS
We assessed the replacement of a laser with a broadband light source to design a traveling Collimated Beam Projector.The critical aspect of designing such a calibration device is the measurement accuracy of the intermediary telescope acting as a light projector.By calibrating the CBP telescope relative to a solar cell itself calibrated at better than 0.1 per cent onto the NIST reference photodiode, we can propagate the calibration to telescope systems.The proof-of-concept setup presented in this paper has shown promising results for its intended use.It delivers average statistical uncertainty below 2 per cent in the 350 to 750 nm range which propagates into zero-point magnitudes uncertainties of ∼ 5 mmag for ztf-g and ztf-r bandpasses.Filters central wavelengths would be determined at the sub-nanometer precision, considering statistical uncertainties only.In its current state, the accuracy of the throughput curve is limited by the solar cell intrinsic dark noise accounting for nearly the entire statistical uncertainty budget.The SC still remains the preferred intermediary calibration sensor because it possesses a large photosensitive surface area adequate for capturing the entirety of photons emitted through the CBP telescope.To alleviate it, multiple efforts of methods and equipment improvements must be undertaken : (i) increasing the photon flux in the IS by using wider and more efficient gratings, upgrading the light source and refining the light injection mechanism; (ii) chopping the light beam at two or three times the current frequency to diminish the contamination of 1/f noise on the SC measurements; (iii) multiple scans as statistical uncertainty is expected to decrease by 1/

√
with  number of measurements; (iv) monitoring sensor temperatures; (v) enhance baffling to reduce ambient/stray light contamination; (vi) install a spectrograph with standard emission lines calibration lamp.Simultaneously, some questions regarding the stability of illumination need to be looked closer in the future.Repeatability of the throughput measurement must be tested by performing repeated scans at specific wavelengths and checking the reliability of the results.Systematic uncertainties can be kept below the statistical uncertainty threshold if concrete steps are taken toward the rigorous mitigation and assessment of each contribution mentioned in the paper.Additional investigation of the impact of beam spectral width on the filter's steep edges is also required.
Overall, the proposed system with detailed improvements has the potential to reach the systematic uncertainty limit of the CBP throughput curve measurement.Future work will focus on building the final refined system, adding the spectral quality evaluation in the analysis, and testing its "traveling" capabilities by calibrating the StarDICE experiment telescope at Observatoire de Haute-Provence.Having a calibration system that is transportable across different observatories is crucial for SNe Ia cosmology in the era of wide-field surveys collecting large amounts of data, where photometric calibration continues to be a significant limiting factor in enhancing cosmological constraints.This paper is dedicated to the memory of our colleague Éric Nuss, who passed away in April 2023.5 .9 9 7 2 5 .9 9 7 6 5 .9 9 8 0 5 .9 9 8 4

Figure 4 .
Figure 4. Top panel: Sketch of the chopper wheel.The mask has 2 blades (black), rotates at 3 Hz, and lets the beam travel through the aperture (colors).Green: fully open aperture state, orange: partially open aperture state; red: fully closed aperture state.Bottom panel: sketch of a time series resulting from the rotation of the chopper wheel; the colors correspond to the situations depicted in the top figure.

Figure 5 .Figure 6 .
Figure5.Semi-log scale illustration of the PD signal (red) superposed with different noise contributions.The orange area represents the total readout noise ±  obtained from the root-mean-square (RMS) of 4 independent dark signals.The green area is the shot noise ± ℎ .The blue line is the dark current    .The chopped aspect of the signal is visible with a period of 1/6 s set by twice the rotating mask frequency of 3 Hz.The dashed black line represents the total noise   .It should be noted that for the sake of visualization, we only display the maximum amplitude of the time-dependent noise time-series in the plot as horizontal lines.

2Figure 7 .
Figure 7. Top panel: PD raw observations are plotted as dotted black points with ±1- uncertainty error bars and the red curve show the fitted model.Time series are folded between a single period.Best-fitting parameters with ±1- uncertainties are shown on the upper right.Bottom panel: residuals between observations and best-fitting model.They appear structured, showing some inadequacy of the model to catch small fluctuations of photocurrent caused by the blurred shape of the light beam.

Figure 8 .
Figure 8. Top panel: SC raw observations are plotted as black crosses and the blue curve depicts the fitted model.Best-fitting parameters with ±1- uncertainties are shown on the upper right.Bottom panel: residuals between observations and best-fitting model.The presence of patterns in the residuals, observed during transition phases (when the photocurrent goes up and down at points 0.25, 0.75, 1.25, and 1.75), is believed to be caused by oscillations in a 50 Hz signal.
Figure9.Results for modeling of the full scan PD data.From top to bottom : frequency  , amplitude , phase , offset , slope , slits widths  1 and  2 , and  2 /  .Bi-modality of  1 and  2 parameters originates from the chopper slits.The linear decreasing trend of the phase offset parameter  is thought to be a constant time lag between consecutive observations at high SNR.The settling time is somewhat constant with a small gap increase caused by the movement of the monochromator's grating turret performed to adjust the grating position that produces the desired wavelength.After reaching the fit lower boundary limit of -0.5, it reaches again a value close to the upper limit of 0.5.

Figure 12 .
Figure 12.Residuals between observations and best-fitting model for the PD signal for all wavelengths (with coherent matching colors).The asymmetric pattern appears strongly defined.The mean curve of residuals is plotted in white.

Figure 13 .
Figure13.Ratio of CBP throughput curves for fitting PD signal with and without residuals (see Fig.12).Individual samples are plotted in black crosses and the spline fit ratio is drawn in orange curve.No bias is visible.Deviations at low and high wavelengths are due to SNR decreasing.

Figure 14 .
Figure 14.Top panel: forecasts of zero-point magnitude uncertainty    as a function of individual observation sample relative uncertainty    ()/   () [per cent] and scan step Δ [nm] for the ztf-g band.Bottom panel: same for central wavelength uncertainty  () forecasts.Similar results are found for the ztf-r filter.

Figure A1 .
Figure A1.Spectrograph data with Gaussian fitting for 5 representative wavelengths.Δ shows the level of agreement between monochromator and spectrograph, while the FWHM gives the spectrograph's resolution.

Figure B1 .
Figure B1.Posterior distributions for the fit of two mocks realizations of the PD simulated signal at  = 550 nm without (blue) and with (orange) residuals from the original observation and best-fitting curve.The colored contours show 1 and 2- joint confidence intervals, the dotted black lines are the true values.The "True" column of the table presents the original best-fitting results on the observation.

Figure B2 .
Figure B2.Posterior distributions for the fit of a mock realization of the SC simulated signal at  = 550 nm.The blue contours show 1 and 2- joint confidence intervals, the dotted black lines are the true values.The "True" column of the table presents the original best-fitting results on the observation.

Table 1 .
SC dark signal best-fit values with their associated uncertainties and relative error.N/A denotes parameters directly extracted from experimental data using the tools from scipy.signal.