Demander un devis

CCCP: A CCD controller for counting photons

Olivier Daiglea,b,c, Jean-Luc Gachb , Christian Guillaumed , Simon Lessardc , Claude Carignana,b,e , Sébastien Blais-Ouellettec
a Laboratoire d’Astrophysique Expérimentale, Département de physique, Université de Montréal, C.P. 6128 Succ. Centre-Ville, Montréal, QC, Canada, H3C 3J7;
b Laboratoire d’Astrophysique de Marseille, Observatoire Astronomique de Marseille-Provence, Technopôle de Château-Gombert, 38, rue Frédéric Joliot-Curie, 13388 Marseille, France;
c Photon etc., 5155 Decelles Avenue, Pavillon J.A Bombardier, Montréal, Québec, Canada, H3T 2B1;
d Observatoire de Haute-Provence, 04870 St-Michel l’observatoire, France;
e Observatoire d’Astrophysique de l’Université de Ouagadougou, BP 7021, Ouagadougou 03, Burkina Faso.

Abstract

CCCP, a CCD Controller for Counting Photons, is presented. This new controller uses a totally new clocking architecture and allows to drive the CCD in a special way. Its design is optimized for the driving of EMCCDs at up to 20MHz of pixel rate and fast vertical transfer. Using this controller, the dominant source of noise of EMCCDs at low flux level and high frame rate, the Clock Induced Charges, were reduced to 0.001 – 0.0018 electron/pixel/frame (depending of the electron multiplying gain), making efficient photon counting possible. CCCP will be deployed in 2009 on the ESO NTT through the 3D-NTT1 project and on the SOAR through the BTFI project.

Plus d'informations sur les applications de Nüvü Camēras en espace et défense ici

1. INTRODUCTION

Electron Multiplying Charge Coupled Devices (EMCCD) allows one to apply a gain to the pixel’s charge before it reaches the noisy output amplifier where the charge-tension conversion is made.2 A gain G in the charge domain affects the effective readout noise by the relation σeff = σreal / G . Sub-electron effective readout noise levels are thus achievable. However, the electron multiplying process is stochastic. This statistical behaviour adds an excess noise factor that reaches a value of 21/2 at high gains.3 The effect on the signal-to-noise ratio (SNR) of the system is the same as if the quantum efficiency (QE) of the CCD would be halved.

Some authors proposed offline data processing to lower the excess noise factor induced by the multiplication register.4, 5 However, the only way to overcome the excess noise factor without any a priori knowledge or stability assumption on the signal is to consider the pixel binary by applying a single threshold to the output signal. This way, only one photon per pixel per frame can be counted and the full QE of the silicon can be recovered, making the EMCCD a theoretically perfect photon counting device. The highest flux rate that can be observed in this mode will thus depend of the frame rate at which the EMCCD is operated. However, charges are generated as the EMCCD is read out. Clock Induced Charges (CIC), a well know source of noise affecting all kinds of CCDs, were typically measured in the range of 0.1 to 0.01 electron per pixel per frame6–8 (for a 512 x 512 CCD97 frame transfer EMCCD from E2V Technologies) and quickly dominate the dark current or even the photon flux as the frame rate in increased. Thus, in order to make photon counting efficient at low flux with an EMCCD, the CIC must be reduced to a minimum. Some techniques were proposed to reduce the CIC8–13 but so far, no commercially available CCD controller was able to implement all of them and get efficient results.

CCCP, a CCD Controller for Counting Photons, has been designed with the aim of reducing the CIC generated when an EMCCD is read out. It is optimized for the driving of EMCCDs at high speed, both vertically and horizontally, but may be used for driving classical CCDs as well. Using this controller, CIC levels as low as 0.001 – 0.0018 event per pixel per frame (as opposed to per transfer) were measured on the 512 x 512 CCD97 EMCCD from E2V Technologies operating in inverted mode (as opposed to non-inverted mode). The impact of this level of CIC on the photon counting efficiency of an EMCCD will be discussed in this article. Data gathered using the controller will also be presented.

2. PHOTON COUNTING WITH AN EMCCD

