7.5 MAMA Operation and Feasibility Considerations

7.5.1 MAMA Saturation—Overflowing the 16 Bit Buffer

The MAMA is a photon-counting detector: as each event is recorded, the buffer memory for the corresponding pixel is incremented by one integer. The buffer memory stores values as 16 bit integers; hence the maximum number it can accommodate is 65,536 counts per pixel in a given ACCUM mode observation. When accumulated counts per pixel exceed this number, the values will wrap, i.e., the memory resets to 0. As an example, if you are counting at 25 counts/s/pix, you will reach the MAMA "accumulation" limit in ~44 minutes.

Keep accumulated counts/pix below this value by breaking individual exposures into multiple identical exposures (see also Section 11.2.4), each of which is short enough that fewer than 65,536 counts are accumulated per pixel. There is no read noise for MAMA observations, so no penalty is paid in lost signal-to-noise ratio when exposures are split. There is only a small overhead for each MAMA exposure (Chapter 9).

7.5.2 MAMA Darks

The STIS MAMA detectors have intrinsically very low dark currents. Dark currents measured during ground testing were less than 10 counts/s for the FUV-MAMA and less than 30 counts/s for the NUV-MAMA over the whole detector. For the FUV-MAMA, this exceptionally low dark current was initially achieved on orbit. For the NUV-MAMA, charged particle impacts on the MgF₂ faceplate cause a higher background that resulted in a dark current of 800–2000 counts/s, varying both with temperature and the thermal history of the detector. This particular phenomenon is not present for the FUV-MAMA, but the dark current for that detector now also varies with time, temperature and position. The different dark current behaviors of the detectors are discussed in more detail below, and up-to-date information can be found on the Monitoring page of the STIS website. We split the discussion of the NUV dark current into before and after servicing mission 4 (SM4) where the electronics were switched and STIS become operational again.

NUV-MAMA Dark Current Before SM4

Most of the dark current in the NUV-MAMA comes from phosphorescence of impurities in the MgF₂ detector faceplate. A simple model of the phenomenon was developed by Jenkins and Kimble that envisions a population of impurity sites each having three levels: (1) a ground state, (2) an excited energy level which can decay immediately to the ground state, and (3) a meta-stable level that is at an energy slightly below the one that can emit radiation. The meta-stable states are initially populated by charged particle impacts that mostly occur during passages through the South Atlantic Anomaly (SAA). Hours or days later, the electrons trapped in these meta-stable states are thermally excited to an unstable upper level and then emit a photon as they decay to the ground state. This thermal excitation rate is proportional to exp(−ΔE/kT), where ΔE is the energy difference between the levels, T is temperature of the detector and k is the Boltzmann's constant. The behavior of the count rate vs. temperature leads to an estimate of 1.1 eV for ΔE.

While this temperature dependent function does a good job of predicting the short term response of the NUV dark current to temperature changes, the longer term response is more complex. In equilibrium, the number of decays will match the number of excitation events. A sudden temperature increase would then result in an initial rapid increase in the dark rate, as the meta-stable states are more easily depopulated at higher temperatures. However, after several days, the population of meta-stable states would reach a new equilibrium, resulting in a dark rate that, while higher than the equilibrium rate was at the cooler temperature, is significantly lower than the short term response to the same temperature increase. The predicted large temporary increases in the dark rate were observed after the initial STIS installation, and also after Servicing Mission 3a (see Figure 7.18).

At a fixed detector temperature of 30° C, the time-constant for the dark current to reach an equilibrium value is about 8 days. In practice because the MAMA high-voltage power supplies have to be shut down during SAA impacted orbits, the detector temperature varies from about 27° C to 40° C on a roughly daily time scale, and the dark current never reaches equilibrium. In addition, the real rate of excitations caused by charged particle impacts is not directly measurable and is expected to vary unpredictably on all time scales. Even the time averaged excitation rate may differ considerably from year to year, depending on the low earth orbit radiation environment.

