HFI Detector Spectral Response

From Planck PLA 2015 Wiki
Revision as of 16:55, 4 October 2012 by Lspencer (talk | contribs)
Jump to: navigation, search

HFI Detector Spectral Response[edit]

Abstract[edit]

This section is comprised of excerpts from the HFI - Spectral Calibration report of the IAS measurements, v3.00, an internal report describing detailed pre-flight ground based measurements of the HFI focal plane. The experimental setup, data collection, and related data processing are described. The official version of the HFI detector spectral transmission profiles is available within the HFI instrument model and the RIMO files in the Planck Legacy archive [ref.] . This data is comprised of broadband Fourier transform spectrometer (FTS) measurements conducted with the HFI focal plane assembly as described above, and includes a waveguide model for the low frequency spectral region, and component-level filter spectra for the remaining out of band spectral regions. Specific attention is given to in-band and near-band spectral regions surrounding CO rotational transitions in order to support the CO extraction component separation effort [ref.] . The spectral transmission profiles are evaluated with parameters such as cut-on, cut-off, centre frequency, effective frequency (including spectral index), and band-width are provided in the analysis. Further evaluation yields band-average spectra ( [ref.] ) and unit conversion / colour correction coefficients ( [ref.] ) and software routines ( [ref.] ).

Objectives[edit]

The main goal of the spectral transmission tests of the HFI instrument is to measure the spectral response of all HFI detectors to a known source of EM radiation individually. This was determined by measuring the interferometric output of all detection channels for radiation propagated through a continuously scanned polarising fourier transform spectrometer (FTS). The required accuracy to which the spectral transmission is to be recovered is 1[math]\%[/math] (see §Req.).

It is important to note that the absolute spectral calibration cannot be achieved solely from the analysis of the FTS data because of uncertainties in the coupling efficiency of the FTS source through the FTS, input optics, and integrating sphere. The relative FTS measurements must be combined with the optical efficiency tests which used internal black-body sources (EFF Test – see §EFF).

Description of setup[edit]

A brief description of the optical setup used for the spectral characterisation of the HFI instrument performed at the FPU level in the Saturne calibration facility at IAS is provided here. Further details are available from the pre-launch calibration papers . The spectrometer configuration is based on a polarizing Fourier transform interferometer, feeding an integrating sphere inside the Saturne tank (Figs. [fig:Optical-layout-of] & and [fig:Top-view-of]). The output of the integrating sphere is conjugated by a mirror with the focal plane of the HFI. Therefore data are acquired simultaneously on all pixels while the moveable mirror in the external FTS is scanned.

Internal optical layout of the FTS.
FTS/ELS assembly with Saturne cryostat.


Optical layout of the FTS with respect to the Saturne cryostat.


Top view of 2K optical setup inside Saturne cryostat.

An external light source system (ELS, see Fig.[fig:FTS/ELS-assembly-with]) is provided in addition to the FTS for time constant measurements. Selection between the FTS and the ELS is done by the M1 flip mirror.

The polarising FTS was provided by Sciencetech Inc. (Canada) under contract with IAS and ESA (Fig.[fig:Internal-optical-layout]). The source can be selected from either a Hg vapour arc lamp or a globar. The 300 mm translation stage allows symmetric interferograms with a maximum theoretical (unapodized) Instrument Line Shape (ILS) Full-Width-at-Half-Maximum (FWHM) resolution of [math]0.020~\mbox{cm}^{-1}[/math]. Filtering of short wavelength is done on the optical path entering the Saturne cryostat as follows:

  • Saturne vacuum window is polyethylene (6mm thick)
  • [math]1^{st}[/math] thermal fiter at 300K
  • [math]2^{nd}[/math] thermal filter at 77K
  • [math]46~\mbox{cm}^{-1}[/math]filter at 77K
  • [math]3^{rd}[/math] thermal filter at 20K
  • filter wheel at 2K (Fig. [fig:2K-filter-wheel]) with selection of: open, closed, [math]35~\mbox{cm}^{-1}[/math] or [math]10~\mbox{cm}^{-1}[/math]filter

FIXME: typo in figure

2K filter wheel on FTS optical path.

FIXME: typo in figure

Saturne integrating sphere and HFI focal plane.

The optical flux from the FTS enters the integrating sphere through an 8mm diameter hole at its bottom (Fig. [fig:FTS-input-aperture]), which is the conjugate of the M1 mirror by M2.

During scanning, encoder pulses from the translation stage are time stamped with a clock synchronized with the HFI signal acquisition clock (itself provided by the spacecraft simulator used during the calibration). In a similar way, a mechanical ZPD (Zero Path Difference) signal is stamped and stored in the database. After verification, the scanning speed was found to be sufficiently stable to be assumed as constant along the section of the interferogram needed for spectral processing. Therefore the data processing described below does not use the reference signal but relies instead on the data acquisition time line to derive the optical path difference.

To provide a reference monitor for the interferometric flux, a dedicated bolometer is used within the integrating sphere. The flux recorded by this detector is directly related to the flux received by the HFI detectors on the focal plane which is coupled directly to the sphere output hole via mirror M1. The reference bolometer was provided by the team of N. Coron at IAS and has a predetermined absolute calibration between the flux measured and its output signal. It is fed by a modified Winston horn from Infrared Lab. Inc. (USA) and operated at 300 mK using a [math]^{3}He[/math] fridge (Torre-Chanin). The reference spectra acquired during the calibration run with this bolometer are used to identify standing wave features in the FTS source and path (lamp, windows,...) and to check the shape of the source spectrum. The diffraction losses at low frequencies due to the horn exit aperture diameter (2.6 mm ) is a limitation taken into account in the data processing of the low frequency channels of HFI.

Data Collection[edit]