Throughout this article, photon counting is referred to as being the process by which the output signal of the EMCCD is thresholded to a single value. Pixels having an output value higher than the threshold are considered having undergone one and only one event. Pixels having an output value lower than the threshold are considered not having undergone an event. This processing is opposed to the analogic processing, where the output signal of the EMCCD is divided by the mean gain of its EM register to allow more than one event per pixel per frame to be considered. Operating an EMCCD in photon counting mode allows the excess noise factor to be reduced to a value of 1.

2.1 Effect of gain and readout noise

Sub-electron readout noise is not necessarily synonym of efficient photon counting. When one takes a look at the histogram of an EMCCD operated under low flux (left panel of figure 1), he realizes that a significant amount of events may be hidden in the readout noise, below the threshold. The proportion t of events lost due to a cut level, cut (expressed in electrons), may be calculated by means of the following convolution:

where f(n, λ) is the Poissonian probability of having n photons during an integration period under a mean flux of λ (in photon/pixel/frame) and p(x, n, G) is the probability of having x output electrons when n input electrons are present at the input of the EM stage at a gain of G. This probability is defined by

Intuitively, the higher the cut level, the more events lost. Since the cut level is determined by the real readout noise, the factor to optimize will thus be the ratio of the gain over the readout noise. Right panel of figure 1 shows the relation between the ratio and the proportion of counted events. In order to count ∼ 90% of the events, a ratio of 50 must be achieved. For a readout noise of 60 electrons (typical value at 10MHz of pixel rate), a gain of 3000 is necessary. However, one can not increase the gain without limit. The CIC will also increase with the gain, as discussed in section 3.2.

The choice of the threshold is important and is resumed by figure 2. Choosing a threshold that is too low yields the counting of pixels whose values above the threshold are due solely to the readout noise, as shown by the left panel of figure 2. However, choosing a threshold that is too high yields the missing of real events (right panel). Given that the CIC level achieved with CCCP is about 0.001 event per pixel per frame, which represents the lowest count rate that will be seen in an image, a threshold of 5 σ will cause a maximum of one pixel out of a million to have a noise event, while allowing nearly 90% of the events to be counted (left panel of figure 1), depending of the gain/noise ratio. Figure 3 summarizes the fate of photons in photon counting mode.

Figure 1. Left: Histogram of an EMCCD operated under low flux, at an EM gain of ∼2750. Only a few pixels underwent more than one event per frame. The vertical dotted line shows the threshold at 5.5σ. The mean event rate is 0.0018 event per pixel per image. Right: Proportion of counted events as a function of the ratio of the EM gain over the readout noise. A cut level of 5σ is used. Values for ratios of 10, 20, 30, 50 and 100 are printed.
Figure 2. Left: Proportion of noise events that are counted as real events in photon counting mode, as a function of the threshold level expressed in multiples of σ, the readout noise. Right: Proportion of real events that are counted in photon counting mode, for a gain/readout noise ratio of 50, as a function of the threshold level expressed in multiples of σ.
Figure 3. Proportion of counted and lost events events in photon counting mode, for a gain/noise ratio of 50 and a 5σ threshold.

2.2 Coincidence losses

The drawback of the photon counting operation is that two events occurring during a single integration time will be counted as only one. Thus, events will be lost. Poissonian statistics allows one to count g, the proportion of counted photons as a function of λ, the mean number of photons expected during the integration period:

From that equation, in order to be able to count more than 90% of the events, the expected flux should not be higher than 0.2 event per pixel per frame.

2.3 CIC: the dominant noise source

Efficient photon counting with an EMCCD requires low CIC. CIC generated during the vertical transfer is dependant, among other things, of the operation mode of the CCD, namely inverted or non-inverted. As specified in Ref. 8, the amount of CIC generated during the vertical transfer could be lowered by a factor of ∼30 by switching from inverted to non-inverted mode, at the price of an increased dark noise. However, the surface dark current, which is suppressed by the inverted operation, is expected to dominate the bulk dark current by a factor of ∼200 at cryogenic temperatures. Thus, a reduction of a factor of 30 in CIC comes at the price of an increase of hundreds in dark signal.