During the last nine months of STIS operations prior to the failure in August 2004, the typical NUV-MAMA dark current ranged between 8 × 10⁻⁴ and 1.6 × 10⁻³ counts/pix/s (Figure 7.19). Changes in the low earth orbit radiation environment affect the rate at which the meta-stable impurity states are populated. For Cycle 17 planning purposes we had guessed that the equilibrium STIS NUV-MAMA dark current would be 10% higher than it was during 2003–2004; which implied a range of dark count rate values between 8.5 ×10⁻⁴ and 1.7 x 10⁻³ counts/pix/s. However, as we will see in the next section (NUV-MAMA Dark Current after SM4), this drastically underestimated the actual dark rate seen after SM4.

To predict the NUV dark current and subtract it as part of standard processing, the calstis pipeline uses a simple temperature dependent relation:

darkrate = norm · 1.805 × 10²⁰ · exp(−12211.8/max(T,Tmin))

Both norm and Tmin are slowly varying functions of time that are empirically adjusted to give a good match to the observed dark rate, and which are tabulated in the temperature dependent dark correction table (TDCTAB) reference file. The temperature for a given observation is taken from the OM2CAT telemetry value, which is included in the extension header of each MAMA observation. This approach usually predicts the dark rate with 5 to 10% accuracy, although the error in individual cases may be substantially larger.

Because 99% of the NUV-MAMA dark current is due to photons hitting the detector, it is appropriate to apply the flat field prior to subtracting the dark current. The dark current varies slowly across the face of the detector, being about 1.25 times higher near the lower left corner (AXIS1, AXIS2 < 300) than at the center. This shape varies with time and temperature enough that subtraction of the scaled dark reference may leave a residual, spatially varying dark current. This is easiest to remove by fitting a low-order two-dimensional function to the background vs. pixel position.

NUV-MAMA Dark Current after SM4

After the Side-2 failure in 2004, STIS was in safe mode with only the survival heaters powered. During this period, the MAMA tubes were at an average temperature of 3.8° C, much lower than the normal operational temperature range of about 29° C to 34° C. We therefore expected to see a temporary increase in the dark current after the post-SM4 recovery of the MAMA tubes due to the increased population of meta-stable states at the colder safe-mode temperature. Detailed modeling led us to expect that the peak dark count rate would not be too much bigger than that seen after SM3a where the MAMAs had been off for about three months. We also expected that any excess would decline with a e-folding time of a week or so, as had been seen previously. This would have allowed the NUV-MAMA dark rate to return to the previous operational range of between 0.0009 and 0.0017 counts/pix/s within about a month of the STIS recovery after SM4.

Instead we found initial dark count rates for the NUV-MAMA as high as 0.016 counts/pix/s (Figure 7.20). As was the case before SM4, the dark current shows short time scale variations that depend exponentially on the detector temperature, with the overall level also varying on longer time scales. However, for the new component of the dark current, the decay time-scale appears to be much longer than the one to two weeks previously seen. Initially the excess dark current appeared to be declining with an e-folding time of about 100 days; however over subsequent months, the rate of decline has continued to slow. By mid-2010 the dark rate was fluctuating between about  0.0025 and 0.0042 counts/pix/s, and it appears that additional significant declines will take years rather than months. For planning purposes and the STIS ETCs, we will adopt a mean  NUV-MAMA dark current of 0.0013 counts/pix/s for mid-Cycle 32 (April 2025). Updates on the current state of the dark current are posted on the STIS MAMA dark monitor pages.

FUV-MAMA Dark Current

The FUV-MAMA dark current is substantially lower than that of the NUV-MAMA. Initially, values as low as 7 counts/s across the face of the detector (7 × 10^–6 counts/pix/s) could be routinely expected. However, there is also an intermittent glow that covers a large fraction of the detector (see Figure 7.21). The source of the dark current is not phosphorescence, but is intrinsic to the micro-channel plate array (it was seen in ground testing). This glow can substantially increase the dark current over a large fraction of the detector, and this leads to count rates of up to 300 counts/s integrated across the face of the detector.

An example of the dark current variation across the detector can be seen in Figure 7.21, which is the sum of a number of 1380-second dark frames taken during periods of high dark current. The dark current in the lower right quadrant appears to be stable to within 10% over time. The dark current in the upper left quadrant varies with time and temperature. The total dark current can be approximated by the sum of a constant dark current plus a "glow" image, scaled to the net rate in the upper left quadrant.