To allow for data sets that have the same instrumental effects the bulk of the sequences have been performed with the same parameters. Also, separate single tests have been performed in different configurations to understand instrument systematics. The following is a list of the sequences adopted for the final spectral calibration data.

  • Max OPD: The Maximum Optical Path Difference adopted during FTS testing is 30 cm. This relates to a nominal maximum resolution for the spectra of [math]\simeq 0.02~\mbox{cm}^{-1}[/math]. This maximum resolution may then be deteriorated by the actual sample of data selected and by the apodization of the data (see §[sec:apod]).
  • Speed: The Speed of the scanning FTS is set to [math]0.1~\mbox{cm}/\mbox{s}[/math] for all of the data. The sampling of the data given by the acquisition frequency allows knowledge of the spectrum up to [math] \simeq 430~\mbox{cm}^{-1}[/math]. This is well above the optical filter cut-off frequency, so aliasing is not a concern. The independent filter measurements used in the spectral transmission outside of the optical pass-band (see §[sec:filter]) have a Nyquist frequency of [math]600~[/math]cm[math]^{-1}[/math].
  • Source: two sources were used in separate measurements. The mercury arc lamp as the main source for the final data and also a GLOBAR[math]^{\circledR}[/math] source to study any systematic effect present in the different measurements. (Although ratio with reference bolometer should remove such artifacts.) The spectra used to generate the Spectral Transmission data products in the IMO all originated from the mercury arc lamp as the optical source.
  • Spectral range: The spectral range of the radiation entering the Saturne cryostat is given by a fixed [math]46~\mbox{cm}^{-1}[/math] LPE filter and a rotating filter wheel with the following positions (1- Open (no added filtering), 2- Dark (no radiation allowed), 3- LPE [math]=10~\mbox{cm}^{-1}[/math], 4- LPE [math]=36~\mbox{cm}^{-1}[/math]). The last two positions have been adopted for the final data of the three lower spectral channels (100, 143, 217 GHz) and the higher three spectral channels (353, 545, 857 GHz), respectively.

Additional Experiments[edit]

There were two significant additional tests that are used in the derivation of the HFI detector spectral transmission profiles beyond the scope of the IAS HFI FTS measurements. Additional filter measurements were recorded in the AIG test facility at Cardiff during filter stack production. These measurements are used to extend the IAS FTS spectral meaurements beyond the HFI spectral passband. The EFF tests are used to obtain optical efficiency parameters for each detector. These parameters, when combined with the respective normalized spectral transmission profiles, provide an estimate of the absolute spectral transmission. The EFF tests are discussed briefly in §EFF, with greater detail provided by A. Catalano . The filter measurements are described in §filter below.

Filter Measurements[edit]

Prior to the IAS measurements with the integrated HFI detectors and filter stacks, FTS measurements of the individual filters comprising the filter stacks for each band were conducted at Cardiff. As will be discussed in §[sec:OOB], the independent measure of the filter stack transmission is used for a portion of the HFI detector spectral transmission for regions of the spectrum where it is deemed to be of better quality then the IAS FTS measurements (i.e. for frequencies outside of the band edge filter cut-off(s)). The filter stacks for each of the frequency bands are comprised of 5 filters. There is an additional low frequency cut-on filter for the 545 and 857 GHz bands as the waveguide cut-on is too low. Figures [fig:Filt100] – [fig:Filt857] show the individual filter transmission measurements as well as the combined filter stack product.

The uncertainty on the combined filter transmission measurement is determined as follows. Let [math]f_i(\sigma)[/math] represent the individual filter transmission. The standard deviation of all values of [math]f_i(\sigma)[/math] below a threshold of [math]0.001[/math] is used as an approximate uncertainty for each individual filter measurement, i.e. [math]\varsigma_{i}[/math]. The individual uncertainty estimates are combined to provide an estimate of the combined filter transmission spectral uncertainty through standard error propagation. For the combined filter transmission represented as [math]F(\sigma)[/math], the associated uncertainty estimate is given as

[math]\label{eq:Ferror} \varsigma_F(\sigma) = \displaystyle\sqrt{[F(\sigma)]^2\Sigma_i [(\frac{\varsigma_i}{f_i(\sigma)})^2] } ~.[/math]

Raw Data[edit]