Figure 4 shows how badly the CIC affects the SNR of an observation in photon counting mode. The left panel shows the simulation of a device in Inverted Mode Operation (IMO), while the right panel show the simulation of a device in Non Inverted Mode Operation (NIMO). The dark current used for the IMO plot is the one measured on a CCD97 at -85C (see temperatures considerations in section 3.4) and the dark current used for the NIMO plots is calculated from equation 1 in Ref. 8 for a 16×16µm pixel at -85C (which is the size of the pixel of the CCD97). From these figures, IMO is clearly the operation mode to favour if one is able to achieve CIC levels in the range of 0.001 event/pixel/image. Very little gain in SNR could be achieved from driving the CCD in NIMO even if that would reduce the CIC further, since the noise would be dominated by the dark noise. It is not expected that the driving of the CCD in NIMO with CCCP would reduce the CIC since the CIC generated during the vertical transfer is very low (see section 3.1 for more details).

Figure 4. Effect of the CIC on the SNR of an observation, compared to a perfect photon counting device having the same QE, whose noise is solely the shot noise. The simulations assume a device running at 30 frames per second, CIC is expressed as event/pixel/frame and coincidence losses are taken into account. Left: IMO operation: dark noise of 0.001 electron per pixel per second. Right: NIMO operation: dark noise of 0.02 electron per pixel per second. See text for details on the values used for the dark noise.

3. RESULTS

This section presents data obtained with CCCP using a scientific grade CCD97 EMCCD from E2V Technologies operated at a pixel rate of 10MHz in IMO. Event rates presented in the figures represents not only the counted events (events above the photon counting threshold), but all the events generated during the read out process, including those buried in the readout noise. In order to achieve this, event rates are calculated by fitting the output histogram with the output probability function (equation 2) of the EM stage. Least square minimization is then made to find the exact parameters of the output signal. This yields the EM gain and the mean event level at the same time and most importantly, this allows the exact event rate to be determined for low gain/readout noise ratios (recall figure 1, right panel).

3.1 CIC in the vertical transfer

Even when operated in Inverted mode, CCCP manages to keep the CIC low during the parallel transfer, as shown by figure 5. The flatness of the plots is something that is expected: the high voltage phase plays no role in the vertical transfer. The small effect of the temperature (left panel) is due to the dark noise that is generated during the read out of the CCD, as this process takes ∼30 milliseconds. Thus, at higher temperatures, more events are seen, even for a 0 second integration. The dark component has been suppressed from the figure in the right panel and curves in this figure should represent only the CIC generated during the vertical transfer.

Figure 5. Measurement of the events generated during the parallel transfer, for a 512 x 512 frame transfer device versus temperature and gain. The label shows the temperature, in Celcius, at which the data was acquired. Obviously, one does not expect the gain to affect the amount of vertical events. Left: This figure has not been corrected for the dark events that are generated during the readout process. Right: This figure is corrected for the dark events and should represent only the CIC events generated by the vertical transfers.

 

3.2 CIC in the horizontal transfer

The high voltage phase of the EMCCD is meant to produce impact ionization and multiply the pixel’s charge. Electrons may however be generated even in the absence of an electron at the input of the EM register. This CIC generated in the horizontal register will produce charges appearing at the output of the EM register that may be above the threshold. Dark events generated in the horizontal register will also undergo the EM amplification and appear as photon events at the output. Thus, the amount of events generated in the serial register will depend mostly on EM gain. Figure 6 shows these relations. Temperature does not play a significant role in the amount of events generated during the horizontal transfer.

Figure 6. CIC generated in the horizontal register as a function of both the EM gain and the temperature.

3.3 Total CIC

The total CIC is considered to be the sum of vertical and horizontal CIC. In fact, the vertical CIC presented in figure 5 is computed from the total CIC minus the horizontal CIC. Thus, the data presented at figure 7 correspond to what is actually seen when the image section of the device is read. Even at EM gains as high as 4000, the total CIC measured is less than 0.002 event/pixel/frame. By comparing this figure with figures 5 and 6, one sees immediately that the horizontal CIC is dominating over the vertical one at all gains (> 800). Given its strong gain dependance, the horizontal CIC is mostly due to the high voltage clock.