During the first two years of STIS operations, this glow was only present intermittently, but since mid-1999 it has been present more often than not. The glow increases with the amount of time the detector high-voltage has been on since the last SAA passage. The rate of increase is greater when the detector is warmer, and it has also been increasing from year to year, even under otherwise comparable conditions. The dark monitor measurements for the first few months of STIS operations after the repair during SM4 are shown in Figure 7.22, while long-term monitor measurements are shown in Figure 7.23.

In Figure 7.22, the FUV-MAMA dark current is shown for 35 dark monitor observations taken between 2009-Jun-09 and 2009-Oct-01. On days that the FUV-MAMA is used, the MAMA high voltage is turned on for the block of HST orbits that is unaffected by SAA passages, and as the detector warms up, the glow increases. Filled circles show the count rate in a region centered on the peak of the dark glow (pixels [200:400,600:800]), while the open circles show the average rate for the entire detector. Note that the lower right corner of the detector (pixels [900:1000,10:110]), denoted by plus signs, shows little or no increase over the initial low dark current values.

Figure 7.23 shows the long-term FUV dark current before and after SM4. The top plot shows the average dark rate over time for the full detector, while the bottom plot is restricted to the glow region. The different colors show how the dark rate trends vary with high voltage on-time. The dark current varies depending upon the detector region, detector temperature, high voltage on-time, as well as a long-term increase. The vertical distribution of points on the graph changes in 2013, but this is merely due to a change in scheduling of calibration observations relative to the high-voltage turn-on. After SM4, the overall level of the glow appears to be similar to what had been seen in 2004, but because the physical origin of the FUV detector glow is poorly understood, it is difficult to predict what the future behavior will be. Past experience would indicate that the overall level of the glow will continue to increase.

For Cycle 32 planning purposes, we recommend assuming that in the center of the glow region the mean dark current will start with a value near 3.0 × 10⁻⁵ counts/pix/s at the beginning of each SAA-free block of orbits, but will then increase at a rate of about 6 × 10⁻⁵ counts/pix/s/HST-orbit, reaching peak values as high as 5 × 10⁻⁴ counts/pix/s. In the dark corner of the detector, the count rate will remain very low, with a rate near 3.0 × 10⁻⁵ counts/pix/s. The STIS ETC currently defaults to 7.5 × 10^-5 counts/pix/s, as representative of the expected mean rate averaged over the whole detector in the first 2-3 orbits. While the ETC can only reflect the average of the dark current rate, the STIS FUV dark current rate in particular exhibits tremendous variations with the position on the detector. STIS FUV-MAMA users whose observations are sensitive to dark current (e.g., faint targets) are strongly encouraged to read the corresponding documentation to assess the feasibility of their observations and better constrain the exposure time needed to achieve the required accuracy.

The count rates discussed above refer to the native MAMA detector pixels. When comparing to ETC results these values need to be scaled by the number of pixels included in an extracted spectral "pixel" or resolution element. For observers of bright targets, even the highest FUV-MAMA dark rates will have a negligible impact on data quality. For observers of very faint targets, there are a number of possible mitigation strategies. For spectra of faint point sources, the observer should consider whether one of the COS FUV channels might be a better choice, as these have both much higher overall throughput and much lower detector background than comparable STIS modes. If STIS spectral observations are needed (e.g., to provide higher spatial resolution in the cross dispersion direction), the dark rate can also be minimized by putting the source near the bottom of the FUV-MAMA detector below the worst of the extra detector glow. The "D1" aperture positions are intended for this purpose (see Section 4.2.3). If neither of these strategies is adequate, the  proposers should contact the HST Help Desk to discuss whether it may be practical to schedule their observations immediately after the detector high voltage has been ramped up, when the intrinsic detector background is lowest. Because of the need to schedule STIS MAMA operations around the South Atlantic Anomaly (SAA), in practice, the FUV-MAMA detector is usually turned on only one time per day, and there is usually only about one orbit per day available where a detector background close to 1 × 10^–5 c/s/pixel can be expected, so it may not be possible to satisfy all such requests. Justification for this scheduling constraint must be provided under the Special Requirements section of the the Phase I proposal.

Time-varying dark current for the FUV-MAMA is not currently subtracted by the pipeline. Some "average glow" images are available for off-line reduction to reduce the effects of the glow region. Users should contact the HST Help Desk to obtain access to these images.

