# Measuring the electrical and optical properties of silicon photomultipliers

*Slawomir Piatek, PhD, Hamamatsu Corporation & New Jersey Institute of Technology*

*February 2, 2014*

## Introduction

The Multi-Pixel Photon Counter (also known as a silicon photomultiplier or SiPM) is a two-prong photodetector consisting of pixels connected in parallel; one pixel is a series combination of an avalanche photodiode (APD) and a quenching resistor, R_{Q} (see Figure 1).

In a typical operation, an MPPC is externally reverse-biased with a *V*_{BIAS} that is up to a few volts greater than the breakdown voltage (*V*_{BR}) of the APD. An overvoltage, Δ*V* ≡ *V*_{BIAS} - *V*_{BR}, is one of the most important adjustable parameters affecting the performance of the detector, controlling its optical and electrical characteristics such as the gain (M), photon detection efficiency (PDE), and the rates of dark counts, afterpulsing, and optical crosstalk.

In the absence of light and ignoring dark counts, there is no current flowing through a pixel; thus, the voltage on the APD equals *V*_{BIAS}. That APD is said to be in Geiger mode – a highly unstable state that can be disturbed even by a single photon, leading to the formation of a self-sustained avalanche. The role of *R*_{Q} in a pixel is to extinguish or quench the avalanche restoring the APD back to Geiger mode. The gain is proportional to the charge released in the avalanche, and its typical values are 10^{5} to 10^{6}. Such a high gain is comparable to that of a photomultiplier tube (PMT), making an MPPC its possible substitute in applications involving detection of very faint light. The advantages of an MPPC compared to a PMT are a smaller size, immunity to magnetic fields, and lower biasing voltage. The detector’s disadvantages are that its performance is sensitive to fluctuations of the ambient temperature and of the biasing voltage, the presence of high rates of thermal noise (~ 1 MHz), and of correlated noise: afterpulsing and crosstalk.

Using an MPPC as a light detector in scientific and clinical research requires that the user of the detector understands its optical and electrical behavior under varied environmental conditions, the types and rates of noise, and the limits on reliable light detection. Although manufacturers of MPPCs provide information on their products, an independent check of the detector characteristics, especially in critical situations where a scientific breakthrough hinges on the validity of the observations, is necessary.

This article is a literature survey of experimental techniques used to characterize optical and electrical properties of an MPPC. Section I describes the basic principles behind the operation of first an APD and then an MPPC. Section II presents and discusses experimental setups to measure the most relevant optical and electrical properties of an MPPC. Section III summarizes the main points of the article.

## Section I: Background

A light-sensitive APD is a PN junction normally operated in reverse bias close to its *V*_{BR}. Unlike a photodiode that has a gain *M* = 1, an APD can have *M* > 1 owing to multiplication of the photo-generated charge carrier in the avalanche section of the depletion region by impact ionization. The values of M depend on *V*_{BIAS} with respect to *V*_{BR}. The following two subsections describe in greater detail the operation of an APD when *V*_{BIAS} < *V*_{BR} and *V*_{BIAS} > *V*_{BR}, respectively, and an MPPC.

### Avalanche photodiode (*V*_{BIAS} < *V*_{BR})

A single pixel in an MPPC is an APD connected in series with a resistor. Figure 2 shows a simplified equivalent circuit of an APD operating at reverse bias that is less than *V*_{BR}>. From left to right, the components are current source, ideal diode, junction capacitance (*C*_{J}), junction resistance (*R*_{J}), and series resistance (*R*_{S}).

The current source outputs current *I* = *MI*_{PH} = *M*σ*P*, where *M* is the gain and *I*_{PH} is the photocurrent. The gain is a function of applied voltage, whereas *I*_{PH} = σ*P*, where σ=σ(λ) is the spectral photosensitivity and P is the power of the incident light. The gain mechanism is avalanche multiplication of charge carriers by impact ionization in the high field region of the APD. If *V*_{BIAS} < *V*_{BR}, the avalanche, once triggered, abates completely in finite time.