Figure 7. Measurement of all the events generated during the read out process. Left: including dark noise. Right: excluding dark noise.

3.4 Temperature and CTE

It is very tempting to lower the operating temperature of a CCD below -100C in order to reduce the dark noise to a minimum. However, the Charte Transfer Efficiency (CTE) in the horizontal register degrades very quickly as the temperature is lowered. The left panel of figure 8 shows the amount of bad events that are seen in an image as a function of both the gain and the temperature. Bad events are defined as being an event (> 5σ) immediately followed (pixel-wise) by another. This kind of event should not occur more than once every event rate, which is only the CIC rate in this case (0.001 - 0.002 event/pixel/image). However, at temperatures below -100C, this can account for more than 10% of the events at high gain. This is due to the events that are ”leaking” into neighboring pixels as they are shifted. This phenomenon is shown by the top right panel of figure 8, where the energy distribution of the events is plotted as a function of temperatures. Each line represents the mean energy contained in a pixel at a given offset from the event (the event is at offset 0). At -105C, less than 80% of the energy is contained in the pixel at offset 0. at -85C, this reaches 95%.

This phenomenon would be hard to correct, since the dispersion of the energy distribution is considerable, as show by the bottom right panel of figure 8. The standard deviation of the value of the neighboring pixels around an event is at least equal to the pixel value. Thus, it would not be possible to use a post-processing based on the mean profile to compensate for the bad CTE. The shape of the profile of an event will hardly be the same as the shape of the mean profile. Thus, at high gain, temperatures below -90C should be avoided. Simulations are needed to find the perfect operating temperature, but -85C seems to be a good compromise.

Figure 8. Effect of temperature on the EMCCD. Left: Fraction of bad events (see text for a definition of bad events) seen in EMCCD images as a function of the EM gain and the operating temperature. Right top: Event mean energy distribution as a function of temperature. Different lines corresponds to different offset from the main (strongest) pixel of the event. Measured gain values were 3195, 2690, 2740, 2900, 2800, 2720 and 2715 for temperatures -105C, -100C, -95C, -90C, -85C, -80C and -75C, respectively. Right bottom: Dispersion of the energy distribution. If all events are normalized to the value of their strongest pixel, the dispersion of the proportion of the energy contained in the neighboring pixels are as plotted.

3.5 Gain stability

The stability of the EM gain is not very critical for photon counting operation. However, for amplified operation, it is mandatory to have a good gain stability to ensure photometric continuity across multiple images. Two factors affects the EM gain: the high voltage phase and the temperature. From data shown in figure 9, left panel, one sees that in order to have a ±1% stability on the gain, one must have no more than ±0.14C temperature variation and ±5mV variation on the high voltage clock (at -85C, 43.48V HV clock, which gives a gain of ∼3900). However, at lower gains, the constraint on the relative gain variation per C and per mV relaxes.

The right panel of figure 9 shows data acquired over time, showing a gain variation of ∼ ±1%. The two dotted lines shows the variation that is expected from the temperature stability of the test dewar, which is about 0.3C peak-to-peak. From this data alone, one can not tell if the stability of the high voltage clock is sufficient; gain variation due to the high voltage clock variation could be hidden in the temperature variations. However, this data shows that precise temperature control is mandatory if less than 1% variation on the gain is expected, at high gain. Measurements of the high voltage clock with an oscilloscope showed no more that 5mV variation through an image and over hours of operation. In fact, the oscilloscope was the limiting factor in this measurement.

Figure 9. Left: Effect of the temperature and the high voltage clock on the gain. At -85C and gain of ∼3900, the gain variation is 7% per C and ∼0.2% per mV. Right: Stability of the EM gain through time. The temperature variation was measured to be 0.3C peak-to-peak around -85C and the expected gain variation due to the temperature is shown by the two dotted lines.

4. CCCP EFFICIENCY