Because the dark current is so low in the MAMA detectors, a typical STIS FUV-MAMA observation will have less than one count/pix from the dark. It is good to keep this in mind when reducing the data, as various standard measures of background (the median for example) are not good estimates when the data are quantized into just a few values. The best way to estimate the background is to identify hot pixels using the standard reference files, then use an unclipped mean for the remaining pixels in a source-free region of the image.

7.5.3 MAMA Signal-to-Noise Ratio Limitations

MAMA detectors are capable of delivering signal-to-noise ratios (S/N) of the order of 100:1 per spectral resolution element or even higher. Tests in orbit have demonstrated that such high S/N is possible with STIS (Kaiser et al., 1998, Proc. SPIE, 3356, 415; Gilliland, STIS ISR 1998-16.)

High S/N observations of several standard stars were obtained during STIS commissioning, and they were reduced with flats obtained during preflight testing of the detectors. Signal-to-noise ratios of 125 and 150 per spectral resolution element (for an 11 pixel extraction height in the cross dispersion direction) were achieved for the FUV- and NUV-MAMA observations, respectively; see Chapter 12 for a more detailed discussion.

For targets observed at a fixed position on the detector, the signal-to-noise ratio is limited by systematic uncertainties in the small-scale spatial and spectral response of the detector. The MAMA flats show a fixed pattern that is a combination of several effects including beating between the MCP array and the anode pixel array, variations in the charge-cloud structure at the anode, and low-level capacitive cross-coupling between the fine anode elements. Intrinsic pixel-to-pixel variations are 3.9% and 2.8% rms for the FUV- and NUV-MAMA, respectively, in 1024 × 1024 pixel format. In the highres 2048 × 2048 format (Section 3.2.1) the intrinsic variations are much larger. This fixed pattern appears to be stable at the 1–2% level.

The structure of the flat may vary slightly for different modes due to different incidence angles of the incoming photons on the microchannel-plate pores. A possible example of this is shown in Figure 7.24, where for an echelle E140H observation, the measured fluxes in the overlapping spectral range of two adjacent orders are displayed. The individual orders show broad features with widths of 10 to 30 pixels and amplitudes of 2 to 4% that are not reproduced in the other order at the same wavelengths. No corresponding features appear in the flat fields at these locations. These features also appear to be very sensitive to the exact placement of the spectrum on the detector, as other high signal-to-noise observations at the same CENWAVE setting do not show them. At wavelengths that are not covered by multiple spectral orders, the only way to distinguish this kind of detector artifact from a weak feature in the target spectrum is to use some kind of dithering to place the feature at multiple locations on the detector (see Section 12.5.1).

Observing strategies for achieving spectral S/N higher than ~100:1 are discussed in Chapter 12. For echelle mode spectra, observers may want to consider the use of the FP-SPLIT slits, or use multiple CENWAVE settings with overlapping wavelength ranges; for first-order mode observations, they may wish to dither the target along the slit.

Since MAMA observations can be binned in post-observation data processing with no additional signal-to-noise price, the option to obtain MAMA observations with unequal binning (e.g., BINAXIS1=1, BINAXIS2=2) was disabled starting in Cycle 8.

7.5.4 MAMA Nonlinearity

Global

The MAMA detectors exhibit nonlinearity (photon incidence rate not equal to photon count rate) greater than 5% at global (across the entire detector) count rates of order 200,000 counts/s. The nonlinearity can be approximately corrected in post-observation data processing for rates up to ~285,000 counts/s.  The MAMA detectors (plus processing software) are not able to count reliably at rates exceeding 285,000 count/s.  For this reason and to protect the detectors from damage, observations at such high rates are not allowed (see Section 7.7).

Local

The MAMA detectors remain linear to better than 1% in their counting up to ~220 counts/s/pix for the FUV-MAMA and 340 counts/s/pix for the NUV-MAMA. At higher rates, they experience local (at a given pixel) nonlinearity. The nonlinearity effect is image dependent—that is, the nonlinearity observed at a given pixel depends on the photon rate affecting neighboring pixels. This property makes it impossible to correct reliably for the local nonlinearity in post-observation data processing. In addition, the MAMA detectors are subject to damage at high local count rates (see Section 7.7).