Open access peer-reviewed chapter

Avalanche Photodiode Focal Plane Arrays and Their Application to Laser Detection and Ranging

Written By

Hai-Zhi Song

Submitted: April 10th, 2018 Reviewed: September 4th, 2018 Published: November 5th, 2018

DOI: 10.5772/intechopen.81294

Chapter metrics overview

1,914 Chapter Downloads

View Full Metrics


Focal-plane avalanche photodiodes (APDs) are being more and more widely and deeply studied to satisfy the requirement in weak light and single photon imaging. The progresses of this worldwide study, especially the distinctive researches and achievements in Southwest Institute of Technical Physics and University of Electronic Science and Technology of China are reviewed in this chapter. We successfully fabricated up to 64 × 1 linear-mode Si APD arrays, and 32 × 32–64 × 64 Si single-photon avalanche detector (SPAD) arrays, and applied them in Laser Detection and Ranging (LADAR) platforms like driverless vehicles. Also, we developed 32 × 32–64 × 64 InGaAsP/InP SPAD arrays, and constructed three-dimensional imaging LADAR using them. Together with the progresses of other groups and other materials, we see a prospective future for the development and application of focal-plane APDs.


  • avalanche photodiode
  • focus plane
  • laser detection and ranging

1. Introduction

Avalanche photodiodes (APDs) have been widely studied and effectively applied in commercial, military, and academic fields [1] for a few decades. Compared with p-i-n photodiodes, APDs provide higher gain, higher sensitivity and lower detection limit [2], so they are mostly well applied in optical communications [3], imaging [4, 5], and single photon detection [6, 7] in recent years. As all-solid-state optoelectronic devices operating at room-temperature or under thermoelectrically-cooled conditions, APDs are scalable to numerous pixels so that they are taking more and more important roles in focal-plane processing and imaging [8]. Owing to the advantages such as internal photoelectric gain, small size, low driving voltages, high efficiency, and fast response, focal-plane APD arrays bring about new three-dimensional (3D) imaging techniques which provide much wealthier and more accurate information for object recognition and identification [9]. Advanced 3D imaging technologies are strongly required in radar systems including laser detection and ranging (LADAR), so the focal-plane APDs and their LADAR applications were widely and deeply studied in recent years [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. For the purpose of more progress in the future, it is necessary to take an overview on the present research and production of APD arrays. Briefly, the most significant progress is made by MIT Lincoln Laboratory. They developed state-of-the art products of Si and InP/InGaAs Geiger-mode focal-plane arrays [10], which have been successfully applied in a few LADAR systems. Princeton Lightwave also succeeded in producing focal-plane single photon avalanche detector (SPAD) arrays and commercializing their single-photon camera based on the SPADs [11]. The research and production of other groups [12, 13, 14] may also be valuable as references for future developments. In this chapter, we review the research and application of the focal-plane APDs in Southwest Institute of Technical Physics and University of Electronic Science and Technology of China [15, 16, 17, 18, 19, 20]. It includes linear mode Si APD arrays, Si SPAD arrays and InGaAsP/InP SPAD arrays, which have been applied in LADAR systems.


2. Linear-mode Si APDs

The detection of weak light is technically significant in many application fields such as single molecule fluorescence, high-speed quantum cryptography, and infrared detection [21, 22, 23]. In all the application fields, APD devices are strongly required to perform photon-counting with high quantum efficiencies, quick optoelectronic response, and low dark counting rates (noise). LADAR imaging systems work in the way of sampling the spatial and/or temporal information of the optical radiation to an array of detectors. Linear-mode (applied bias slightly lower than the breakdown voltage) APDs are often desired by LADAR systems because their dead-time is normally much shorter than that of Geiger-mode (applied bias slightly higher than the breakdown voltage) APDs so that they can measure sequential pulse returns from closely spaced multiple objects. In extreme cases, linear-mode APDs can even detect a few photons or a single photon, which adds an extra dimension to LADAR scene data [21]. Generally, in the near-infrared spectral band, especially at 905 nm, Si APDs might be applied for ultra-weak light detection, and can be used in linear-mode at gains up to about 500 or greater [23]. Therefore, linear-mode Si APD arrays were developed and applied in LADAR systems.

2.1. Fabrication of the linear-mode Si APD chips

A basic linear-mode APD detector, as shown schematically in Figure 1, consists of the APD element and the readout integrated circuit (ROIC) [24, 25]. The ROIC is composed of a trans-impedance amplifier (TIA), a stabilivolt source circuit and a comparer. The APD element converts incident light signal into primary photo-generated carriers and photocurrent, then amplifies the resulting photocurrent through internal avalanche gain, i.e. the impact ionization. The TIA converts the amplified APD current into a voltage signal, which is proportional to the total multiplied charge delivered by the APD.

Figure 1.

Schematic of the linear-mode Si APD.

The design and simulation of an APD device were carried out using a full-band Monte Carlo (MC) device model. For each APD geometry, the MC model incorporates realistic band structures [26, 27]. The basic reach-through APD model with separate layers of absorption, charge and multiplication (SACM) is shown in Figure 2(a). In particular, it is important to know that, in general, electrons can be much more ionized than holes in silicon. Electrons rather than holes should be swept by the electric field into the high field region where the multiplication takes place. Thus, there should be a π-type absorbing region of suitable width for absorbing the incident radiation, and the radiation should be able to enter this region with no loss in any n+-type layer.

Figure 2.

(a) Schematic cross section (not to scale), and (b) schematic profile of electric field of a typical APD structure.

The basic design of a reach-through APD consists of a narrow high-field region where the multiplication takes place, with a much wider low field region in which the incoming radiation is absorbed. As schematically shown in Figure 2(b),an avalanche process occurs as the electric field in a p-n junction is higher than the critical field (Ecr) at which impact ionization of carrier starts. The electric field in the p-n junction of a Si APD should be some 2–5 × 105 V/cm. It should not be more than 106 V/cm at which the Zener effect may take place [28, 29, 30].

For satisfactory operation of the APD, the high resistivity π-type substrate must be fully depleted by the applied bias voltage. Generally, it works well provided the substrate wafer is not too thick and the required response times are not less than ~10 ns. However, fabrication of Si APDs on 6-inch or 8-inch wafers, as is now usually the case, will often mean that the absorption region is thick (~700 μm), operating voltages are high, and response times are slow. These problems can be avoided with the use of an epitaxial version of the design in Figure 2(a). In this approach, a high-resistivity π-type layer is epitaxially grown on the surface of a low resistivity p-type Si substrate. The absorption region (epitaxial layer) may be of any thickness (typically chosen to be in the range of 30–50 μm), and its resistivity is chosen to be low enough so that it does not introduce a significant series resistance. When bias voltage is applied, the depletion layer stops at the interface between the substrate and the epitaxial layer [31, 32, 33, 34]. While fast response is a requirement, the narrow active region of this APD is normally the best option.

Linear-mode Si APD arrays are fabricated by adopting Si planar manufacturing process on a high-resistivity π-type layer, which is grown epitaxially on the top of a low resistivity p-type Si substrate. The initial material developed is the Si layer 35–40 μm thick, a highly-resistive epitaxial layer on a p+-type Si <111> substrate.

2.2. Design of the TIA

As mentioned above, the ROIC chips consist of a voltage-stabilized source and a TIA. The voltage-stabilized source effectively reduces external noise jamming and increases voltage suppression ratio of the power source. Here we have the structure of a new-type regulated cascode circuit configuration which is compatible with the APD chips. The bandwidth, the parallel negative feedback and the trans-impedance gain of the TIA are improved by using regulated cascode circuit [35, 36].

Figure 3 schematically illustrates the TIA with the regulated cascode circuit configuration. The regulated cascode circuit consists of common gate amplifier input stage (including R1, R2and M2) and common source amplifier stage (including M1 and R3), forming partial current parallel negative feedback. The primary photocurrent generated by APD is imported to the source electrode of the cathode-input amplifier. Then the APD’s current signal is converted into a voltage signal. The source follower, consisting of M3 and M4, is used for isolation. The secondary common-source amplifier is composed of M3 and R4, which play a part role in amplifying signal again. To improve the output drive capability, the output stage contains two-stage source followers, made up of M6 + R5 and M7 + M8 respectively.

Figure 3.

The TIA with the regulated cascode circuit configuration.

2.3. Properties of the Si APD array

Developed at SITP, Si APD arrays were characterized at UESTC. The fabricated devices exhibit high primary photoelectric sensitivity (about 0.5 A/W @905 nm at gain M = 1) and high speed of operation (about 10 ns). Figure 4 shows an example of typical dependences of the gain on the reverse bias. As the bias arises up to the reach-through voltage Vrt, it depletes the π-type avalanche region. For the APDs, Vrt attains values of 60–70 V, over which not much more regions are depleted. Further increasing the bias voltage mainly leads to higher electric fields in the structure. As the highest electric field reaches the critical value Ecr, multiplication of carriers starts to occur. More rising in the reverse voltage makes the steady current density go up to in principle infinity, where actually the avalanche breakdown takes place [37, 38, 39]. The corresponding voltage here is thus named avalanche breakdown voltage Vbr, about 110 V for Si APDs.

Figure 4.

Dependence of dark current, photocurrent (with/without multiplication) and multiplication gain on reverse bias voltage of the designed SACM Si APD.

At operating voltage (Vbr × 98%), the multiplication region of a Si APD has an electric field as high as about 3.7 × 105 V/cm and an impact generation rate as high as about 2.8 × 1025 s−1 cm−3. As a result, the avalanche gain (M) and the sensibility (S) of the linear-mode Si APDs are observed to be up to about 600 and 300 A/W @905 nm respectively.

The ROIC chips were developed on the 0.18-μm CMOS platform of SMIC, Shanghai. The voltage-stabilized source effectively reduces external noise jamming and increases voltage suppression ratio of the power source. TIA shows trans-impedance of 120 dBΩ, the equivalent input noise is about 6 pA/Hz1/2, the rise time is 7.3 ns, and the bandwidth is BW ≥ 35 MHz.

Arrays of 64 × 1 Si APDs and ROIC chips were integrated to form the photodetector device by performing bonder-leading welding techniques. Together with packaging processing, the devices of 64 × 1 Si APD focal-plane arrays were successfully fabricated, one of which is shown in Figure 5. The power of input signal light is 0.9 nW (the duty cycle is 1/1000), and the maximum output voltage amplitude is 1.04 V. The devices present pulse responsivity R ≥ 1 × 106 V/W, noise equivalent power NEP ≤ 5 pW/Hz1/2, rise time tr ≤ 3 ns, and inhomogeneity of responsivity of each pixel ≤10%, under 905 nm, 100 ns and 10 kHz of laser irradiation.

Figure 5.

The device of 64 × 1 linear-mode Si APD focal-plane array.

2.4. Application of the linear mode Si APD array

As we constructed linear-mode Si APD focal-plane detectors, the 64 × 1 array devices are tested for possible applications. One example is that, the device is effective in running a driverless platform. Using this APD array, an obstacle-avoidance LADAR, as shown in Figure 6(a), is successful with detection distance of 110 m, distance resolution of 5 cm and angle resolution of 0.5o. This LADAR can effectively detect the obstacles on the way, as shown in Figure 6(b). Compared with traditional technique, in which a single detector was used, the image is much clearer (10 times of pixels) and the imaging speed is much faster (35 versus 15 Hz), so this newly developed obstacle-avoidance LADAR is more accurate and better to be used in driverless vehicles.

Figure 6.

(a) An obstacle-avoidance LADAR, using 64 × 1 linear-mode Si APD array as the focal-plane detector, is installed on a driverless vehicle. (b) The imaging effect of the LADAR, where the red pattern shows the existence of obstacles.


3. Si SPAD focal-plane arrays

A Geiger-mode APD can detect a signal as weak as a single photon. In recent years, it is very active and effective as a single-photon detector and usually termed SPAD. Organized into arrays, SPAD can be used in many systems such as LADAR, mobile laser imaging and viewing instrument. By using some special processing, we developed typical Si-based SPAD arrays working at 905 nm. The key techniques are described as follows.

3.1. Fabrication of Si SPAD array chip

3.1.1. Design

According to the requirement of a Si SPAD, there would be a depletion region as thick as 30 μm. Using usual single-sided abrupt p-n junction, to get such a large depletion while remaining avalanche gain, one need to apply a voltage as high as 500 V, which is not realistic enough. We design a reach-through structure, containing n+-π-p-π-p+ layers, as shown in Figure 7. The electric field distribution under bias near breakdown is similar to that in Figure 2(b). The light-generated carriers is multiplied in the region with highest electric field, so called multiplication region. This region is very thin compared to the whole depletion region. The other parts in the depletion region can have electric field as weak as possible but sufficient to ensure carrier drifting at the saturated speed. As a result, the operating voltage can be greatly decreased.

Figure 7.

The designed structure of Si SPAD.

There could be two types of host materials. One is high-resistive Si wafer (p-type with 1014 cm−3 boron doped), the other is epitaxial lowly-doped p-type Si on a p+-doped Si wafer. To fabricate the reach-through structure on bulk Si, the whole wafer has to be depleted as the device works. To get low breakdown voltage of 150 V or so, we must polish the wafer down to a thickness of 50 μm, which is difficult and brings a lot of unstable features to the devices. Therefore, we use epitaxial wafer as the material for fabricating Si SPAD arrays.

Decreasing the dark count rate (DCR) is conflicting with choosing lower avalanche electric field Emax and lower operating voltage. There should be tradeoff between these two to get optimized structure parameters. Our previous product has DCR of several 10 kHz with light receiving area of Φ500 μm. Here we try making a device area of Φ50 μm to get much lower DCR. By further well designing and optimizing the internal structure, improving the pixel uniformity and surface passivation effect, we succeeded in controlling the DCR under 10 kHz.

3.1.2. Precise control and uniformity

Precise control of the device structure is a decisive step. The multiplication layer and charge layer are most important because they greatly influence the key parameters such as quantum efficiency, response time and gain. To precisely control the charge quantity and multiplication length, we adopt the following process. Boron ion-implantation is firstly performed to accurately control the dopant amount. Due to smaller diffusion coefficient of boron in Si, the broadening of boron distribution after thermal treatment is weaker. Then, second epitaxy is carried out to grow n-type layer. This is the way to precisely control the multiplication length.

It is the critical process of a Si SPAD array to make all the pixels controllably consistent in characteristics, e.g. the avalanche gain, the response time, and the breakdown behavior. The pixels uniformity is influenced by four factors: epitaxial structure, ion-implantation of the charge layer, diffusion of the p-n junction, heat-induced doping redistribution in the device process. The most important is that, the epitaxial layers should be grown as uniform as possible. Usually, a 50-μm thick epitaxial layer should have thickness uncertainty of less than 50 nm.

In an SPAD array, one of the critical structures is the guard-ring. The designed SPAD structure shows that, the n+ contact and light-incident layer is so small that the p-n junction depth is about 0.5 μm. In order to prevent the device from being lowly broken down, a guard-ring around the device can be fabricated by doping at the edge of the n+-doped area. It uses more deeply diffused n-doping (Phosphorus has big diffusion coefficient in silicon) to decrease the curvature rate of the edge of the p-n junction with the π-region, to reduce the electric-field of the junction edge, to improve the breakdown voltage at the edge of the device. Properly controlling diffusion depth, doping level and ring width, the edge breakdown voltage can be improved to be about 1.5 times that in the avalanche region. As an example, usually, 2–3 μm wide guard-ring is suitable for the above purpose.

3.1.3. Isolation process

In the avalanche procedure, there come photons at wavelength shorter than 1 μm while creating multiplication carriers. These photons may enter other pixels nearby to generate unintentional count. In order to suppress cross-talk between pixels, design of the SPAD structure is optimized. Further by processing the light-hiding belt and with the aid of decreasing the reflectivity at the interface, the optical cross-talk is well controlled. In the Si SPAD array, there exist big shunt capacitances in between the adjacent pixels, so the electric cross-talk could be of high possibility. This is resolved by decreasing capacitances related to the wires in the later interconnection process.

With a long absorption region (~30 μm), it needs to keep background doping level being lower than 1014 cm−3 in order to remain lower electric field in this layer. In the meantime, the inversion layer near the surface is a major factor causing device failure. Thus, it is necessary to set a p+-doped area around the surficial active region, i.e. the channel-resisting region. It can cut the inversion layer at the SiO2-p-Si interface and stop the surficial expansion of the depletion region. Its final role is to suppress the surficial leakage current and to prevent the device from being broken down by a low bias. By calculation and experience, a doping level higher than 1016 cm−3 in the surface layer is sufficient to form the channel-resisting effect.

By performing the above processes, chips of 32 × 32–64 × 64 Si SPAD arrays are fabricated. One sample is shown in Figure 8.

Figure 8.

One chip of Si SPAD array.

3.2. ROIC optimization

To realize highly accurate timing of the photon arrival, we use a time-digit-conversion (TDC) circuit with the aid of phase-shifting technique. This approach satisfies the requirements of 2 ns in the time resolution and 20 kHz in frame frequency while decreasing the power consumption of the whole chip. An active-quenching design is used to reach an extinction time of less than 50 ns.

High precision timing and time-reading circuit is composed of TDC and memory readout unit, as shown in Figure 9. After gate/signal conversion circuit transforms the high voltage output of SPAD into low voltage signal, it firstly needs to read the photon flying time of every pixel through TDC, to change times into digits, and to read out the digits through the memory readout circuit. For the purpose of 2 ns of time resolution, highly frequent, highly stable main clock is supplied to the 12-digits counter. Frequency more than 500 MHz is not easy to be realized at every pixel at the same time. Due to the processing limit, shunt resistors and capacitors may significantly contribute to the power consumption (To an 32 × 32 array, the power consumed in clock lines would be 100 mW). Without using phase-locking loop (PLL), here we design TDC circuit with the aid of phase-shifting technique, to satisfy the requirement of high time resolution while decreasing the requirement for clock frequency. The TDC consists of counter and time-delay chain. When the external timing signal (a rising edge) comes, the counter starts to work; when a photon is detected, the circuit generates a signal to stop the counter and remain the present count data. Via the time-delay chain composed of 8 units, the external clock creates 8 clock signals with different phases. As the starting signal comes, every clock signal is sampled and the time-delay chain outputs an 8-digit signal, which will be coded and saved into the data process module. A similar process happens when a stop signal comes. Difference calculation between the start and stop data gives a 4-digit signal, constructing a 12-digit information together with the data from the counter. Then, the time interval between start and stop is converted into a 12-bit digital signal. When the counter output clock is effective, its locked state is transferred bit by bit into a 12-digit register. Controlled by a 25 MHz clock, the register transfers its digits into the register of the neighboring pixel. This serial transform mode gives at last a frame frequency of 25 kHz.

Figure 9.

Structure of the high precision time-digit conversion circuit designed for Si SPAD array.

Avalanche-quenching circuit is another important aspect in SPAD array. To realize a dead-time less than 50 ns, we design actively-quenching circuit as shown in Figure 10. When avalanche photocurrent is detected, the voltage at point a jumps down and forces the quenching circuit run. After a while, the voltage at point b comes higher, switching on the transistor M1 via the feedback branch, and quickly pulling down the voltage at point Vapd to make the APD bias lower than breakdown voltage (quenching the APD). After more a while, the charge–discharge circuit gradually decreases the voltage at point b, M1 is turned off via feedback circuit, and Vdd charges the APD through R0 and parasitic capacitances to restore the Geiger mode.

Figure 10.

Structure of the quenching circuit for Si SPAD array.

The above designs are realized by performing CMOS processing, and thus the ROIC chips are produced.

3.3. Interconnection and package

The next key processing is interconnecting the SPAD chip and the ROIC chip. The fabrication is as follows. After some degree of thinning processing, the backside of the wafer is treated to have a light-incident window for every pixel. As shown in Figure 11, the window area is formed by etching off the p-type substrate, and the layer under the etched window is made to be ~35 μm thick. Dry etching such as SF6 + O2 ICP is used to fabricate windows with straight sidewalls. In addition, the uniformity of this processing must remain ≤ ±2% in the window depth and ≤ ±0.5 μm in the window diameter. Then, it is realized that negligible light is absorbed by the p-type substrate. Using standard Indium-shot interconnection processing, the SPAD chip is connected exactly with the ROIC, as shown in Figure 11. Integrating the interconnected chip with TEC cooling cells, and packing these all into a vacuum can, a Si SPAD focal-plane device is developed.

Figure 11.

Schematic of the interconnection process for fabricating Si SPAD arrays.

3.4. Application of Si SPAD focal-plane arrays

Measurements show that the developed SPAD array have DCR lower than 5 kHz, average photon detection efficiency (PDE) ~25%, time resolution <1 ns, and frame speed ~25 kHz. The 32 × 32 SPAD array exhibits pixel uniformity < ±5% (e.g. in counting rate). It can thus be applied in a practical imaging system. The above fabricated SPAD array device is installed in a LADAR 3D imaging instrument. The instrument can three-dimensionally detect and recognize multi-hided objects in forests and mountains, and be of small size, light weight, high resolution and rapid imaging.


4. InGaAsP/InP SPAD focal-plane arrays

Like InGaAs/InP SPADs [40], InGaAsP/InP SPADs are also extensively studied and practically explored for their utility in many fields including single photon imaging [41] and quantum information processing [42, 43] in the near-infrared wavelength range. Thanks to the advanced epitaxial techniques, this kind of SPADs has been well developed and applied in, e.g., LADAR in recent years [41, 44]. Nevertheless, many critical problems are still open to be resolved. One of them is the device homogeneity, such as the reproducibility and uniformity of the SPAD performance [45], which are strongly required to be precisely controlled by refining the structure parameters in epitaxial growing and device processing. One can, of course, estimate the effects of some individual physical parameters on the performance homogeneity using some analytical method [46], but it is not easy to obtain the knowledge of many parameters at the same time. It is even unlikely to make clear the collective influence of multiple parameters and to quantitatively take a balance between various parameters. The quantitative association between the device inhomogeneity and structure uncertainty should thus be necessarily established. Therefore, we firstly carry out a statistical analysis on InGaAsP/InP SPAD characteristics by randomizing the structure parameters, and then figure out the necessary control accuracies in a few significant structure parameters, which are required for nice device homogeneity. Accordingly, we fabricate InGaAsP/InP SPAD focal-plane arrays and demonstrate their applications.

4.1. SPAD array homogeneity versus material controllability

An InGaAsP/InP APD structure is designed as an example SPAD object. It is of a hetero-structure comprising SACM layers, as can be seen schematically in Figure 12(a). By using conventional APD theory [28, 47] and lately advanced approaches [48, 49, 50], citing material parameters from previous reports [49, 51, 52], and neglecting the dead-space effect [53], the device characteristics of this structure was calculated.

Figure 12.

(a) The structure diagram of a designed InGaAsP/InP Geiger-mode APD (or SPAD) studied in this work; (b) calculated I-V relationship in dark and (c) calculated DCR versus PDE for the InGaAsP/InP SPAD device at 230 K.

Figure 12(b) illustrates the calculated current–voltage (I-V) characteristics in dark. Here the breakdown voltage Vb is principally the applied reverse bias V at infinite avalanche current. The simulated DCR rd versus PDE η is shown in Figure 12(c). Both of them are dependent on the applied excess bias Vex = VVb. As the SPAD is running under a middle excess bias Vex0 = 5 V, rd is found below 10 kHz and η appears some 0.50. It suggests that 5 V of Vex is an optimal operating condition at 230 K, so Vex = 5 V will be the reference point in this study.

As a SPAD array is used, there is usually a common bias V0 generally applied to more than thousands of pixels [44]. Provided structure fluctuations exist among the numerous pixels, the effective Vex will vary from this to that pixel so that device performance exhibits inhomogeneity. To clarify this effect, we first set any structure parameter t randomly changing in a way as

t i = t 0 1 + W σ i , E1

where the subscription i = 1, 2, 3 is the calculation sequence number in a series of simulations, ti is the ith value of parameter t, t0 is the designed value of t, W is the fluctuation degree of t relative to t0, and σi is the ith value of the random variable σ, distributed in a normal mode with full width at half maximum (FWHM) of unity. Similarly we can set

n c i = n c 0 1 + W nc σ nc i , t m i = t m 0 1 + W tm σ tm i , E2
and so forth, where nci(tmi) is the ith value of charge density nc (multiplication width tm), nc0(tm0) is the designed value of nc(tm), Wnc(Wtm) is the relative FWHM of nc(tm) with respect to nc0(tm0), and σnci(σtmi) is the ith value of σ for changing nc(tm). All the variables are defined in a similar way to the above. The structure parameters are changing independently because each has its own FWHM and σ values. With a set of structure parameters (nci, tmi,…), one set of device performance data is calculated. With thousand sets of device performance data, distributions of Vb, Vex, rd, and η are obtained through statistics, and then the correlation between the performance fluctuations and device structures is figured out.

Our simulations show that nc, tm and tc (charge layer thickness) have strong effects, while absorption layer doping level na, thickness ta and multiplication layer doping level nm have weak effects on Vex. The strong tm effect can be easily understood since the width of the multiplication region is crucial to determine the characteristics of a SPAD [48, 54]. It means that the charge quantity should be controlled most precisely in design and epitaxy process. In addition, we see that every single structure parameter leads to Vex fluctuation in a linear manner.

Extending the above simulation to more structure parameters, the homogeneity of device performance can be obtained in terms of varying parameter numbers. The result of Vex is displayed in Figure 13(a), where a sublinear change of Vex fluctuation is seen to happen with increasing parameter number. Taking all of the six parameters into account, we get that Vex varies with a FWHM of 24%, far less than a simple summation of the effects of individual parameters. The Vex distribution arisen by six independently varying parameters is demonstrated in Figure 13(b). With the referred excess bias Vex0 = 5 V, the practical value of Vex varies mainly in the range of 4.4–5.6 V, quite good for many applications. It may be also worthy to get the effects on other performance characters. Figure 13(c) presents the variation of DCR rd, which exhibits a roughly normal distribution with a wide relative FWHM (54%). In detail, the worst DCR is some 30% higher than the designed value, which is acceptable in applications. Figure 13(d) displays the distribution of PDE η, normal with a narrow relative width (17%). It suggests a PDE change within 0.46–0.54, which is homogeneous enough in many applications. The reason why the DCR fluctuation is much larger than PDE is that DCR depends almost exponentially on Vex and PDE, as Figure 12(c) shows.

Figure 13.

(a) Simulated distribution FWHM of the excess bias Vex dependent of the number of simultaneously varying parameters; distribution histograms of (b) Vex, (c) DCR and (d) PDE brought about by six simultaneously varying parameters including the doping level and the width of the absorption region na and ta, of the charge layer nc and tc, and of the multiplication region nm and tm, respectively. The dashed lines on the histograms indicate the fitted results to normal distributions.

In realistic epitaxial growth, controlling the thickness or the doping level may usually have a common precision in different layers, although a non-dope layer would have worse carrier density fluctuation than the intentionally doped regions. The following thus shows a way more effective to investigate the uncertainty correlation between the epitaxy growth and the device performance. Let us study the device characteristics varying with common fluctuation width Wnd of the residual carrier densities in the absorption layer and the multiplication layer, common fluctuation width Wt of the widths of the absorption, charge and multiplication regions, and fluctuation width Wn of the doping level in the charge layer, that reads

n a i = n a 0 1 + W nd σ na i , n m i = n m 0 1 + W nd σ nm i , n c i = n c 0 1 + W n σ nc i , t a i = t a 0 1 + W t σ ta i , t m i = t m 0 1 + W t σ tm i , t c i = t c 0 1 + W t σ tc i . E3

One typical trial is to examine the dependence on two fluctuating parameters while remaining the others fixed. With Wnd fixed, Figure 14 shows the Vex contours as functions of Wn and Wt. From these data, the precision range of epitaxy growth required for definite performance homogeneity can be clearly seen. At first, two conditions with 10 and 20% of Wnd appear close to each other, especially for high Vex variations, owing to the weak effect of the carrier density in non-dope layer. The relationship between Wn and Wt is far away from a linear curve but more like a circle for a finite Vex fluctuation. To constraint Vex fluctuation below a certain value, the fluctuations in thickness and charge control should roughly follow

W n 2 + W t 2 < W x 2 , E4

where Wx is a certain precision control value of thickness and charge. Quantitatively speaking, Vex relative fluctuations below 50, 40, 30 and 20% need Wx values of about 4, 3, 2 and 1%, respectively.

Figure 14.

Contours of Vex fluctuation width as a function of the thickness and doping level fluctuation widths Wt and Wn, under different fixed fluctuation widths of the residual carrier density in the non-dope layers Wnd.

In conventional growth, it is more difficult to control doping than thickness. Based on the result of Figure 14, as the thickness control can be better than 1–2%, Vex homogeneity of 50, 40 and 30% could be realized by constraining the charge precision within 4–4.5, 3–3.6 and 2–2.7%, respectively. Viewed from another angle, the result is suggestive of a large space to tradeoff between the controls in thickness and charge. The example of Vex fluctuating below 50% with Wnd = 10% suggests that the thickness (charge) precision Wt(Wn) is better to be as small as 1% if Wn(Wt) just satisfies 4.5%(4%), while Wt(Wn) could be roughened to 3%(3.5%) if Wn(Wt) weakly decreases to be about 3.5%(3%). In general, limiting the device inhomogeneity (in term of Vex) below 50, 40, 30, and 20% needs the thickness and charge be controlled to a precision degree better than 3–3.2, 2.4–2.6, 1.7–1.9, and 0.9–1.2%, respectively. Since these degrees of control accuracies are easy or possible in epitaxy growth, InGaAsP/InP SPAD arrays are now producible in many laboratories [55, 56, 57] including our group, as will be described below. In order to finely limiting the device homogeneity, such as with Vex fluctuation less than 10%, the thickness and charge should be controlled better than 0.5% in fluctuation, together with non-dope carrier density controlling within 10%. This degree of epitaxy precision is quite a challenging technique. It is possibly one of the reasons why it is presently still difficult to prepare 512 × 512 or larger scale SPAD arrays. Obviously, the above method is very helpful and effective to quantitatively correlate the controllability of multiple structure parameters with the SPAD device homogeneity.

4.2. Fabrication of InGaAsP/InP SPAD arrays

InP based APDs must use epitaxial materials. Firstly, we prepared APD materials with the main structure as shown in Figure 12(a) by using metal organic chemical vapor deposition (MOCVD). MOCVD growth is performed to satisfy the material uniformity requirements described above. On this epitaxial wafer, SPAD device structure as shown by Figure 15 will be fabricated. As an array, there are isolating grooves (channels) between the pixels, Indium shots on the front side for interconnection and micro-lenses on the backside for light collection.

Figure 15.

Structure of InGaAsP/InP SPAD array chip.

The key processes to fabricate the InGaAsP/InP SPAD arrays are as follows. First the active p-n junction is formed by selective diffusion. The diffusion process includes, thermally evaporating one layer of solid Zn3P2, depositing one layer of SiNx to thoroughly cover Zn3P2, rapid thermal annealing to diffuse Zn into the chip, and etching off the SiNx and the resident Zn3P2. To make the response time of a SPAD device shorter than 5 ns, the active area (the p-n junction area) of a pixel is made to be less than ϕ100 μm. Considering the guard-rings, the lateral depletion width and the diffusion length of electrons and holes, the distance between neighboring pixel centers is taken to be 300 μm. To suppress the cross-talk between pixels, the isolation between pixels is, besides the deep grooves, aided by highly resistive p-n junction. To increase the filling factor, light is incident on the backside, where there fabricated micro-lenses for each pixel. The microlenses here are not bonded onto the backside, but directly fabricated on the backside by specific dry-etching.

The fabricated SPAD array is characterized as shown in Figure 16. Measurements on material properties show that residual carrier density, layer thickness, and doping level fluctuations in a 10 × 10 mm2 area appear about 8, 0.8 and 1.5%. On such a chip, 32 × 32–64 × 64 arrays of SPADs were developed and characterized at low temperatures. Under gated mode with gate repetition rate of 500 kHz and gate width of 10 ns, DCR and PDE were measured using a single-photon laser at 1.06 μm. The afterpulsing probability is controlled below 2% by setting the dead time to be ~2 μs. As presented by the inset in Figure 16, various pixels have dark I-V curves with Vb (defined to be the bias at 10 μA) weakly changing but around 75 V. Figure 16 shows that, the fluctuation in the excess bias distributes in a normal way with FWHM of 14%, which is consistent with the simulated value 18%, and compatible with an estimation based on Figure 16. The DCR and PDE vary normally with FWHM of 31 and 10%, and this is also consistent with the simulated values 37 and 12%, respectively.

Figure 16.

An experimental result of Vex fluctuation distribution of InGaAsP/InP SPAD array. The dashed line represents the fitting to a normal distribution. The inset exhibits the I-V curves in dark of a few typical InGaAsP/InP SPAD devices.

The ROIC is designed in a way similar to that of Si SPAD arrays. The interconnection is a standard indium-shot inversion-bonding process. After packaging into a vacuum can with a transparent window, the InGaAsP/InP SPAD array is developed and can be used in an imaging system.

4.3. Application of the InGaAsP/InP SPAD arrays

A 64 × 64 InGaAsP/InP SPAD array device is installed onto the focal plane of a LADAR system. Under 1.06 μm laser irradiation, the scence 1–3 km away was successfully imaged with 3D information, as shown in Figure 17.

Figure 17.

Imaging by a LADAR instrument with 64 × 64 InGaAsP/InP SPAD focal-plane array detectors.


5. Summary

APDs are being more and more widely and deeply studied to satisfy the requirement in weak light and single photon imaging. The progresses of this worldwide study, especially the distinctive researches and achievements in SITP and UESTC are reviewed. We successfully fabricated up to 64 × 1 linear-mode Si APD array, and 32 × 32–64 × 64 Si SPAD arrays, and applied them in LADAR platforms like driverless vehicles. Also, we developed 32 × 32-64 × 64 InGaAsP/InP SPAD arrays, and constructed 3D imaging LADAR using them. Together with the progresses of other groups and other materials, we see a prospective future for the development and application of focal-plane APDs.



This work was partially supported by 1000 Foreign Experts Program, China, and by the 1000 Talents Plan of Sichuan Province. We would like to thank Profs. L.B. Yu, Z. Shi, O. Wang, X. Li, Q. Dai, Y. Yang and Drs. J. Chen, X. Xie and Q. Xu at SITP for their practical help, and Profs. Zh. M. Wang, Q. Zhou and G.W. Deng at UESTC for their theoretical assistance.


  1. 1. Campbell JC. Recent advances in avalanche photodiodes. Journal of Lightwave Technology. 2016;43(2):278-285
  2. 2. Woodson ME, Ren M, Maddox SJ, Chen YJ, Bank SR, Campbell JC. Low-noise AlInAsSb avalanche photodiode. Applied Physics Letters. 2016;108:081102
  3. 3. Campbell JC. Advances in photodetectors. In: Kaminow I, Li T, Wilner AE, editors. Optical Fiber Telecom-munications, Part A: Components and Subsystems. 5th ed. Vol. 5. San Francisco, CA, USA: Academic; 2008
  4. 4. Bertone N, Clark WR. APD arrays—avalanche photodiode arrays provide versatility in ultra-sensitive applications. Laser Focus World. 2007;43(9):69-73
  5. 5. Mitra P, Beck JD, Skokan MR, Robinson JE, Antoszewski J, Winchester KJ, et al. Adaptive focal plane array (AFPA) technologies for integrated infrared microsystems. In: Proceedings of SPIE 6232. SPIE - The International Society for Optics and Photonics. 2006. pp. 62320G-1-62320G-11
  6. 6. Tosi A, Calandri N, Sanzaro M, Acerbi F. Low-noise, low-jitter, high detection efficiency InGaAs/InP single-photon avalanche diode. IEEE Journal of Selected Topics in Quantum Electronics. 2014;20(6):3803406
  7. 7. Jiang XD, Itzler M, O’Donnell K, Entwistle M, Owens M, Slomkowski K, et al. InP-based single-photon detectors and geiger-mode APD arrays for quantum communications applications. IEEE Journal of Selected Topics in Quantum Electronics. 2015;21(3):3800112
  8. 8. Aull B. Geiger-mode avalanche photodiode arrays integrated to all-digital CMOS circuits. Sensors. 2016;16:495. DOI: 10.3390/s16040495
  9. 9. Othman MA, Taib SN, Husain MN, Napiah ZAFM. Reviews on avalanche photodiode for optical communication technology. ARPN Journal of Engineering and Applied Sciences. 2014;9(1):35-44
  10. 10. Diagne M, McIntosh A, Donnelly J, Myszka M, Frechette J. Advances in InP/InGaAs Geiger-mode APD focal plane arrays. In: Proceedings of SPIE 10659. SPIE - The International Society for Optics and Photonics. 2018. p. 1065904
  11. 11. Itzler M, Salzano G, Entwistle M, Jiang XD, Owens M, Piccione B, et al. Asynchronous Geiger-mode APD cameras with free-running InGaAsP pixels. In: Proceedings of SPIE 10212. SPIE - The International Society for Optics and Photonics. 2017. p. 102120K
  12. 12. Cai X, Wu C, Lu H, Chen Y, Zhou D, Liu F, et al. Single photon counting spatial uniformity of 4H-SiC APD characterized by SNOM-based mapping system. IEEE Photonics Technology Letters. 2017;29(19):1603-1606
  13. 13. Steindl B, Hofbauer M, Schneider-Hornstein K, Brandl P, Zimmermann H. Single-photon avalanche photodiode based fiber optic receiver for up to 200 Mb/s. IEEE Journal of Selected Topics in Quantum Electronics. 2018;38013:3801308. DOI: 10.1109/JSTQE.2017.2764682
  14. 14. Ceccarelli F, Acconcia G, Labanca I, Gulinatti A, Ghioni M, Rech I. 152-dB dynamic range with a large-area custom-technology single-photon avalanche diode. IEEE Photonics Technology Letters. 2018;30(4):391-394. DOI: 10.1109/LPT.2018.2792781
  15. 15. Chen J, Yang Y, Guo CD, Kong FL, Kou X, Wang O, et al. Linear-mode linear arrays 16 pixel silicon avalanche photodiodes with high gain and low noise readout. In: Proceedings of SPIE 10462. SPIE - The International Society for Optics and Photonics. 2017. p. 104625E
  16. 16. Song HZ, Deng J, Dai Q, Shi Z, Yu LB. Quantitative correlation between fabrication precision and device homogeneity of single-photon avalanche diodes. IEEE Transactions on Electron Devices. 2016;63(12):4845-4851. DOI: 10.1109/TED.2016.2618222
  17. 17. Dai M, Li X, Shi Z, Dai Q, Song HZ, Tang Z, et al. Gain and noise properties of multi-gain-stage superlattice InGaAs avalanche photodiode. Infrared and Laser Engineering. 2016;45(5):0520009
  18. 18. Dai Q, Deng J, Shi Z, Yu LB, Song HZ. Multi-parameter fluctuation effects on InGaAsP/InP Geiger-mode avalanche photodiodes. Advances in Computer Science Research. 2015;10:965-969
  19. 19. Chen J, Dai Q, Shen B, Deng J, Xu Q, Shi Z, et al. Multi-parameter study on the performance fluctuation of InGaAsP/InP single photon avalanche diodes. In: Proceedings of 2016 IEEE International Conference on Integrated Circuits and Microsystems (ICICM 2016). IEEE - The Institute of Electrical and Electronics Engineers. 2017. pp. 186-190
  20. 20. Hao X, Chen J, Kong F, Ke Z, Wang O, Yu LB, et al. Design and photoelectric properties of silicon Geiger-mode APD array. In: ALAT 2017, Shenzhen, China: Laser Technology Press; May 20-21, 2017
  21. 21. Williams GM, Huntington AS. Probabilistic analysis of linear mode vs Geiger mode APD FPAs for advanced LADAR enabled interceptors. In: Proceedings of SPIE 6220. SPIE - The International Society for Optics and Photonics. 2006. p. 622008
  22. 22. Kishimoto S, Yonemura H, Adachi S, Shimazaki S, Ikeno M, Saito M, et al. A fast X-ray detector using silicon avalanche photodiodes of 64-pixel linear array. Journal of Physics: Conference Series. 2013;425:062007
  23. 23. Sultana A, Kamrani E, Sawan M. CMOS silicon avalanche photodiodes for NIR light detection: a survey. Analog Integrated Circuits and Signal Processing. 2012;70:1-13
  24. 24. Niclass C, Rochas A, Besse P, Charbon E. Design and characterization of a CMOS 3-D image sensor based on single photon avalanche diodes. IEEE Journal of Solid-State Circuits. 2005;40(9):1847-1854
  25. 25. Brian FA. Silicon Geiger-mode avalanche photodiode arrays for photon-starved imaging. In: Proceedings of SPIE 9492. SPIE - The International Society for Optics and Photonics. 2015. p. 94920M
  26. 26. Silvaco International Software. Atlas User’s Manual [Z]. Santa Clara, CA, USA: Silvaco; 2012
  27. 27. Tsang WT. Semiconductors and Semimetals. New York: Academic Press; 1985. pp. 247-328
  28. 28. McIntyre RJ. On the avalanche initiation probability of avalanche diodes above the breakdown voltage. IEEE Transactions on Electron Devices. 1973;ED-20(7):637-641
  29. 29. McIntyre RJ. Multiplication noise in uniform avalanche diodes. IEEE Transactions on Electron Devices. 1966;ED-20(7):637-641
  30. 30. Thomas O, Yuan ZL, Shields AJ. Practical photon number detection with electric field-modulated silicon avalanche photodiodes. Nature Communications. 2012;3(644):1-5
  31. 31. Thomas O, Yuan ZL, Dynes JF, Sharpe AW, Shields AJ. Efficient photon number detection with silicon avalanche photodiodes. Applied Physics Letters. 2010;97:031102
  32. 32. Wang W, Feng Q, Wu W, Xie YT, Wang Z, Feng SJ. Analysis and simulation of process and performance of silicon avalanche photodiode. Infrared and Laser Engineering. 2014;43(1):140-144
  33. 33. Wang W, Yan LS, Wang C, Du CY, Wang T, Wang GY, et al. Analysis of separate absorption charge multiplication Ge/Si-APD. Infrared and Laser Engineering. 2015;44(4):1349-1353
  34. 34. Wegrzecka I, Wegrzecki M, Grynglas M, Bar J, Uszynski A, Grodecki R. Design and properties of silicon avalanche photodiodes. Opto-Electronics Review. 2004;12(1):95-104
  35. 35. Mukherjee PP. Design, performance analysis and application of a modified improved MOS regulated cascode circuit. Procedia Technology. 2012;4:477-481
  36. 36. Ben-Esmael M, Mathew M, Hart BL, Hayatleh K. Technique for increasing the output impedance of CMOS regulated cascode circuits. Journal of Circuits, Systems and Computers. 2016;25(10):1650118
  37. 37. Wegrzecka I, Grynglas M, Wegrzecki M. Spectral dependence of the main parameters of silicon avalanche photodiodes. In: Proceedings of SPIE 4516. SPIE - The International Society for Optics and Photonics. 2001. pp. 187-193
  38. 38. Kunar R, Chamberlain SG, Roulston DJ. Two-dimensional computer simulation of the breakdown characteristics of a multi-element avalanche photodiode array. IEEE Transactions on Electron Devices. 1984;ED-31(7):928-933
  39. 39. Webb PP, Mcintyre RJ. Multi-element reach through avalanche photodiodes. IEEE Transactions on Electron Devices. 1984;ED-31(9):1206-1212
  40. 40. Boso G, Korzh B, Lunghi T, Zbinden H. Low noise InGaAs/InP single-photon negative feedback avalanche diodes: Characterization and applications. In: Proceedings of SPIE 9492. SPIE - The International Society for Optics and Photonics. 2015. pp. 94920Q-1-94920Q-11
  41. 41. Kutteruf MR, Lebow P. 1541 nm Gm APD LADAR system. In: Proceedings of SPIE 9080. SPIE - The International Society for Optics and Photonics. 2014. pp. 908006-1-908006-7
  42. 42. Tosi A, Zappa F, Cova S. Single-photon detectors for practical quantum cryptography. In: Proceedings of SPIE 8542. SPIE - The International Society for Optics and Photonics. 2012. pp. 85421U-1-85421U-1-8
  43. 43. Zhang J, Itzler MA, Zbinden H, Pan JW. Advances in InGaAs/InP single-photon detector systems for quantum communication. Light: Science & Applications. 2015;4:e286-1-e286-13
  44. 44. Itzler MA, Krishnamachari U, Entwistle M, Jiang XD, Owens M, Slomkowski K. Dark count statistics in Geiger-mode avalanche photodiode cameras for 3-D imaging LADAR. IEEE Journal of Selected Topics in Quantum Electronics. 2014;20:3802111-1-3802111-11
  45. 45. Itzler MA, Krishnamachari U, Chau Q, Jiang XD, Entwistle M, Owens M, et al. Statistical analysis of dark count rate in Geiger-mode APD FPAs. In: Proceedings of SPIE 9250. SPIE - The International Society for Optics and Photonics. 2014. pp. 925003-1-925003-12
  46. 46. Savuskan V, Brouk I, Javitt M, Nemirovsky Y. An estimation of single photon avalanche diode (SPAD) photon detection efficiency (PDE) nonuniformity. IEEE Sensors Journal. 2013;13:1637-1640
  47. 47. McIntyre RJ. Multiplication noise in uniform avalanche diodes. IEEE Transactions on Electron Devices. 1966;13:164-168
  48. 48. Ramirez DA, Hayat MM, Itzler MA. Dependence of the performance of single photon avalanche diodes on the multiplication region width. IEEE Journal of Quantum Electronics. 2008;44(12):1188-1195
  49. 49. Donnelly JP, Duerr EK, McIntosh KA, Dauler EA, Oakley DC, Groves SH, et al. Design considerations for 1.06-μm InGaAsP–InP Geiger-mode avalanche photodiodes. IEEE Journal of Quantum Electronics. 2006;42(8):797-809
  50. 50. Javitt M, Savuskan V, Merhav T, Nemirovsky Y. Revisiting single photon avalanche diode current-voltage modeling and transient characteristics. Journal of Applied Physics. 2014;115:204503-1-204503-9
  51. 51. Tan LJJ, Ng JS, Tan CH, David JPR. Avalanche noise of electrons and holes in (100)-oriented Ga1−xInxAsyP1−y. IEEE Journal of Quantum Electronics. 1985;21(9):1326-1337
  52. 52. Itzler MA, Jiang XD, Ben-Michael R, Slomkowski K. InGaAsP avalanche photodetectors for non-gated 1.06 μm photon-counting receivers. In: Proceedings of SPIE 6572. SPIE - The International Society for Optics and Photonics. 2007. pp. 65720G-1-65720G-10
  53. 53. Hayat MM, Saleh BAE, Teich MC. Effect of dead space on gain and noise of double-carrier-multiplication avalanche photodiodes. IEEE Transactions on Electron Devices. 1992;39(3):546-552
  54. 54. Lee K, Yang K. Analysis of InGaAs/InP single-photon avalanche diodes with the multiplication width variation. IEEE Photonics Technology Letters. 2014;26(10):999-1002
  55. 55. Bronzi D, Villa F, Tisa S, Tosi A, Zappa F, Durini D, et al. 100000 frames/s 64 × 32 single-photon detector array for 2-D imaging and 3-D ranging. IEEE Journal of Selected Topics in Quantum Electronics. 2014;20(6):3804310-1-3804310-10
  56. 56. Smith GM, McIntosh KA, Donnelly JP, Funk JE, Mahoney LJ, Verghese S. Reliable InP-based Geiger-mode avalanche photodiode arrays. In: Proceedings of SPIE 7320. SPIE - The International Society for Optics and Photonics. Vol. 7320. 2009. pp. 73200R-1-73200R-10
  57. 57. Itzler MA, Entwistle M, Owens M, Patel K, Jiang XD, Slomkowski K, et al. Comparison of 32 × 128 and 32 × 32 Geiger-mode APD FPAs for single photon 3D LADAR imaging. In: Proceedings of SPIE 8033. SPIE - The International Society for Optics and Photonics. 2011. pp. 80330G-1-80330G-10

Written By

Hai-Zhi Song

Submitted: April 10th, 2018 Reviewed: September 4th, 2018 Published: November 5th, 2018