It is of interest to compare the expected efficiency in photon counting mode of an EMCCD driven by CCCP with the photon counting systems actually in operation. In order to have zero readout noise, GaAs photocathode-based image amplifiers are often placed in front of a fast read-out CCD.14, 15 Photons hitting the photocathode are amplified several hundreds of thousands times and produce a bright spot of a few pixels wide on the imaging CCD. Centering must then be made on these spots to recover the exact location of the incident photon. Coincidence losses on these Image Photon Counting Systems (IPCS) is thus higher than on an EMCCD in photon counting mode: two incident events that are located near from one another will be counted as only one. However, the IPCS do not suffer from CIC and their dark noise is typically an order of magnitude lower than the one of a CCD. The biggest drawback of these systems is their limited Quantum Efficiency (QE): the photocathode itself has a peak QE of ∼25%. QE of EMCCDs are the same as classical CCDs and can be higher than 80% on a wide spectral range (450-800nm) and peak at >95%.

Figure 10 compares the efficiency of such a photocathode-based system and the CCD97 driven by CCCP. The characteristics of each system is given in the caption of the figure. When taking into account all the sources of noise (dark noise, CIC, readout noise) and losses (coincidence, threshold), the left panel of the figure shows that en EMCCD will outperform an IPCS for incident fluxes higher than 0.022 photon per pixel per second (1 photon per 45 seconds). Thus, for the same pixel size, the gain in observing efficiency will be as shown in the right panel. The observing efficiency is defined as being the time it takes to reach a given SNR at a given flux. The efficiency of the EMCCD in IMO at low flux could be raised by lowering the frame rate since the noise is dominated by the CIC in this flux regime. NIMO operation of the EMCCD with CCCP would not benefit from the lowering of the frame rate as the image would be quickly dominated by dark noise (figure 4).

Figure 10. Comparison between CCCP in photon counting mode and photocathode-based Image Photon Counting Systems (IPCS). Parameters for CCCP: gain of 3000, readout noise of 60 electrons, CIC of 0.0015 electron/pixel/frame, dark noise of 0.001 electron/pixel/second (IMO), threshold of 5σ, 30 frames per second, quantum efficiency of 80%. Parameters for the IPCS: readout noise of 0 electron, dark noise of 0.0001 electron/pixel/second, 60 frames per second, quantum efficiency of 20%, coincidence losses are expected to occur if two events fall in the same 3x3 pixels box. The dotted line marks the flux of equal SNR (0.022 photon/pixel/second). Left: SNR of both systems compared to a perfect photon counting system (QE 100%, only shot noise). Right: Compared observing time efficiency of CCCP vs the GaAs IPCS. The dotted line show a relative efficiency of 1.

5. CONCLUSIONS

Tests made with CCCP on a CCD97 EMCCD shows that the CIC can be greatly reduced without having to resort to Non Inverted Mode Operation. The low level of CIC of 0.001 – 0.0018 event/pixel/frame, depending on the gain, allows one to use an EMCCD in photon counting mode and be more efficient than a GaAs IPCS at fluxes higher than ∼0.02 photon/pixel/second.

Analogic (amplified, no threshold) and conventional (no amplification) operation of the EMCCD with CCCP is also possible: the controller has the possibility to read the two outputs of the EMCCD. So far, only the CCD97 was tested, but the controller is technically able to drive other EMCCD as well as classical CCD. CCCP is designed to operate the vertical and horizontal clocks at the maximum speed specified by manufacturers, allowing fast read-out. Lower speed operation is also possible, down to a few kHz of pixel rate.

REFERENCES