The behavior of an ideal diode is described by the Shockley equation *I*_{D} = *I*_{0}[exp(*V*_{D}/*V*_{T}) – 1]. In this equation, *I*_{D} is the current through the diode, *I*_{0} is the reverse saturation current, *V*_{D} is the voltage on the diode, and *V*_{T} ≡ *kT*/*e* is the thermal voltage, where, *k* is the Boltzmann constant (*k* = 8.6 x 10^{−5} *eV*/*K*), *T* is temperature, and *e* is the magnitude of the charge of the electron (*e* = 1.6×10^{−19} *C*). In a typical operation of an APD, *V*_{D} has a large negative value and, thus, *I*_{D} ≈ - *I*_{0} (the minus sign indicates that this is a reverse current). Note that *V*_{D} is not an external bias voltage; it is the voltage across the diode and, also, across the current source I, *C*_{J}, and *R*_{J}.

The capacitor *C*_{J} represents the capacitance of the depletion region of the APD. Because the width of the region depends on the voltage across it (which equals *V*_{D}), *C*_{J} = *C*_{J}(*V*_{D}). *C*_{J} affects the performance of an MPPC in an important way: it controls the shape of the output current pulse. Therefore, measuring its value as a function of *V*_{BIAS} constitutes a crucial electrical characterization of the device.

The resistor *R*_{J} has a value equal to the resistance of the depletion region. The value depends on the width of the region and on the concentration of charge carriers there. *R*_{J} is a complicated function of *V*_{D} because the latter controls both the width and the degree of impact ionization in the region. When the APD is in Geiger mode and develops an avalanche, *R*_{J} becomes negligible.

The resistance *R*_{S} equals the combined resistance of the undepleted regions in the APD. The value of *R*_{S} depends on *V*_{D} since the latter controls the total length of the undepleted regions. *R*_{S} is on the order of several hundred ohms at zero bias and decreases to tens of ohms around *V*_{BR}.

### Avalanche photodiode (*V*_{BIAS} > *V*_{BR})

When a *V*_{BIAS} applied to an APD exceeds its *V*_{BR}, the APD is said to be in Geiger mode. McIntyre (1961) and later Haitz (1964) developed an electrical model (see Figure 3) for an APD operating in this state. Geiger-mode APD is represented by *C*_{J}, the switch S, voltage source *V*_{BR}, and series resistance *R*_{S}. In the absence of light and ignoring dark events, the switch S is open and *C*_{J} is charged to *V*_{BIAS}. The voltage on the APD is Δ*V* = *V*_{BIAS} – *V*_{BR} above *V*_{BR}, and the APD is in a highly unstable state. Suppose that the APD absorbs a photon and the resulting charge carrier, either an electron or a hole, triggers an avalanche. At this instant, the switch *S* closes and *C*_{J} begins to discharge through *R*_{S}, which tends to lower the voltage across the APD. However, since the APD is biased by a constant voltage source *V*_{BIAS}, current begins to flow through the terminals of the APD reaching a steady state value of Δ*V*/*R*_{S}. This current will persist unless *V*_{BIAS} is reduced to, or below, *V*_{BR}; doing so restores the APD to the light-sensitive state.

Suppose that a resistor *R*_{Q} >> *R*_{S} is now connected in series with the Geiger-mode APD, as shown in Figure 4.

In the absence of light and ignoring dark events, the switch *S* is open and *C*_{J} is charged to *V*_{BIAS} – no current flows through the APD. The APD is in a photosensitive "ready" state. The switch S closes the instant a photo-generated charge carrier triggers an avalanche. *C*_{J} begins to discharge through *R*_{S}; this process tends to lower the voltage across the APD and causes a voltage drop across *R*_{Q}. The latter implies that a current flows through the terminals of the APD, and Figure 5 depicts its time dependence.

The avalanche begins at *t* = *t*_{i}. The leading edge of the current pulse increases with time constant *R*_{S}*C*_{J} and reaches a maximum value of *i*_{max} = (*V*_{BIAS} – *V*_{BR})/(*R*_{Q} + *R*_{S}) ≈ Δ*V*/*R*_{Q}at *t*_{max}. At around this time the avalanche stops or is quenched because of *R*_{Q}. After *t*_{max}, *C*_{J} recharges to the nominal voltage of *V*_{BIAS} while the current decreases with a time constant *R*_{Q}*C*_{J}. Because *R*_{Q} >> *R*_{S} the leading edge of the pulse is steeper than the trailing (or declining) edge; the pulse is asymmetric around *t* = *t*_{max}. For a more detailed discussion of the electrical behavior of Geiger-mode APD see Haitz (1964).

MPPCs are designed so that *i*_{max} (referred to as the "latch" current) is about 20 μ*A*. The condition for proper quenching is 0 < Δ*V* < *R*_{Q}*i*_{max}. For Δ*V* = 1 *V*, the minimum *R*_{Q} is 50 kΩ. If *R _{Q}* is set below this limit, the quenching time increases, which would show as flattening of the current pulse around

*t*

_{max}. If

*R*

_{Q}is substantially below the minimum limit, the quenching never occurs and the current becomes constant with an amplitude of (

*V*

_{BIAS}–

*V*

_{BR})/(

*R*

_{Q}+

*R*

_{S}).

Integrating the current pulse with respect to time yields the total charge *Q* that transferred between the terminals of the APD. Since the transfer was triggered by a single charge carrier, the gain *M* = *Q*/*e*. A detailed analysis shows that *Q* = *C*_{J}Δ*V* and, thus, *M = C_{J}ΔV/e.*

*
*### Measuring

### MPPC

The MPPC is a rectangular array of light-sensitive elements – pixels – that are all connected in parallel and externally biased by a single voltage source, *V*_{BIAS}. Each pixel is a series combination of an APD and a resistor *R*_{Q}. By design, all pixels are identical. The pixels operate in Geiger mode, that is, *V*_{BIAS} > *V*_{BR}, or Δ*V* > 0.

Figure 6 shows n equivalent circuit of an MPPC; the dotted rectangle delineates a single pixel. In the absence of light and ignoring dark counts, all of the switches are in the OFF position, and the voltage is *V*_{BIAS} on each APD and is zero on each *R*_{Q}. Suppose that a single APD absorbs a photon and the process triggers an avalanche. At this instant, the switch in the equivalent circuit of that APD goes to the ON position, and the current pulse begins to flow through the terminals of the MPPC as shown in Figure 5. *R*_{Q} of that pixel quenches the avalanche, and the pixel is restored back to the "ready" state. If two (or more) photons simultaneously trigger avalanches in two (or more) distinct pixels, the current pulse flowing through the terminals of the MPPC is a superposition of the current pulses generated by the two (or more) pixels. In contrast, if a single pixel absorbs simultaneously two or more photons, the ensuing current pulse is identical to the one produced by a single photon.

When a pixel recovers from the avalanche, its *C*_{J} recharges to the nominal voltage *V*_{BIAS}. At the outset of the recharge process, the current flows from *V*_{BIAS} but also from the junction capacitors of all other pixels. The latter current keeps lowering Δ*V* on the other pixels until the value equals to the Δ*V* on the recovering pixel. At this instant, *V*_{BIAS} becomes the only source of current, and all of the pixels recharge together back to *V*_{BIAS}.

## Section II: Characterization of an MPPC

This section describes experimental methods of measuring electrical and optical characteristics of an MPPC.

### Measuring *R*_{Q} and *C*_{J}

The forward-bias I-V characteristic of an MPPC (e.g., see Figure 10) can be used to measure the approximate value of *R_{Q}. The I-V characteristic can be obtained with a circuit shown in Figure 7. Assuming that R_{Q} is the same for each pixel and that R_{Q} >> R_{S}, the slope of the linear region in the first quadrant equals N/R_{Q}, where N is the number of pixels. Thus, R_{Q} = N/slope.*

*
$Z=\frac{1}{N}\left({R}_{Q}+{R}_{S}+\frac{{\displaystyle \frac{{R}_{j}}{jw{C}_{j}}}}{{R}_{J}+{\displaystyle \frac{1}{jw{C}_{j}}}}\right)=\frac{{R}_{Q}+{R}_{s}}{N}+\frac{1}{{\displaystyle \frac{N}{{R}_{J}}}+j\omega N{C}_{J}}$
*

Figure 8 shows an equivalent circuit of an MPPC that illustrates how *C*_{J} of a single pixel can be measured. There are N pixels connected in parallel, where one pixel is represented by *C*_{J} in parallel with *R*_{J} and this combination is in series with *R*_{S} + *R*_{Q}. The terminals of the MPPC are at the right.

The impedance Z of the network as seen from the terminals of the MPPC is given by the following equation.

Equation 1 implies that the effective capacitance seen at the terminals of the MPPC is *NC*_{J}, the capacitance of a single pixel multiplied by the number of pixels. Because *C*_{J} is a function of reverse voltage, the Keithley C-V 590 analyzer is well suited for determining this dependence. Figure 9 shows the block diagram of the experimental setup.

The MPPC is inside a temperature-controlled lightproof Faraday shield. The terminals of the detector are connected to the input and output of the C-V 590 using BNC cables. The "ground" of the cables is connected to the shield. The polarity of the connection is such that when the C-V 590 performs the voltage sweep, the MPPC is reverse-biased.

Figure 10 is an example of forward and reverse I-V and C-V characteristics for an MPPC with 667 pixels. The forward I-V characteristic becomes linear at a bias voltage of around 0.7 V. The inverse of the slope in the linear region is about 225 Ω implying that *R*_{Q} ≈ 150 kΩ. The reverse-bias I-V characteristic shows that the breakdown voltage *V*_{BR} ≈ 70 *V* (reverse).

The combined junction capacitance *N* x *C*_{J} decreases steeply when the reverse voltage increases from 5 *V* to 20 *V*. Above 20 *V*, the decrease is much less pronounced; at around *V*_{BR}, *N* x *C*_{J} ≈ 60 *pF*, or *C*_{J} ≈ 90 *fF*.

*
$\xi \left(\u2206t\right)=1-exp\left(-\frac{\u2206t}{{T}_{r}}\right)$
**
*### Observing noise

An MPPC generates three noise components: dark noise (or dark counts), afterpulses, and optical crosstalk. Thermally-generated charge carriers that trigger an avalanche are responsible for dark noise. Dark noise is random and is always present regardless of the light level illuminating the MPPC. Its rate (number of pulses per unit time) increases with temperature and Δ*V*. The presence of the other two noises is correlated to the occurrence of an avalanche, which has been triggered either by a photo-generated or thermally-generated (dark) charge carrier. A release of a trapped charge in a pixel experiencing an avalanche can trigger a secondary avalanche while the pixel is recovering from the primary avalanche. This is afterpulsing. Optical crosstalk arises when an avalanche in a pixel emits photons that trigger one or more avalanches in the adjoining pixels. The resulting output signal is a superposition of the signals from the "fired" pixels.

Figure 11 illustrates an experimental setup to observe and record the three noise components. An MPPC is inside a temperature-controlled dark box. The variable high voltage source (e.g., Keithley 487) reverse-biases the MPPC and sets the value of Δ*V*. The output signal from the MPPC is AC coupled to a "fast" amplifier (e.g., Hamamatsu C5594) using a 10 *k*Ω load resistor and 100 nF coupling capacitor. The output of the amplifier connects to a digital oscilloscope (e.g., LeCroy Wave Pro 7300A). The values of the load resistor and the coupling capacitor can be adjusted to optimize the fidelity of the displayed signals.

Figure 12 is an example of an oscilloscope display showing an accumulation of waveforms (the persistence was set to 10 s) resulting from the three types of noise. The most numerous waveforms are due to random dark noise. Such a waveform is indistinguishable from the one that would have been produced had the avalanche been triggered by a photo-generated charge carrier. A waveform due to afterpulsing trails (in time) the primary waveform – a waveform due to dark noise. The delay is on the order of tens of ns. The signals due to optical crosstalk are practically in phase with the primary signal and, thus, the resulting waveform is the superposition of the signals. Those that are labeled 2 p.e. (3 p.e.) result when one (two) neighboring pixel is (are) affected.

*
*

### Measuring crosstalk

In the absence of light, 2 p.e. (or greater) waveforms originating from an MPPC can be either due to random dark noise occurring in two or more pixels simultaneously or crosstalk. The probability of the former is very small and will be neglected. Therefore, counting the number of waveforms above the 1 p.e. threshold per unit time yields the crosstalk rate, while dividing this rate by the rate of 1 p.e. waveforms yields the crosstalk probability.

Figure 13 is an experimental setup to measure the rates (and probabilities) of crosstalk and afterpulsing. The MPPC is inside a temperature-controlled and light-leak-proof dark box reversed-biased by a high voltage source. The amplifier amplifies the output signal from the MPPC and feeds the resulting signal to the discriminator. If the input crosses a user-specified threshold level, the discriminator issues a logical output pulse of some duration (several ns), which the Scaler module, with its gate set to 1 s, counts. Figure 14 shows an example of a resulting count rate as a function of the threshold level for the Hamamatsu S10362-11-050C MPPC operating at *ΔV = 1.3 V*.

Figure 14 shows that the noise rate, ν, decreases with the discriminator threshold level showing a characteristic step-like dependence. The horizontal red dashed lines delineate the values of ν for the threshold levels that correspond to 0.5 p.e., 1.5 p.e., and 2.5 p.e., respectively. Setting the threshold level to 0.5 p.e., which is above the electronic noise but below 1 p.e., gives the rate of the total dark noise ν_{0.5 p.e.}. Setting the threshold to 1.5 p.e. gives the rate of the events due to crosstalk ν_{1.5 p.e.}. The probability of crosstalk P_{C} is therefore given by *P*_{C} = ν_{0.5 p.e.}/ ν_{1.5 p.e.}.

*
The crosstalk probability P_{C} depends on the overvoltage ΔV, and the experimental setup depicted in Figure 13 can also be used to determine this dependence. The procedure is to measure P_{C} as outlined above for different values of ΔV, which is controlled by the high voltage source. Figure 15 is an example of the resulting dependence (note that the gain and ΔV are linearly proportional to each other). The figure shows that P_{C}* steeply increases with Δ

*V*for all of the tested MPPCs.

### Measuring afterpulsing

When a pixel develops an avalanche, some charge carriers can be trapped by lattice defects. Their subsequent release, while the pixel recovers from the primary avalanche, can trigger a secondary correlated avalanche – an afterpulse. The amount of charge released in an afterpulse depends on the time delay δ*t* between the primary and secondary avalanches, or the recovery state of the pixel. The ratio of charge released in an afterpulse and that released in the primary avalanche can be modeled by the following equation.

In this equation, τ_{r} is the recovery time (from an avalanche) of a pixel. If *δ t is small compared to τ_{r}, then ξ is small implying that an afterpulse has amplitude less than 1 p.e. The amplitude of an afterpulse increases as δt increases and in the limit δt >> τ_{r} it is 1 p.e. In this limit, the afterpulse is indistinguishable from a signal that was generated either thermally or optically.*

*
${\eta}_{tp}\left(\u2206t\right)=\frac{{N}_{tp}}{{{\rm T}}_{tp}}exp\left(-\frac{\u2206t}{{{\rm T}}_{tp}}\right)$
${\eta}_{\alpha p}\left(\u2206t\right)=\frac{{N}_{\alpha pf}}{{{\rm T}}_{tpf}}exp\left(-\frac{\u2206t}{{{\rm T}}_{\alpha pf}}\right)+\frac{{N}_{\alpha ps}}{{T}_{\alpha ps}}exp\left(-\frac{\u2206t}{{{\rm T}}_{\alpha ps}}\right)$
${P}_{\alpha p}=\frac{{\int}_{\mathrm{O}}^{\infty}\xi {\eta}_{\alpha \rho}d\u2206t}{{\int}_{\mathrm{O}}^{\infty}\xi \left({\eta}_{\alpha \rho}+{\eta}_{t\rho}\right)d\u2206t}$
*

An experimental setup depicted in Figure 13 can be used to measure the probability of afterpulsing *P*_{ap}. Here, the Scaler is replaced with Time to Digital Converter (TDC), which measures δ*t* for consecutive pulses coming from the discriminator (threshold level at 0.5 p.e.). The solid black curve in Figure 16 is an example of the measured dark noise time-difference distribution. For δ*t* < 20 ns, the distribution is unreliable because to detect an afterpulse, the primary signal must fall below the discriminator threshold level; the time to do so is on the order of 10 – 15 ns (see Figure 12). For δ*t* greater than the pixel recovery time, the dark-noise time-difference distribution can be modeled as the sum of two probability density functions (what follows is from Eckert et al. 2010).

Equation 3 is a model probability density function for events due to thermal pulses. Here, *N*_{tp} is the integrated number of thermal events and 1/τ_{tp} is the rate (or frequency) of these events.

Equation 4 is a model probability distribution function for events due to afterpulsing. The function is a sum of two exponentials: the first exponential term models "fast" events with a characteristic time τ_{apf}, whereas the second models "slow" events with a characteristic time τ_{aps}. The parameters *N*_{apf} and *N*_{aps} are the integrated number of "fast" and "slow" events, respectively.

The probability density functions given by Equation 3 and Equation 4 are simultaneously fitted to the data, as shown in Figure 16. The fit yields the values of the adjustable parameters in the probability density functions.

With the adjustable parameters determined from the fit, Equation 5 can be used to calculate the probability of afterpulsing. Equation 2 gives the charge ratio function ξ. To evaluate the integral it is necessary to adopt the appropriate pixel recovery time τ_{r}.

The probability of afterpulses *P*_{ap} is a function of overvoltage Δ*V* for two reasons. As Δ*V* increases: (1) the amount of charge in a primary avalanche increases making charge trapping more likely and (2) the probability that a charge carrier triggers an avalanche in the depletion region increases. The same experimental setup and analysis as described above can be used to determine how *P*_{ap} depends on Δ*V* by varying the bias voltage on the MPPC. Figure 17 is an example of such dependence for three MPPCs manufactured by Hamamatsu.

*
$P\left(\eta ,{\eta}_{pe}\right)=\frac{\eta {}_{pe}{}^{\eta}e^{-{n}_{pe}}}{\eta !}$
${\eta}_{pe}=-1n\frac{{{\rm N}}_{0}}{{{\rm N}}_{{\rm T}}}+1n\frac{{{\rm N}}_{0}^{D}}{{{\rm N}}_{{}_{{\rm T}}}^{D}}$
$PDE=\frac{\eta \rho e}{{\eta}_{\u22ce}}$
$M=\frac{QD{C}_{cal}}{e{K}_{amp}}{\u2206}_{pp}$
${N}_{f}=N\left(1-e\frac{-{N}_{\gamma \xb7PDE}}{N}\right)$
*

### Measuring photon detection efficiency

The photon detection efficiency, PDE, of an MPPC is defined as the product of quantum efficiency η, the geometrical fill factor *F*_{g}, and a combined probability *P*_{a} of electrons and holes initiating Geiger-mode avalanche; that is, *PDE* = η*F*_{g}P_{a}. The value of η is a probability that a photon generates a charge carrier that enters the avalanche region in a pixel. It is a function of the wavelength of light λ and temperature T, but it is independent of Δ*V* for full depletion. The geometrical fill factor *F*_{g} is a ratio of the combined photosensitive area to the total area of an MPPC. Due to edge effects, *F*_{g} exhibits a weak dependence on Δ*V*. The avalanche probability *P*_{a} is a function of λ, *T*, and Δ*V*.

The PDE is a crucial parameter in photon counting experiments; therefore, its value must be known for a given wavelength and overvoltage if absolute photon fluxes are to be measured.

### PDE versus overvoltage

Figure 18 shows an experimental setup to measure the PDE of an MPPC as a function of Δ*V* using a pulsed-LED method discussed in Musienko et al. (2006) and Vacheret et al. (2011).

An LED with a narrow emission spectrum emits pulses of light with a duration and repetition rate set by the pulse generator. The neutral density filter controls the number of photons per pulse, and the collimator allows a uniform illumination of the photosensitive section of an MPPC. The optical depth of the filter is adjusted so that 2 to 5 pixels are activated per pulse. The number of photons per pulse *N*γ arriving at the plane of the MPPC must be measured using a calibrated photodetector, such as the XP2020 or R7899 PMTs. The high-voltage source, such as the Keithley 487, reverse-biases an MPPC and determines the value of Δ*V*. A fast transimpedance amplifier amplifies the output signal, and the picoscope, such as the PicoScope 5023, digitizes and integrates the output to obtain a charge amplitude spectrum such as the one shown in Figure 19.

The probability *P* of detecting *n* photons by an MPPC is given by the Poisson distribution.

In this equation, *n*_{pe} is the average number of photoelectrons. The events for *n* = 0 (no photon detected) correspond to the "pedestal" on the charge amplitude spectrum; the subsequent peaks are for *n* =1, 2, … (see Figure 19). The values of *P*(*n*, *n*_{pe}) for *n* > 0 are distorted by the presence of crosstalk and afterpulsing, but the value for *n* = 0 is not. Measuring *P*(0, *n*_{np}), which is the ratio of the number of events in the pedestal *N*_{0} and the total number of events *N*_{T}, allows unbiased determination of *n*_{pe}. Both *N*_{0} and *N*_{T} are determined from the charge amplitude spectrum. Since *P*(0, *n*_{pe}) = *N*_{0}/*N*_{T} = exp(-*n*_{pe}), then *n*_{pe} = -l*n*(*N*_{0}/*N*_{T}). Taking into account dark pulses, *n*_{pe} is given by the following.

Here, *N*_{0}^{D} and *N*_{T}^{D} are the number of events in the pedestal and the total number of events, respectively, in a charge amplitude spectrum taken under dark conditions (see Figure 20). Thus, the PDE is given by Equation 8.

### Measuring gain and breakdown voltage

A charge amplitude spectrum of an MPPC, such as the one depicted in Figure 19, can be used to determine the gain *M* of the device using the following equation.

In the equation, *QDC*_{cal} is the charge per one channel of *QDC* (charge to digital converter), δ_{pp} is the number of channels separating two adjacent peaks, *K*_{amp} is the amplification factor of the electronics, and e is the value of the fundamental charge. If charge amplitude spectra are collected for various settings of *V*_{BIAS}, *M* as a function of *V*_{BIAS} can be determined, as shown in Figure 21.

As expected, the relationship between *M and V _{BIAS}* is linear with the slope and vertical-intercept equal to

*C*

_{J}/

*e*and

*C*

_{J}

*V*

_{BR}/

*e*, respectively. Fitting a straight line to the data yields the values of

*C*

_{J}and

*V*

_{BR}.

### Measuring time jitter

Suppose that a photodetector is illuminated with δ-function pulses of light that are equally spaced in time with period *T* and, in response, the detector produces output signal pulses. If the time interval between the consecutive output pulses varies, the output does not have the same period *T* and the detector is said to exhibit "time jitter."

Figure 22 shows a possible experimental setup (Nagano et al. 2012) to measure the time jitter of an MPPC. The detector is inside a temperature-controlled dark box. A pico-second laser (λ = 650 nm, 70 ps pulse width, 1 kHz repetition rate, and time jitter less than 10 ps) is the source of light. The beam splitter divides the beam into two beams. One illuminates a photodiode and the other - after passing through a neutral density filter and collimator - the MPPC. The filter and collimator combination can be adjusted to vary the number of photons per pulse striking the MPPC. The output from the photodiode, after passing through a constant fraction discriminator (CFD), provides the "start" pulse to the time-to-amplitude (TAC) converter. The role of the CFD is to filter out noise: only those signals for which the value of the leading edge exceeds some fraction (e.g., 30%) of the difference between the peak value and the base level are permitted. After being amplified by a factor of 30 to 50, the output from the MPPC passes through the CFD and delay circuit (or line) and arrives at the TAC providing the "stop" pulse. The TAC outputs a signal whose amplitude is proportional to the time difference between the start and stop pulses. The multi-channel analyzer (MCA) stores and sorts the pulses according to their height, and displays their histogram, or spectrum, such as that displayed in Figure 23.

For a normal distribution, the full-width at half maximum (FWHM) is a measure of the total time jitter: *FWHM* = 2.355σ_{T}. If the time jitters of the individual components of the setup are independent of one another, the sum of their squares equals σ_{T}^{2}. From the time jitters of the laser pulses, photodiode, and the combined electronics, the time jitter of the MPPC can be calculated.

### Measuring linearity and dynamic range

An MPPC consists of N pixels where each pixel can exist in only two states: ON or OFF. A pixel in the OFF state switches to the ON state if it absorbs one or more photons and the resulting photo-generated charge carriers trigger an avalanche. The maximum possible output signal from an MPPC is N p.e. in response to a δ-pulse of light containing *N*_{γ} photons. Thus, the total number of pixels determines the dynamic range of the detector.

Equation 10 is a model that gives the number of pixels that have turned ON (or that have "fired") in response to *N*_{γ} of incident photons on the surface of an MPPC with N pixels and the photon detection efficiency PDE.

If the product of *N*_{γ} and PDE is much smaller than N, Equation 10 reduces to *N*_{f} = *N*_{γ}*PDE*. Here, the response of the detector is linear with *PDE* being the constant of proportionality between *N*_{f} and *N*_{γ}. However, in the limit *N*_{γ}*PDE* → ∞, *N*_{f} → *N*, the saturation limit.

Gruber at al. (2014) have investigated how a variety of MPPCs respond to pulses of light (32 ps FWHM duration and λ = 404 nm) containing from one to thousands of photons using a setup depicted in Figure 24. The PIN photodiode can be calibrated at low light levels (where *N*_{γ}*PDE* << *N*) so that its output is proportional to *N*_{seed} ≡ *N*_{γ}*PDE* at any light level. Here, *N*_{seed} is the average number of photons incident on the photosensitive area of an MPPC. Thus, the output from the PIN photodiode yields *N*_{seed} and the scope yields the height of the output signal, which is proportional to *N*_{f}. The setup, thus, allows a measurement of *N*_{f} as a function of *N*_{seed} for different settings of overvoltage and of different MPPCs. The results can be normalized to allow a direct comparison among the MPPCs used in the tests and comparison to the predictions from Equation 10.

*
*

The results are surprising. Figure 25 shows the response curves of four different MPPCs normalized to the total number of pixels in each MPPC for high light levels (left panel) and medium to low light levels (right panel). The figure shows that the outputs from the MPPCs are fairly linear and consistent with Equation 10 up to *N*_{seed}/*N* ≈ 0.5. Above this limit, the outputs are larger than the predictions from Equation 10 (Equation 3 in the original study) and the saturation level may be more than 2N (Hamamatsu 050U MPPC), in disagreement with Equation 10. The study finds no acceptable explanation for this "oversaturation" and questions out current understanding of how the device operates.

*
*

## Section III: Summary

This article is a review of the experimental methods to characterize optical and electrical properties of an MPPC. The following are the main points of the article:

- An MPPC is an array of pixels connected in parallel where one pixel is a series combination of an APD and a quenching resistor. The APD is externally biased to operate in Geiger mode.
- The junction capacitance and quenching resistance are key electrical parameters of an MPPC. Both can be measured from forward and reverse I-V characteristics of the device.
- Three major sources of noise exist in an MPPC: thermally-induced random dark counts, correlated optical crosstalk, and afterpulsing. The rates of these noises as a function of temperature and overvoltage can be measured, and these rates will enter into the calculations of the photon flux on the detector.
- Photon detection efficiency is a ratio of the detected number of photoelectrons to the number of photons in a pulse of light. Knowing how this parameter depends on the wavelength of light and overvoltage is crucial in all photon counting experiments. This article describes a method of measuring this parameter, correcting for the presence of noise.
- The time interval between the instant when a photon strikes a pixel of an MPPC and the instant the device produces an output varies in a sequence of events. The rms of the variation is known as time jitter - a crucial characteristic in experiments involving time-of-flight or coincidence. Time jitter is a function of the wavelength of light, overvoltage, and temperature. This article describes a method of measuring this parameter.
- The total number of pixels of an MPPC limits its dynamic range.
- The output of an MPPC is linear in a regime where the product of the number of the incident photons and the photo detection efficiency is much smaller than the total number of pixels. MPPCs exhibit a phenomenon of "oversaturation" where the output signal may exceed N p.e. (N is the number of pixels) in response to a pulse of light containing N
_{γ}of photons. This phenomenon is not well understood.

**References**

- P. Eckert, H.C. Schultz-Coulon, W. Shen, R. Stamen and A. Tadday, "Characterization studies of silicon photomultipliers,"
*Nucl. Instrum. and Meth.*A 620 (2010): 217-26. - R. H. Haitz, "Model for the Electrical Behavior of a Microplasma,"
*J. of Appl. Phys.,*vol. 35, no. 2 (1964): 1370-1376. - L. Gruber, S. E. Brunner, J. Marton, K. Suzuki, "Over saturation behavior of SiPMs at high exposure,"
*Nucl. Instrum. Meth.*A 737 (2014): 11-18. - R. J. McIntyre, "Theory of Microplasma Instability in Silicon"
*J. of Appl. Phys.,*vol. 32, no. 6 (1961): 983-995. - T. Murase, H. Oide, H. Otono, S. Yamashita, "Development of PPD: characterization and simulation," International Workshop on New Photon Detectors PD09, June 24-26 2009, Shinshu University, Matsumoto, Japan.
- Y. Musienko, S. Reucroft and J. Swain, "The gain, photon detection efficiency and excess noise factor of multi-pixel Geiger-mode avalanche photodiodes,"
*Nucl. Instrum. Meth.*A 567 (2006): 57. - T. Nagano, K. Sato, A. Ishida, T. Baba, R. Tsuchiya, K. Yamamoto, "Timing Resolution Improvement of MPPC for TOF-PET Imaging," 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference Record (NSS/MIC).
- A. Vacheret, G. J. Barker, M. Dziewiecki, P. Guzowski, M. D. Haigh, B. Hartfiel, A. Izmaylov, W. Johnston, M. Khabibullin, A. Khotjantsev, Y. Kudenko, R. Kurjata, T. Kutter, T. Lindner, P. Masliah, J. Marzec, O. Mineev, Y. Musienko, S. Oser, F. Retiere, R. O. Salih, A. Shaikhiev, L. F. Thompson, M. A. Ward, R. J. Wilson, N. Yershov, K. Zaremba and M. Ziembicki, "Characterization and simulation of the response of multi-pixel photon counters to low light levels,"
*Nucl. Instrum. Methods Phys. Res.*A vol. 656, no. 1 (2011): 69-83.