Table 1 contains the details for the sets of data adopted for data analysis in the generation of the HFI detector spectral transmission profiles. The first two columns report the date and time of the starting point of the data set. The third column indicates the approximate scan length (in minutes). The fourth column is the maximum number of valid interferograms that are present in the data file. The column labeled RFP (Rotating Filter Wheel Position) indicates the filtering condition under which the data was taken: 1-Totally open (no filtering other than the ever present [math]46~\mbox{cm}^{-1}[/math] LPE; 2-Closed (no radiation is allowed to pass into the cryostat from the FTS); 3-[math]10~\mbox{cm}^{-1}[/math] LPE filter (for the low frequency channels); 4-[math]36~\mbox{cm}^{-1}[/math] LPE filter (for the high frequency channels). The column immediately following indicates, within a milliKelvin, the average reading of the bolometer plate thermometers during the scan time. The data sets with the RFP set to [math]3[/math] were used for the lower frequency channels, i.e. the 100, 143, and 217 GHz bands. The data corresponding to an RFP of [math]4[/math] were used for the higher frequency channels, i.e. the 353, 545, and 857 GHz bands. Table [tab:BrefData] indicates the data sets used to obtain the reference bolometer spectra. In most instances the data were acquired synchronously, however, there were cases when the HFI data was of good quality and the reference bolometer data was not, and vice-versa. The recorded data corresponding to the entries in Tables [tab:RawData] & [tab:BrefData] were used in the determination of the HFI spectral transmission profiles. There is also additional data that was recorded as part of the test campaign but not used in this specific work that is not listed here.


Table n:
Summary of recorded data used to generate the spectral transmission profiles.

Date (DD/MM/YY) Time length (min) #Ifgm. RFP Source T(mK)
25/06/06 07:52 25 4 4 Hg 140
25/06/06 08:33 25 4 4 Hg 140
25/06/06 09:02 65 12 4 Hg 140
25/06/06 10:05 65 12 4 Hg 140
25/06/06 11:15 65 12 4 Hg 140
26/06/06 09:44 25 4 4 Hg 140
27/06/06 17:48 148 28 3 Hg <130
27/06/06 20:20 335 50 3 Hg <120
28/06/06 21:39 163 28 3 Hg 95
05/07/06 06:23 82 16 4 Hg 100
Summary of recorded data used to generate the reference bolometer spectra.
Date (DD/MM/YY) Time length (min) #Ifgm. RFP Source
22/06/06 17:50 50 9 3 Hg
22/06/06 18:40 50 9 3 Hg
22/06/06 19:30 50 9 3 Hg
22/06/06 20:20 50 9 3 Hg
22/06/06 21:10 50 9 3 Hg
22/06/06 22:50 50 9 3 Hg
22/06/06 23:40 50 9 3 Hg
25/06/06 07:52 25 4 4 Hg
25/06/06 10:05 65 12 4 Hg
25/06/06 11:15 65 12 4 Hg
26/06/06 09:44 25 4 4 Hg
26/06/06 17:04 38 7 3 Hg
26/06/06 22:23 41 8 3 Hg
27/06/06 00:04 86 16 3 Hg
27/06/06 17:48 148 28 3 Hg
27/06/06 20:20 335 50 3 Hg
28/06/06 21:39 163 28 3 Hg
05/07/06 06:23 82 16 4 Hg

Data Processing[edit]

<1-- [sec:DP] -->

Having identified the data sections that are of interest, here follows the data processing sequence used to obtain the resulting spectra. The processing steps taken are as follows:

  1. Selection and extraction of time sampled data sets
  2. Conversion of time sampled data to arrays of OPD sampled interferogram data sets
  3. Fourier transformation and averaging of interferogram data sets
  4. Discrimination of poor quality spectra by standard deviation comparison (see §[sec:FTavg])
  5. Division of detector spectra by reference bolometer spectra to obtain normalized spectral transmission profiles
  6. Combination of relative transmission spectra with filter measurements and waveguide models
  7. Determination of optical efficiency through evaluation of normalized spectral transmission in the context of the EFF experiments
  8. Addition of over-sampled data into spectrum for the CO transition regions (see §[sec:CO])
  9. Identification of common frequency sampling per channel and interpolation of spectra onto the common sampling

Selection and extraction of data sets[edit]

The selected arrays of data are imported into a KST client from the QLA server

plck-op-3.ias.u-psud.fr

and subsequently exported as text files. Separate files are generated for the various bolometer signals and sample times for each detector. The exported text files were transferred from the remote IAS machines to the dedicated Planck computer at Cardiff,

planck01.astro.cf.ac.uk, 

for further processing and analysis. A sample bolometer signal timeline is shown in Fig. [fig:BolSig].


The bolometer signal is stored within the initial database in several formats: raw ADU, signal voltage, resistance, current, temperature, total power, electrical power, and radiant power, all of which may be exported to a KST session as a function of sample time. Although all of the above formats originated from the recorded detector ADU signal and instrument model parameters, it was not possible to trace these refined database vectors to the specific bolometer parameters used in their conversion. Also exported from a KST session were the set balance voltages and bias currents for each detector so that a bolometer model using the current IMO parameters (IMO ver. 2_50) could be applied to the raw data as well. The motivation for the application of a bolometer model to the data was to reduce the in-band effects of detector nonlinearities as evidence of detector operation in a non-linear regime was identified within the raw data.

For the interferograms within each of the data sets in the above tables (Tab. [tab:RawData] & [tab:BrefData]), a course estimate of the temporal location of each position of Zero Optical-Path Difference (ZPD) was determined using the characteristic local maxima associated with a broad-band fringe for detectors with good S/N. This course ZPD location was then used as a starting point to identify a unique ZPD location within each interferogram through the fitting of a low order polynomial to the ZPD fringe.

The interferogram boundaries are determined as the mid-points between subsequent interferograms, less a small number of buffer points ([math]\sim[/math]10) to ensure that the extracted interferograms include regions associated with the FTS stage travel having constant velocity while excluding the acceleration regions.

Visual verification of the extracted interferograms is performed to ensure that each ZPD was identified correctly and to remove any low-quality interferograms. The overlap of the extracted interferograms is also verified visually, with any anomolies either corrected (through a revised ZPD location) or removed. These interferograms are then saved to separate files, to be combined as a large dataset further on in the processing.

Once interferograms have been extracted from each data set, they are combined into a single array for each detector. The Maximum Optical Path Difference (MPD) is determined for each interferogram, and the minimum for each detector within a combined data set is used to generate an evenly sampled Optical Path Difference (OPD) grid onto which each interferogram in the combined data is then interpolated. This ensures that each individual interferogram is sampled at ZPD and that each spectrum has identical frequency sampling and can thus be averaged together. An example of a combined interferogram data set is illustrated in Figure [fig:CombinedIFGMex] where the central portions of the recorded interferograms are shown. A similar plot for each HFI detector is included in the Appendices (see §[sec:allIFGMs]). The MPD value for each detector, and corresponding spectral resolution, is shown in Fig. [fig:MPD].


/figures/LSfigs/

Fourier Transformation and Averaging of Interferogram Data Sets[edit]

In preparation for Fourier transformation, a low-order polynomial baseline removal is performed on the individual interferograms. Consequently, no information can be recovered from the spectrum below [math]\sim 0.1~\mbox{cm}^{-1}[/math], but this is of no concern as this region of the spectrum is replaced by a waveguide fit in the final data product (see §[sec:WG]). An average interferogram is determined and used to identify glitches for removal from the interferogram data. If possible, the exponentially decaying tail of a glitch is identified, fit, removed, and replaced by the average of the remainder of the data set (which does not contain a glitch). If this is not possible then the entire interferogram is removed from the data set.

Although an average interferogram is determined, it is not directly used in obtaining the average spectrum. Once the data set has been deglitched, individual interferograms may be apodized (see §[sec:apod]), and then undergo Fourier transformation and phase correction . To allow phase correction to take place at the individual interferogram level, data averaging is then performed in the spectral domain. The uncertainty for every spectral data point is determined statistically through the standard deviation at a given frequency. A check for poor quality spectra is performed by comparing the overall standard deviation including and excluding any given spectrum. An example of the individual spectra and uncertainty for bc00 is shown in Fig. [fig:CombinedSpecEx]; similar plots for all of the detectors are shown in §[sec:allSpecs].

Table [tab:Igstats] contains information on the individual detector data sets for each of the HFI detectors. Columns three and four contain the number of interferograms and spectra, respectively, used in the generation of subsequent data. These numbers may disagree due to the discrimination of poor quality data discussed above. There will never be more spectra than interferograms for a given data set. The fifth column indicates the minimum MPD for the detector data set. The FTS Instrument Line Shape (ILS) Full-Width-at-Half-Maximum (FWHM) for the data set is provided in the sixth column. An estimate of the S/N is obtained using the average spectrum and its statistical uncertainty, averaged across the in-band region of the spectrum.

The averaged spectra are then normalized and divided by a (normalized) reference bolometer spectrum (see §[sec:setup] & §[sec:Bref]). This is done to remove spectral features from the data beyond those due to the HFI optics and detectors. The terminology ‘pre-ratio’ is used to refer to the spectra before division by the reference bolometer spectra.

|c|c|c|c|c|c|c|



<tbody> </tbody>
Band
(GHz)

& bc & # Ifgm. & # Spec. & MPD (cm) & ILS[math]_{\mbox{FWHM}}[/math] (cm[math]^{-1}[/math]) & avg. S/N


<tbody> </tbody>
Band
(GHz)

& bc & # Ifgm. & # Spec. & MPD (cm) & ILS[math]_{\mbox{FWHM}}[/math] (cm[math]^{-1}[/math]) & avg. S/N
100 & 00 & 96 & 93 & 29.648639 & 0.020355066 & 104.70555
100 & 01 & 96 & 95 & 29.648639 & 0.020355066 & 203.34226
100 & 20 & 96 & 96 & 29.649801 & 0.020354268 & 271.38262
100 & 21 & 95 & 94 & 29.650962 & 0.020353471 & 264.65385
100 & 40 & 96 & 94 & 29.652124 & 0.020352674 & 432.40921
100 & 41 & 96 & 96 & 29.652124 & 0.020352674 & 962.98064
100 & 80 & 95 & 95 & 29.650962 & 0.020353471 & 262.71036
100 & 81 & 96 & 96 & 29.650962 & 0.020353471 & 227.35782
143 & 02 & 96 & 95 & 29.652124 & 0.020352674 & 434.93136
143 & 03 & 96 & 96 & 29.652124 & 0.020352674 & 458.86741
143 & 10 & 95 & 95 & 29.652124 & 0.020352674 & 565.21311
143 & 30 & 92 & 92 & 29.653286 & 0.020351876 & 424.53824
143 & 31 & 92 & 92 & 29.652124 & 0.020352674 & 447.50343
143 & 42 & 92 & 92 & 29.652124 & 0.020352674 & 536.23423
143 & 50 & 96 & 96 & 29.652124 & 0.020352674 & 441.68983
143 & 51 & 96 & 96 & 29.652124 & 0.020352674 & 476.22684
143 & 60 & 96 & 96 & 29.652124 & 0.020352674 & 517.43269
143 & 70 & 95 & 94 & 29.652124 & 0.020352674 & 326.06542
143 & 82 & 96 & 96 & 29.652124 & 0.020352674 & 454.23293
143 & 83 & 96 & 96 & 29.652124 & 0.020352674 & 439.50612
217 & 04 & 96 & 96 & 29.652124 & 0.020352674 & 625.92713
217 & 11 & 96 & 96 & 29.652124 & 0.020352674 & 580.32904
217 & 12 & 67 & 67 & 29.663741 & 0.020344703 & 549.24264
217 & 22 & 96 & 96 & 29.652124 & 0.020352674 & 550.65242
217 & 43 & 95 & 95 & 29.652124 & 0.020352674 & 616.57750
217 & 44 & 96 & 96 & 29.652124 & 0.020352674 & 543.78288
217 & 52 & 95 & 94 & 29.652124 & 0.020352674 & 684.91039
217 & 61 & 92 & 92 & 29.652124 & 0.020352674 & 594.70296
217 & 62 & 95 & 95 & 29.652124 & 0.020352674 & 579.69385
217 & 71 & 96 & 96 & 29.652124 & 0.020352674 & 631.30065
217 & 72 & 96 & 96 & 29.652124 & 0.020352674 & 687.10844
217 & 84 & 92 & 92 & 29.652124 & 0.020352674 & 636.36790
353 & 05 & 55 & 55 & 29.651916 & 0.020352816 & 354.15815
353 & 13 & 57 & 57 & 29.651916 & 0.020352816 & 319.09858
353 & 23 & 55 & 55 & 29.651916 & 0.020352816 & 218.05772
353 & 24 & 64 & 64 & 29.651916 & 0.020352816 & 236.44464
353 & 32 & 56 & 56 & 29.651916 & 0.020352816 & 238.40018
353 & 33 & 55 & 55 & 29.651916 & 0.020352816 & 219.60443
353 & 45 & 64 & 64 & 29.653078 & 0.020352019 & 246.97889
353 & 53 & 57 & 57 & 29.651916 & 0.020352816 & 225.27621
353 & 54 & 57 & 57 & 29.651916 & 0.020352816 & 228.47058
353 & 63 & 57 & 56 & 29.651916 & 0.020352816 & 140.16991
353 & 64 & 55 & 55 & 29.651916 & 0.020352816 & 158.85037
353 & 85 & 57 & 57 & 29.653078 & 0.020352019 & 246.18460
545 & 14 & 63 & 63 & 29.654240 & 0.020351221 & 448.67624
545 & 34 & 64 & 64 & 29.654240 & 0.020351221 & 481.87111
545 & 55 & 57 & 57 & 29.654240 & 0.020351221 & 328.07835
545 & 73 & 57 & 57 & 29.654240 & 0.020351221 & 447.49901
857 & 25 & 64 & 64 & 29.653078 & 0.020352019 & 501.25106
857 & 35 & 64 & 64 & 29.654240 & 0.020351221 & 567.55125
857 & 65 & 42 & 42 & 29.653078 & 0.020352019 & 482.41906
857 & 74 & 64 & 64 & 29.654240 & 0.020351221 & 491.89133


Apodization[edit]

Apodization of FTS interferograms results in reduced sidelobes of the FTS spectral ILS, at the cost of broadening the central feature of the ILS itself. Another perspective is that the interferogram is prefferentially weighted at certain OPD regions (i.e. ZPD) over others prior to Fourier transformation. This also, ideally, results in an improvement in the spectral S/N for a broad-spectrum source, again at the cost of reduced spectral resolution. Figure [fig:apodKernel] illustrates a variety of interferogram apodization kernels that have been used in this work . The modified Norton-Beer 1.5 apodization function has been selected to be used for the final spectral transmission profile data set as it represents a good compromise between the desired ILS sidelobe reduction and improved S/N with marginal reduction in spectral resolution. Figure [fig:ApodCompare] illustrates the difference between apodizations on the same sample data.




Reference Bolometer Spectra[edit]

The reference bolometer spectra are obtained in a fashion similar to that used for the HFI detectors. Wherever possible, the same data processing is applied to the reference bolometer data as was applied to the HFI detector data, including RFP setting, scan speed, scan length, source intensity, apodization, phase correction, etc. Table [tab:Brefstats] contains the reference bolometer data set properties corresponding those listed for the HFI detectors in Table [tab:Igstats]. Figures [fig:BrefL] - [fig:Bref857] illustrate the resultant spectra and S/N from the reference bolometer data sets. Figure [fig:SNall] compares the approximate S/N of the average spectrum for each detector against the reference bolometer average spectrum S/N over the same spectral region.


|c|c|c|c|c|c|c|

<tbody> </tbody>
Band
(GHz)

& RFP & # Ifgm. & # Spec. & MPD (cm) & ILS[math]_{\mbox{FWHM}}[/math] (cm[math]^{-1}[/math]) & avg. S/N
100 & 3 & 164 & 164 & 29.644698 & 0.020357772 & 11.478057
143 & 3 & 164 & 164 & 29.644698 & 0.020357772 & 37.361317
217 & 3 & 164 & 164 & 29.644698 & 0.020357772 & 126.14280
353 & 4 & 24 & 24 & 29.652306 & 0.020352549 & 136.06159
545 & 4 & 24 & 24 & 29.652306 & 0.020352549 & 314.58093
857 & 4 & 24 & 24 & 29.652306 & 0.020352549 & 388.23050











Out-of-Band Spectral Transmission Content[edit]

The HFI detector spectral transmission profiles have been extended beyond the optical pass-band of the detectors. This is done by using a combination of a waveguide model and external filter measurements for the out-of-band regions of the detector spectral response. An uncertainty estimate for these additional spectral regions is also provided, however, it should be noted that the spectral uncertainty for the waveguide and filter spectra are determined indirectly (as described above/below).

There is a transition from IAS FTS data to filter data for every band edge which is defined by an optical filter. For the 100, 143, 217, and 353 GHz bands this is the high frequency cut-off band edge. For the 545 and 857 GHz bands a separate filter is used to define each of the high and low frequency band edges. For the spectral regions outside of the HFI detector optical bands, first the IAS FTS data is used to qualitatively verify that there are no spectral leaks or features, and then the external filter measurements are grafted into the ratioed spectra where they better represent the relative spectral transmission. The spectral transmission profiles switch from the ratioed IAS FTS data to the filtered data as follows:

  1. The normalized ratioed spectrum is scaled by the optical efficiency (see §[sec:EFF]) whereas the filter spectrum remains normalized. This results in a conservative over-estimation of the out-of-band transmission.
  2. A lower threshold frequency, [math]\nu_{l}[/math], is defined for each band, below which the FTS data must be used.
  3. An upper threshold frequency, [math]\nu_{u}[/math], is defined for each band, above which the filter data must be used. The [math]\nu_{u}[/math] is defined to be twice the band cut-on frequency. This is done to avoid introducing any detector nonlinearity residuals into the final spectral transmission data products.
  4. The region between the lower and upper threshold frequencies is defined as the transition region.
  5. The transition point is defined as the first data point in the transition region where the ratioed FTS spectrum is greater in amplitude than the filter spectrum.
  6. If no such transition point exists, the filter spectrum is temporarily scaled to be smaller until a transition point is identified, after which the filter spectrum is restored to its original amplitude. If the temporary filter data scaling becomes less than 0.01 then the upper threshold frequency is defined as the transition point.
  7. The filter spectrum is used for the transmission profile for all frequencies greater than and equal to the transition frequency.
  8. The above method resulted in a mild distortion to the spectrum within the transition region for the 143 and 857 GHz spectra. This was a result of nonlinearity features in the FTS spectra just before the transition region cut-off. In these cases, a comparison between the FTS data derivative and filter data derivative was used to determine the transition point, instead of using the spectral amplitudes.

Additionally, a similar technique is used, with decreasing frequency instead of increasing, for the 545 and 857 GHz bands with a filter-induced frequency cut-on. Figure [fig:OOBex] illustrates an example of both the FTS and filter spectra used in extending the transmission profiles beyond the HFI optical bands. Similar plots for every detector are shown in §[sec:stitch].

Waveguide Model[edit]

A waveguide model is used to provide the data for the lowest frequency portion of the HFI detector spectral transmission. For the 100, 143, 217, and 353 GHz bands the waveguide model is transitioned (with increasing frequency) to the FTS ratioed spectra directly. There is an intermediate transition to the filter data, and then the ratioed spectra for the 545 and 857 GHz bands. For each detector, the waveguide transmission, [math]W(\sigma)[/math], is given by the following relation

[math]\label{eq:WGmodel} W(\sigma) = \left\{ \begin{array}{lcl} \exp{\left[ -\displaystyle\left(\sqrt{(c_{nm}/r_{w})^2 - (2\pi\sigma)^2}\right) \left( l_w \right) \right]} & , & \mbox{ for }\sigma \leq c_{nm}/(2\pi r_w) \\ 1 & , & \mbox{ for }\sigma \gt c_{nm}/(2\pi r_w) \end{array}\right.~,[/math]

where [math]\sigma[/math] is the frequency in cm[math]^{-1}[/math], [math]r_w[/math] is the waveguide radius in cm, [math]l_w[/math] is the waveguide length in cm, and [math]c_{nm}[/math] is a waveguide specific constant; 1.841 for the TE[math]_{11}[/math]/TM[math]_{11}[/math] (HE[math]_{11}[/math]) hybrid mode in this case . Table [tab:WG] lists the waveguide radii and lengths resultant from the waveguide model fit to the ratioed spectra. As all of the feedhorns for a given band are meant to be identical (i.e. within mechanical manufacturing tolerances) the uncertainty of the waveguide transmission is estimated statistically using all of the waveguide models for each band. I.e. for [math]n[/math] bands, the uncertainty at each spectral data point is determined by the standard deviation of [math]n[/math] transmission values at that frequency. As a result of each detector in a given band having a unique cut-on frequency, this method begins to over-estimate the uncertainty for frequencies approaching the cut-on; for regions very near the waveguide cut-on, the uncertainty is extrapolated from the ratioed spectrum as a more accurate representation.

|c|c|c|c|c|


Band (GHz) & bc & Det. & [math]r_w[/math] (mm) & [math]l_w[/math] (mm)

Band (GHz) & bc & Det. & [math]r_w[/math] (mm) & [math]l_w[/math] (mm)
100 & 00 & 1a & 1.039705 & 12.1000
100 & 01 & 1b & 1.038040 & 12.5813
100 & 20 & 2a & 1.038375 & 14.3000
100 & 21 & 2b & 1.041035 & 12.3750
100 & 40 & 3a & 1.042370 & 12.9250
100 & 41 & 3b & 1.042370 & 12.5125
100 & 80 & 4a & 1.033050 & 13.2000
100 & 81 & 4b & 1.033050 & 12.6500
143 & 02 & 1a & 0.740345 & 9.4500
143 & 03 & 1b & 0.737550 & 9.5625
143 & 10 & 5 & 0.736155 & 9.4500
143 & 30 & 2a & 0.740345 & 9.4500
143 & 31 & 2b & 0.739880 & 9.2250
143 & 42 & 6 & 0.741275 & 9.4500
143 & 50 & 3a & 0.741275 & 9.4500
143 & 51 & 3b & 0.741275 & 9.4500
143 & 60 & 7 & 0.729635 & 9.7875
143 & 70 & 8 & 0.739415 & 9.4500
143 & 82 & 4a & 0.739415 & 9.4500
143 & 83 & 4b & 0.738480 & 9.4500
217 & 04 & 1 & 0.4740750 & 7.8000
217 & 11 & 5a & 0.4814250 & 7.8000
217 & 12 & 5b & 0.4817310 & 8.0000
217 & 22 & 2 & 0.4749935 & 7.8000
217 & 43 & 6a & 0.4832625 & 8.0000
217 & 44 & 6b & 0.4832625 & 8.0000
217 & 52 & 3 & 0.4753000 & 8.0000
217 & 61 & 7a & 0.4820375 & 8.0000
217 & 62 & 7b & 0.4820375 & 8.0000
217 & 71 & 8a & 0.4838750 & 8.0000
217 & 72 & 8b & 0.4838750 & 8.1000
217 & 84 & 4 & 0.4740750 & 7.6000
353 & 05 & 1 & 0.2921850 & 6.47500
353 & 13 & 2 & 0.2914300 & 7.00000
353 & 23 & 3a & 0.2891650 & 7.70000
353 & 24 & 3b & 0.2893535 & 7.52500
353 & 32 & 4a & 0.2869000 & 7.35000
353 & 33 & 4b & 0.2870890 & 7.39375
353 & 45 & 7 & 0.2876550 & 6.12500
353 & 53 & 5a & 0.2936950 & 7.35000
353 & 54 & 5b & 0.2936950 & 7.61250
353 & 63 & 6a & 0.2929400 & 6.73750
353 & 64 & 6b & 0.2925625 & 6.62500
353 & 85 & 8 & 0.2889765 & 6.73750
545 & 14 & 1 & 0.3861000 & 3.600
545 & 34 & 2 & 0.3865875 & 3.525
545 & 55 & 3 & 0.3800065 & 3.525
545 & 73 & 4 & 0.3787875 & 3.525
857 & 25 & 1 & 0.2985 & 2.4
857 & 35 & 2 & 0.2985 & 2.4
857 & 65 & 3 & 0.2985 & 2.4
857 & 74 & 4 & 0.2985 & 2.4


Optical Efficiency[edit]

The EFF tests involved exposing the HFI detectors to a known blackbody source and observing the response. Sufficient details for the HFI detector spectral transmission profiles are provided here while full details of the EFF experiments and results are provided in a separate technical report . A blackbody source internal to the Saturne cryostat was set to a variety of temperatures ([math]\sim[/math]1 – 6 K) and the bolometer detector response was recorded. A bolometer model was applied to the recorded response in order to obtain the radiative optical power absorbed by the detector, in units of W, i.e. [math]P_{\mbox{abs}}(T_i)[/math] where [math]T_i[/math] represents the blackbody source temperature. Using the measured source temperature, the theoretical radiative optical power incident on the detector is also calculated using the Planck function. The ratio of the received power and the theoretical power provides the optical efficiency term. To remove any offsets in the measurement, a ratio of differences between unique temperature settings is used. The measured absorbed optical power difference is given by

[math]\label{eq:EFFdPabs} \Delta P_{\mbox{abs}} = P_{\mbox{abs}}(T_j) - P_{\mbox{abs}}(T_i) ~,[/math]

where [math]T_j[/math] and [math]T_i[/math] represent two unique source temperature settings. The theoretical incident power is determined using the HFI detector spectral transmission profiles. Let [math]\tau(\nu)[/math] represent the normalized detector transmission spectrum (i.e. it has been ratioed and had the waveguide model and filter data appropriately grafted). The spectral transmission is scaled for [math]\lambda^2[/math] throughput and then re-normalized as follows

[math]\label{eq:EFFrenorm} \tau'(\nu) = \mbox{Norm}\left[ \tau(\nu) \left( \frac{\nu}{c}\right)^2 \right]~,[/math]

where [math]\mbox{Norm}\left(f(x)\right)[/math] is the division of [math]f(x)[/math] by its maximum value, and [math]c[/math] is the speed of light. The normalized spectral transmission is then used with the Planck function at the temperature setting to determine the theoretical power, [math]P_{\mbox{th}}(T_i)[/math] , as follows

[math]\label{eq:EFFPth} P_{\mbox{th}}(T_i) = \displaystyle \int_{\nu_1}^{\nu_2}{ \left\{ \tau'(\nu) \left[ \displaystyle \frac{2 h \nu^3}{c^2 (\exp{(\frac{h\nu}{k T_i})} - 1)} \right] \left( \displaystyle \frac{c^2}{\nu^2}n_m \right) d\nu \right\}} ~,[/math]

where [math]h[/math] is the Planck constant, [math]k[/math] is the Boltzmann constant, the integration limits are given by [math]\nu_1[/math] and [math]\nu_2[/math], and [math]n_m[/math] is the expected mode content of the frequency band. Table [tab:modes] lists the [math]n_m[/math] values used for each band . In this case the integration is performed over the range [math]\nu \in [67~\mbox{GHz},1142~\mbox{GHz}][/math]. The difference between the theoretical power loading is given by

[math]\label{eq:EFFdPth} \Delta P_{\mbox{th}} = P_{\mbox{th}}(T_j) - P_{\mbox{th}}(T_i) ~,[/math]

which allows the optical efficiency term to be determined as follows

[math]\label{eq:EFFoptEff} \epsilon = \displaystyle \frac{\Delta P_{\mbox{abs}}}{\Delta P_{\mbox{th}}} ~.[/math]

Thus, if [math]\epsilon \tau'(\nu)[/math] were used in Equation [eq:EFFPth] in place of [math]\tau'(\nu)[/math], the resultant optical efficiency would be unity, indicating that the transmission losses have already been taken into account.

The uncertainty estimate of the optical efficiency is statistically based on the results from the multiple temperature settings used in the EFF test sequences.


<thead> </thead> <tbody> </tbody>
Band (GHz) 100 143 217 353 545 857
[math]n_m[/math] 1 1 1 1 3.4 8.3

CO interpolation[edit]

The IAS FTS/Saturne data taken with the HFI detectors is limited to a spectral resolution of [math]\sim0.017~\mbox{cm}^{-1}~(\sim0.5~\mbox{GHz})[/math] by the mechanical travel of the FTS translation stage. This corresponds to an unapodized FTS ILS FWHM of [math]\sim0.020~\mbox{cm}^{-1}~(\sim0.61~\mbox{GHz})[/math]. The previous iteration of the spectral transmission (IMO ver. 2_28), in addition to a triangular apodization kernel, has a factor of 3 sub-sampled spectral resolution to improve the signal-to-noise ratio (S/N) of the output spectra. In light of the CO contribution to the 100 GHz spectra, the spectral transmission has been reprocessed at the original full spectral resolution with a different apodization function which provides better noise reduction and less spectral resolution degradation, as discussed above.

In order to provide an improved estimate of the spectral transmission near the CO features, an interpolation of the spectra by a factor of [math]\sim[/math]10 has been performed. This over-sampling was accomplished by zero-padding the FTS interferograms prior to Fourier transformation , and subsequently incorporated into the Nyquist-sampled spectral data near the CO transitions (see Table [tab:CO]). Although the data are presented at higher resolution, the resolution of independent data points is not improved, i.e. the ILS line-width remains the same. A flag column has been added to the spectral transmission profile data files to indicate whether a given data point originates from the actual data, or is a result of the ILS-based interpolation. The region of the over-sampled, interpolated, data has been extended to also include other CO isotopes; the COJ1-0 – J9-8 transitions should be oversampled for CO, C[math]^{13}[/math]O, CO[math]^{17}[/math], and CO[math]^{18}[/math]. The original data points within the over-sampled region have been preserved (i.e. every tenth data point – the data point that is not an interpolated one – is flagged with a zero rather than a one), as is indicated by the data flag, so a flag filter on the data will restore the independent data points easily. Figure [fig:COflag] illustrates the regions where the over-sampling has been incorporated into the spectral transmission profiles. An example of the over-sampled spectra is shown in Figure [fig:COzoom100] for the 100 GHz detectors. Examples for the other detector bands are shown in §[sec:COappend]


|c|c|c|c|

<tbody> </tbody>
Band
(GHz)

& CO transition ([math]J_{\mbox{\tiny{upper}}} - J_{\mbox{\tiny{lower}}}[/math]) & [math]\nu_o[/math] (GHz) & over-sampled region (GHz)
100 & 1 – 0 & 115.2712018 & 109.67 – 115.39
217 & 2 – 1 & 230.5380000 & 219.34 – 230.77
353 & 3 – 2 & 345.7959899 & 329.00 – 346.15
545 & 4 – 3 & 461.0407682 & 438.64 – 461.51
545 & 5 – 4 & 576.2679305 & 548.28 – 576.85
857 & 6 – 5 & 691.4730763 & 657.89 – 692.17
857 & 7 – 6 & 806.6518060 & 767.48 – 807.46
857 & 8 – 7 & 921.7997000 & 877.04 – 922.73
857 & 9 – 8 &1036.9123930 & 986.57 – 1037.95


|c|c|c|c|c|



<tbody> </tbody>
Band
(GHz)

& BC & Det. & CO transition ([math]J_{\mbox{\tiny{upper}}} - J_{\mbox{\tiny{lower}}}[/math]) & Transmission


<tbody> </tbody>
Band
(GHz)

& BC & Det. & CO transition ([math]J_{\mbox{\tiny{upper}}} - J_{\mbox{\tiny{lower}}}[/math]) & Transmission
100 & 00 & 1a & 1 – 0 & 0.1495 [math]\pm[/math] 0.0054

100 & 01 & 1b & 1 – 0 & 0.1812 [math]\pm[/math] 0.0052
100 & 20 & 2a & 1 – 0 & 0.2577 [math]\pm[/math] 0.0119
100 & 21 & 2b & 1 – 0 & 0.1944 [math]\pm[/math] 0.0086
100 & 40 & 3a & 1 – 0 & 0.2798 [math]\pm[/math] 0.0133
100 & 41 & 3b & 1 – 0 & 0.1660 [math]\pm[/math] 0.0093
100 & 80 & 4a & 1 – 0 & 0.2538 [math]\pm[/math] 0.0114
100 & 81 & 4b & 1 – 0 & 0.2173 [math]\pm[/math] 0.0105

217 & 04 & 1 & 2 – 1 & 0.3781 [math]\pm[/math] 0.0026

217 & 11 & 5a & 2 – 1 & 0.4530 [math]\pm[/math] 0.0029
217 & 12 & 5b & 2 – 1 & 0.4174 [math]\pm[/math] 0.0030
217 & 22 & 2 & 2 – 1 & 0.3773 [math]\pm[/math] 0.0026
217 & 43 & 6a & 2 – 1 & 0.3149 [math]\pm[/math] 0.0023
217 & 44 & 6b & 2 – 1 & 0.3609 [math]\pm[/math] 0.0030
217 & 52 & 3 & 2 – 1 & 0.3924 [math]\pm[/math] 0.0025
217 & 61 & 7a & 2 – 1 & 0.3428 [math]\pm[/math] 0.0022
217 & 62 & 7b & 2 – 1 & 0.2770 [math]\pm[/math] 0.0017
217 & 71 & 8a & 2 – 1 & 0.4623 [math]\pm[/math] 0.0031
217 & 72 & 8b & 2 – 1 & 0.4340 [math]\pm[/math] 0.0031
217 & 84 & 4 & 2 – 1 & 0.3506 [math]\pm[/math] 0.0025

353 & 05 & 1 & 3 – 2 & 0.4487 [math]\pm[/math] 0.0036

353 & 13 & 2 & 3 – 2 & 0.5461 [math]\pm[/math] 0.0044
353 & 23 & 3a & 3 – 2 & 0.3443 [math]\pm[/math] 0.0030
353 & 24 & 3b & 3 – 2 & 0.4706 [math]\pm[/math] 0.0037
353 & 32 & 4a & 3 – 2 & 0.3099 [math]\pm[/math] 0.0024
353 & 33 & 4b & 3 – 2 & 0.2801 [math]\pm[/math] 0.0027
353 & 45 & 7 & 3 – 2 & 0.2923 [math]\pm[/math] 0.0022
353 & 53 & 5a & 3 – 2 & 0.3150 [math]\pm[/math] 0.0024
353 & 54 & 5b & 3 – 2 & 0.3181 [math]\pm[/math] 0.0023
353 & 63 & 6a & 3 – 2 & 0.2059 [math]\pm[/math] 0.0014
353 & 64 & 6b & 3 – 2 & 0.2113 [math]\pm[/math] 0.0017
353 & 85 & 8 & 3 – 2 & 0.3509 [math]\pm[/math] 0.0028
545 & 14 & 1 & 4 – 3 & 0.0747 [math]\pm[/math] 0.0003
545 & 34 & 2 & 4 – 3 & 0.0731 [math]\pm[/math] 0.0003
545 & 55 & 3 & 4 – 3 & 0.0521 [math]\pm[/math] 0.0002
545 & 73 & 4 & 4 – 3 & 0.0473 [math]\pm[/math] 0.0002

545 & 14 & 1 & 5 – 4 & 0.3306 [math]\pm[/math] 0.0012

545 & 34 & 2 & 5 – 4 & 0.3183 [math]\pm[/math] 0.0011
545 & 55 & 3 & 5 – 4 & 0.2428 [math]\pm[/math] 0.0009
545 & 73 & 4 & 5 – 4 & 0.2597 [math]\pm[/math] 0.0009

857 & 25 & 1 & 6 – 5 & 0.0280 [math]\pm[/math] 0.0001

857 & 35 & 2 & 6 – 5 & 0.0241 [math]\pm[/math] 0.0001
857 & 65 & 3 & 6 – 5 & 0.0292 [math]\pm[/math] 0.0001
857 & 74 & 4 & 6 – 5 & 0.0159 [math]\pm[/math] 0.0001

857 & 25 & 1 & 7 – 6 & 0.1636 [math]\pm[/math] 0.0005

857 & 35 & 2 & 7 – 6 & 0.1427 [math]\pm[/math] 0.0004
857 & 65 & 3 & 7 – 6 & 0.2176 [math]\pm[/math] 0.0007
857 & 74 & 4 & 7 – 6 & 0.1168 [math]\pm[/math] 0.0003

857 & 25 & 1 & 8 – 7 & 0.2554 [math]\pm[/math] 0.0009

857 & 35 & 2 & 8 – 7 & 0.2218 [math]\pm[/math] 0.0007
857 & 65 & 3 & 8 – 7 & 0.2744 [math]\pm[/math] 0.0009
857 & 74 & 4 & 8 – 7 & 0.1119 [math]\pm[/math] 0.0004

857 & 25 & 1 & 9 – 8 & 0.0053 [math]\pm[/math] 0.0000

857 & 35 & 2 & 9 – 8 & 0.0060 [math]\pm[/math] 0.0000
857 & 65 & 3 & 9 – 8 & 0.0085 [math]\pm[/math] 0.0000
857 & 74 & 4 & 9 – 8 & 0.0001 [math]\pm[/math] 0.0000

Conformity with requirements[edit]

The two defining requirements indicated in the HFI calibration plan are the acquisition of the spectral transmission of the single pixels with a prescribed accuracy and spectral resolution. The desired accuracy is 3% for the low frequency channels (CMB – 100, 143, and 217 GHz) and 1% for the high frequency channels (353, 545, 857 GHz). The spectral resolution requirement is for a resolution superior to 0.1 cm[math]^{-1}[/math].

The spectral resolution requirement has been exceeded by more than a factor of five. It is also possible to degrade the spectral resolution to the 0.1 cm[math]^{-1}[/math] requirement to gain an improvement in the S/N.

No quantitative number is present on the document regarding the blocking of high frequency (near IR, visible, UV) radiation outside the range of the instrument. Checks in order to quantify the rejection have been performed at subsystem level and estimates of the out-of-band transmission profiles have been incorporated into the data products.

The requirement on the spectral transmission accuracy for the CMB channels is of [math]3\%[/math]. There is also a requirement that the high frequency channels effective frequencies be known to within [math]1\%[/math]; knowledge of the transmission to 1% accuracy should allow this.

Considering statistical fluctuations in the determination of the spectra, these goals have been achieved. There are, however, caveats regarding the nature of error bars when dealing with frequency space. The nature of uncertainties in spectra determination is less obvious than when dealing with timestream data. Systematic effects produced from instrumental setup, but also by data reduction can in some cases exceed the actual statistical oscillation in the determination of the final spectra. This is the case for the high-frequency data for instance where the statistical fluctuation of the different determinations of the spectra in some case are better than 1 part in [math]10^3[/math].

A second caveat regards the method of data analysis of the calibration test data, for which the ratio with the reference bolometer data (of which the relative error is a function of frequency) introduces an error that increase with wavelength. With the spectral resolution of the data provided being much higher that the stated 0.1 cm[math]^{-1}[/math], the transmission accuracy requirement is not met for the 100 GHz detectors. It is possible to degrade the spectral resolution to the 0.1cm[math]^{-1}[/math] level to allow the accuracy to achieve the required level, but the higher spectral resolution data has been provided to better assist with the CO contamination removal from the 100 GHz signal.

(Planck) High Frequency Instrument

Focal Plane Unit

European Space Agency

Instrument Line Shape

Full-Width-at-Half-Maximum

[LFI meaning]: absolute calibration refers to the 0th order calibration for each channel, 1 single number, while the relative calibration refers to the component of the calibration that varies pointing period by pointing period.

LFI Waveguide

Cosmic Microwave background