[1] M. Marcelin, “3D-NTT: A New Instrument for the NTT Based on Versatile Tunable Filter Technology,” in Science Perspectives for 3D Spectroscopy, M. Kissler-Patig, J. R. Walsh, and M. M. Roth, eds., pp. 15–+, 2007.
[2] P. Jerram, R. Pool, R. Bell, D. Burt, S. Bowring, M. Spencer, M. Hazelwood, I. Moody, N. Catlett, and P. Heyes, “The l3ccd: Low-light imaging without the need for an intensifier,” tech. rep., E2V Technologies, http://www. e2vtechnologies.com/secure/datasheets/l3vision_ccds/4306a_20.pdf, 2001.
[3] M. Stanford and B. Hadwen, “The noise performance of electron multiplying charge coupled devices,” IEEE Transactions on Electron Devices 1488R, Dec. 2002.
[4] E. Lantz, J.-L. Blanchet, L. Furfaro, and F. Devaux, “Multi-imaging and Bayesian estimation for photon counting with EMCCDs,” MNRAS 386, pp. 2262–2270, June 2008.
[5] A. G. Basden, C. A. Haniff, and C. D. Mackay, “Photon counting strategies with low-light-level CCDs,” MNRAS 345, pp. 985–991, Nov. 2003.
[6] S. Tulloch, “Monte Carlo Modeling of L3 Detectors in High Time Resolution Applications,” in High Time Resolution Astrophysics: The Universe at Sub-Second Timescales, American Institute of Physics Conference Series 984, pp. 148– 161, Feb. 2008.
[7] Y. Wen, B. J. Rauscher, R. G. Baker, M. C. Clampin, P. Fochie, S. R. Heap, G. Hilton, P. Jorden, D. Linder, B. Mott, P. Pool, A. Waczynski, and B. Woodgate, “Individual photon counting using e2v L3 CCDs for low background astronomical spectroscopy,” in High Energy, Optical, and Infrared Detectors for Astronomy II. Edited by Dorn, David A.; Holland, Andrew D.. Proceedings of the SPIE, Volume 6276, pp. 62761H (2006)., Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 6276, July 2006.
[8] “Low-light technical note 4 dark signal and clock-induced charge in l3vision ccd sensors,” tech. rep., E2V Technologies, http://www.e2vtechnologies.com/secure/datasheets/l3vision_ccds/ low_light_tn4.pdf, June 2004.
[9] J. R. Janesick, Scientific charge-coupled devices, Scientific charge-coupled devices, Bellingham, WA: SPIE Optical Engineering Press, 2001, xvi, 906 p. SPIE Press monograph, PM 83. ISBN 0819436984, 2001.
[10] J.-L. Gach, C. Guillaume, O. Boissin, and C. Cavadore, “First Results of an L3CCD in Photon Counting Mode,” in Scientific Detectors for Astronomy, The Beginning of a New Era, P. Amico, J. W. Beletic, and J. E. Belectic, eds., Astrophysics and Space Science Library 300, pp. 611–614, 2004.
[11] O. Daigle, J.-L. Gach, C. Guillaume, C. Carignan, P. Balard, and O. Boisin, “L3CCD results in pure photoncounting mode,” in Optical and Infrared Detectors for Astronomy., J. D. Garnett and J. W. Beletic, eds., pp. 219–227, Sept. 2004.
[12] C. Mackay, A. Basden, and M. Bridgeland, “Astronomical imaging with L3CCDs: detector performance and highspeed controller design,” in Optical and Infrared Detectors for Astronomy. Edited by James D. Garnett and James W. Beletic. Proceedings of the SPIE, Volume 5499, pp. 203-209 (2004)., J. D. Garnett and J. W. Beletic, eds., Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 5499, pp. 203–209, Sept. 2004.
[13] S. Tulloch, “L3 CCD Wavefront Sensor Developments at the ING,” in Scientific Detectors for Astronomy 2005, J. E. Beletic, J. W. Beletic, and P. Amico, eds., pp. 303–+, Mar. 2006.
[14] J.-L. Gach, O. Hernandez, J. Boulesteix, P. Amram, O. Boissin, C. Carignan, O. Garrido, M. Marcelin, G. Ostlin, ¨ H. Plana, and R. Rampazzo, “Fabry-P´erot Observations Using a New GaAs Photon-counting System,” PASP 114, pp. 1043–1050, Sept. 2002.
[15] O. Hernandez, K. Fathi, C. Carignan, J. Beckman, J. . Gach, P. Balard, P. Amram, J. Boulesteix, R. L. M. Corradi, M. de Denus-Baillargeon, B. Epinat, M. Rela˜no, S. Thibault, and P. Vall´ee, “GHαFaS : Galaxy H-alpha Fabry-Perot System for the WHT,” ArXiv e-prints 805, May 2008. Accepted for publication in PASP.

Télécharger le communiqué de presse
Contactez-nous pour bénéficier de notre expertise CCD & EMCCD
Nous joindre
Produits
Applications
Ressources
À propos de nous
crossmenuchevron-downarrow-right