|
1.IntroductionThere is a growing need in science and industry for large-scale detector arrays capable of imaging with high sensitivity and high speed near room temperature, over a wide range of natural illumination conditions from direct sunlight to a cloudy moonless night. Existing solid-state camera technology does not enable imaging in the visible spectrum over such a wide dynamic range near room temperature. In this design study, we demonstrate through detailed calculation means that a novel back-illuminated, wide dynamic range avalanche photodiode (APD) device design will enable imaging with high sensitivity near room temperature at , over the full range of natural illumination conditions from the AM 0 in space orbit around Earth, to ultralow light levels when the illuminance at the camera lens , as on a cloudy moonless night when far from city lights. The high pixel sensitivity is achieved using a compact readout circuit, capable of controlling detector operation in linear and Geiger-mode regimes. The proposed wide dynamic range APD detector can operate in full depletion mode in the linear regime, when biased below the reverse breakdown voltage to support high quantum efficiency and speed. When charged to reverse bias, the detector design was shown in Ref. 1 to support Geiger-mode, and the efficiency of collecting photogenerated electrons to produce detection events was evaluated by Monte Carlo simulation to be above 80% for wavelengths between 300 and .1 The silicon APD array was designed with square pixels of side length for high resolution imaging with 100% sensitive-area-fill-factor, on a mechanically rugged and optically transparent sapphire substrate. When (100) silicon is properly etched with potassium hydroxide:isopropyl alcohol:water or tetramethyl ammonium hydroxide (TMAH) solution through a thermally grown oxide mask, square-based pyramidal frustum or mesa arrays result with the four mesa sidewalls of the APD formed by (111) silicon planes that intersect the (100) planes at a crystallographic angle . The APD would be fabricated in the mesa using conventional silicon processing technology. The proposed APD design enables operation over a wide dynamic range in linear mode with or without internal gain for passive imaging in daylight conditions, and in single-photon sensitive Geiger-mode for active or passive imaging under low-light level conditions.1, 2 The detector is back-illuminated through a light focusing spherical microlens fabricated in the thinned sapphire substrate that directs light toward the device interior and away from the silicon mesa sidewalls, thereby compensating the optical dead space between pixels resulting from the mesa isolation etch and also from the internal device guard diode geometry. The sapphire microlenses can be coated with an antireflective layer, for example, to improve the optical power transmittance from air into the sapphire. In addition, a lattice-matched (AlN) antireflective crystal layer can be grown at the sapphire-silicon interface to further increase the optical power transmittance into the device silicon, as shown in Fig. 1 and described in Ref. 2. The APDs in the array share a common, front-side anode contact, made locally at the base of each device mesa. A low resistance aluminum (Al) or copper (Cu) metal anode grid fills the space between pixels and also functions to inhibit optical crosstalk by shielding adjacent detectors from stray photons generated during the avalanche process.3 SOS-APD arrays are indium bump-bonded to complementary metal oxide semiconductor (CMOS) readout ICs to produce hybrid focal plane arrays (FPAs), as shown in Fig. 2 , which might be or larger.4 The camera lens focal length and focal ratio maintain 100% sensitive-area-fill-factor over the full FPA active area by limiting the maximum incidence angle of the extreme rays shown in red in Fig. 2, incident at the microlens edges of the pixel, located farthest from the camera optic axis at the corner of the FPA.2 The silicon APD is designed with a carefully optimized high electric field multiplication region, fabricated by high temperature diffusion of boron, described in Ref. 1. The boron impurity profile determines the spatial electric field in the device, and the boundary between the low electric field absorption-collection region and the high electric field multiplication region, as well as the APD guard diode geometry described by the guard diode boundary function , are shown in Fig. 1. The narrow, high electric field multiplication region and the much wider low electric field collection-absorption region, help to provide stable, uniform gain with minimum variance, for nearly all of the photogenerated electrons in the device. Since most electron-hole pairs (EHPs) are generated through light absorption outside of the multiplication region, this results in a low excess noise factor for the APD. 5, 6, 7 2.Device Characteristics of Wide Dynamic Range Avalanche PhotodiodeThe silicon-on-sapphire material system was originally developed for integrated circuit applications; however, it also has many ideal attributes for use as a substrate, enabling back-illuminated, solid-state, silicon detector arrays. Sapphire is an anisotropic, dielectric crystal of the negative uniaxial type that is weakly birefringent and possesses broadband optical transmittance ranging from the deep ultraviolet to the midwave IR . Sapphire is extremely resilient, supporting thinning below , which is an important requirement for high resolution back-illuminated detector arrays. Sapphire can be optically polished to better than an 80 to 50 (scratch-dig) surface finish and can be etched optically using excimer laser ultraviolet light to fabricate light focusing microlenses for the silicon detectors. Sapphire is chemically resistant to most liquid etchants at room temperature, and therefore functions as an ideal etchstop material during liquid crystallographic etching with or TMAH solution to define the silicon pixel mesa arrays. To enable high quantum efficiency back-illuminated silicon detector arrays, the refractive index mismatch between air, sapphire, and silicon has to be corrected. The wide bandgap semiconductor material aluminum nitride (AlN) is closely lattice matched and refractive index matched to both sapphire and silicon, and offers the prospect of enabling fabrication of high transmittance (100) silicon-on-(AlN)-sapphire substrates for back-illuminated silicon imagers. 8, 9, 10 The APD in Fig. 1 has been designed for wide dynamic range operation by allowing the polysilicon cathode to extend beyond the edges of the mesa top over the thermally grown passivation oxide to form an MOS capacitor that acts as a charge-coupled device (CCD) gate. When the APD is reverse biased, the combination of the APD -junction and the MOS capacitor formed by the overhanging polysilicon operates similarly to the input diode and input gate of a CCD.11 When the wide dynamic range APD in Fig. 1 is reverse biased at voltages well below the reverse breakdown voltage of the detector, the -junction cannot deplete fully to the silicon-sapphire interface due to blocking by the high electric field boron impurity diffusion. However, the polysilicon MOS capacitor can deplete into the -type silicon fully to collect any photogenerated electrons in the device. The APD -junction prevents an inversion layer from forming under the polysilicon CCD gate by collecting the photoelectrons under the gate, producing an unmultiplied photocurrent. Therefore, at low reverse bias voltages, the polysilicon MOS capacitor depletes deep into the silicon to collect photogenerated electrons, which are in turn collected by the cathode of the APD -junction, bypassing the high electric field region of the APD without being multiplied. At low reverse bias, the photogenerated electrons are mainly directed toward the MOS capacitor at the periphery of the device and relatively few photogenerated electrons in the APD reach the high electric field multiplication region. Therefore, the internal gain of the device remains low, even though the electric field in the multiplication region might reach the critical value of in silicon required for impact ionization. The region at the periphery of the device is indicated as a guard diode, because its purpose is to collect surface generated electrons. The electrons collected at the periphery will have reduced probability of drifting to the high field region to be multiplied. When the APD is biased above the reverse breakdown voltage, a much larger fraction of photogenerated electrons are directed toward the multiplication region of the device than at lower bias voltages, resulting in larger internal gain. For both linear and Geiger-mode operation, the device silicon will be fully depleted to enable efficient collection of photogenerated electrons from the entire silicon volume. These results were confirmed through extensive device simulations using MEDICI software of mesa APD structures created using the TSUPREM 4 fabrication simulation software as described in Refs. 1, 2. 2.1.Calculated Photocurrent of Wide Dynamic Range Mesa Avalanche PhotodiodeThe sensitivity of the wide dynamic range APD in Fig. 1 can be evaluated by considering that the silicon detector images primarily over the visible wavelength range. Manufacturers of commercial consumer cameras that detect primarily in the visible spectral range often describe the camera sensitivity by specifying the minimum illuminance condition at the entrance aperture of the camera in photometric units of lux required to capture an image with clarity. Table 1 lists illuminance values for different imaging conditions.12 Table 1Natural illumination conditions and the corresponding illuminance EI in photometric units at the entrance aperture or camera lens of the imaging system.
The natural illuminances in Table 1 represent values at the entrance aperture of an imaging system of the type shown in Fig. 2. The range of illuminance values for direct sunlight listed in Table 1, existing at the camera lens aperture, are produced by solar spectral irradiance at airmass AM 0 and other airmass conditions arising for different angles of the sun above the Earth’s horizon. These spectral irradiance curves can be used to calculate the total photon flux incident on the wide dynamic range mesa APD pixel that produces a photocurrent in response to the illumination. Under AM 0 illumination, when the wide dynamic range APD camera stares directly at the sun from Earth’s orbit in space, the AM 0 spectral irradiance is given by the plot in Fig. 3 .13 Figure 4 shows the geometry of imaging the sun directly using the APD-FPA camera with a lens focal length and focal ratio setting , required for 100% sensitive-area-fill-factor over the full FPA area, as calculated in Ref. 2. The dashed blue curve in Fig. 3 represents the calculated spectral irradiance of the solar black body obtained from Eq. 1, which is derived from the Planck black body law. In Eq. 1, the average distance of the Earth from the sun is given by and the radius of the sun . The G173-03 standard described in Fig. 3 is produced by the American Society for Testing and Materials (ASTM) for photovoltaic performance evaluation and also defines in addition to AM 0 (for orbit around the Earth), two standard terrestrial solar spectral irradiance curves both at AM 1.5 conditions.14To calculate the AM 0 and AM 1.5 illuminances , the radiometric AM 0 and AM 1.5 spectral irradiances need to be converted to photometric spectral illuminances using the conversion formula in Eq. 2. The photopic response curve of the human eye is shown in Fig. 5 and Fig. 6 shows the spectral illuminance curves for solar AM 0 and AM 1.5 conditions. In Eq. 2, and is the radiometric spectral irradiance. Integrating the AM 0 spectral illuminance curve in Fig. 6 over wavelengths from yields a value of total illuminance that closely matches the maximum value of in Table 1, as expected. Integrating the AM 1.5 Global Tilt and AM 1.5 Direct curves in Fig. 6 yields and , respectively. These values represent illuminances in units of at the entrance aperture of a camera staring directly into the sun. When imaging the sun directly, as shown in Fig. 4 using a focal plane array consisting of , wide dynamic range mesa pixels shown in Fig. 1, the photon flux incident on a pixel and the resulting photocurrent produced by the pixel can be calculated. Table 2 lists important camera parameter values for calculating the photon flux incident on the wide dynamic range APD pixel when imaging the sun directly, as shown in Fig. 4. Table 2Parameters of the APD focal plane array camera.
Using the parameter values in Table 2, the optical power and photon flux incident on a single APD pixel in the array can be calculated for different illuminance conditions at the camera lens aperture. If the camera is orbiting Earth in space and staring directly into the sun at the AM 0 condition, the irradiance at the camera lens will be .13 The total optical power admitted by the lens is given by Eq. 3. In Eq. 3, is the camera aperture area. According to Eq. 3, the total optical power admitted by the lens is calculated to be . For the camera in Fig. 4, imaging the sun directly, the angle and the corresponding solid angle subtended by the sun at the camera lens is given by . The total optical power , projected by the sun into under AM 0 conditions will be focused onto the on the FPA that record the sun’s projected image. The total optical power incident per pixel is therefore given by , however, not all of this incident optical power is actually detected by the mesa APD. As the sun descends toward the horizon, higher airmass conditions apply that produce less illuminance, approaching the lower limit shown in Table 1.Imaging a scene is different from direct imaging of the sun, because the scene is subtended by the field of view of the camera and the scene is projected onto the full light sensitive area of the FPA. Using the illuminance values at the camera aperture provided in Table 1, it is possible to calculate the total integrated photon flux and also the spectral photon flux curves with units [photons/sec-nm] as a function of wavelength that will be incident on a single mesa APD pixel shown in Fig. 1. The spectral photon flux density , incident on the mesa APD for the different scene illuminations, is used together with the previously calculated APD detector quantum efficiency from Ref. 2 to calculate the average APD photocurrent according to Eq. 4.15 The quantity represents the APD detector quantum efficiency and depends on the optical wavelength, while is the average internal avalanche gain of the APD, is the charge of a single electron, and represents the incident photon flux. In Fig. 7 , the average unmultiplied APD photocurrent is calculated for different illumination conditions at the camera lens aperture for the mesa APD pixel in Fig. 1, without a microlens. In Fig. 7, the APD photocurrent corresponding to a total detected photon flux of photons/sec represents the threshold photon flux for enabling Geiger-mode silicon APD operation near room temperature.16 If the mean linear mode photocurrent for the APD operating with unity gain is at or below , as shown in the figure, then Geiger-mode operation of the APD will be supported. This statement is based on the assumption that all photogenerated electrons would be multiplied in the APD to produce detection events or photon counts and also that dark count rates would be much lower than the photogenerated counts by three or more orders of magnitude to not saturate the count rate limit of the APD. The calculation in Fig. 7 shows that Geiger-mode operation of the mesa APD is possible only at night after late dusk. The signal detected by the APD is calculated from the mean photocurrent in Fig. 7, according to Eq. 5. To image with high sensitivity over the normal range of illumination conditions shown in Fig. 7, the APD optical receiver electronics should be capable of supporting linear and Geiger-mode operation. 2.2.Calculated Photocurrent Noise and Dark Current Noise of the Wide Dynamic Range Mesa Avalanche PhotodiodeThe photoelectron and dark current noise are contributed by the mesa APD detector to the overall noise of the pixel or optical receiver, that also includes a contribution from the electronic readout noise of the pixel control circuit. The photoelectron noise can be calculated for the two terminal mesa APD with guard diode shown in Fig. 1, with parameters and describing the probability of a photogenerated electron being collected in the active part of the device and multiplied or collected by the guard diode and not multiplied, respectively. The parameters and are therefore random variables. The instantaneous mesa APD photocurrent with a guard diode structure is given by Eq. 6. Equation 6 expresses the instantaneous APD photocurrent where , , and are all random variables with assumed . The quantity represents the optical power reflectance, is the absorption coefficient of silicon, and is the silicon mesa height. The quantity represents the number of photogenerated electrons produced during integration time interval and describes the avalanche gain. The variables and represent the fractions of that undergo gain multiplication or not, respectively. The mean photocurrent and the variance of the photocurrent for the mesa APD with a guard diode can be calculated. The instantaneous photocurrent in the APD with guard diode is composed of multiplied and unmultiplied photoelectrons and both are expressed by Eqs. 7, 8, respectively.15The total instantaneous photocurrent can be written as the sum of Eqs. 7, 8 as given by Eq. 9.The mean value of the random variable is given by Eq. 10.In Eq. 10, represents the optical power transmittance into the device and represents the absorption efficiency of light in the silicon. The expression for the mean photocurrent can therefore be written as shown in Eq. 11.Equation 11 can be solved in the limit of to yield the mean photocurrent in Eq. 12.To calculate the variance of the photocurrent in the presence of random gain , random , and , the mean square value of the photocurrent must be calculated first according to Eq. 13 where the two currents are uncorrelated.Calculating the mean square value of the photocurrent yields the result in Eq. 14.The variance of the photocurrent in the APD device with a guard diode can now be calculated according to Eq. 15 using Eqs. 12, 14.Since , , and are random variables, a relationship between their mean square values and the square of the means can be defined using excess noise factors as formulated in Eq. 15.6, 17 The excess noise factors for , , , and are provided in Eq. 16.Equation 15 describes the photoelectron noise for the mesa APD structure shown in Fig. 1, with two terminals, a guard diode region and random gain. In Eq. 15, the mean values for and are calculated as a function of the optical wavelength for a particular guard diode boundary function , seen in Fig. 1. The excess noise factors associated with , , and can be determined experimentally in a fabricated device or by Monte Carlo simulation.18In addition to the photoelectron noise contribution by the mesa APD device, the thermally generated dark current in the detector also contributes to the overall noise. The contribution to the detector dark current in the mesa APD structure of Fig. 1 with a guard diode will be composed of multiplied and unmultiplied dark current. The thermally generated dark current in the mesa APD device can be defined and measured under low reverse bias conditions with the APD -junction fully depleted to the silicon-(AlN)-sapphire interface, with a peak electric field below that needed for avalanche gain in the device.19 According to this definition, the mean APD dark current is given by Eq. 17. In Eq. 17, the quantity describes the average number of thermally generated electrons per second and is the charge of an electron. The quantity describes the average number of dark electrons generated during integration time interval . The thermally generated mean dark current can be calculated to yield an estimate for the reverse biased silicon mesa APD -junction. Assuming the dark current is only generated thermally in the silicon bulk, and , where is the silicon intrinsic carrier concentration and only a single type of trap with an energy located inside the bandgap participates; then the recombination rate is described by Eq. 18.20In Eq. 18, the quantities and represent the electron and hole capture cross section probabilities, is the carrier thermal velocity, is the trap density, is the trap energy level, is the intrinsic Fermi level, and is the intrinsic carrier density. Many of the quantities in Eq. 18 depend on temperature as does the recombination rate itself. The thermally generated dark current density from the silicon bulk is calculated by integrating Eq. 18 over the full width of the reverse biased mesa APD -junction depletion region, as given by Eq. 19.In Eq. 19, the quantity represents the mean electron life time in -type silicon material. An estimate of the mean APD dark current from Eq. 19 is obtained by multiplying the result by the cross sectional area of the mesa APD device to yield Eq. 20.The quantity in Eq. 20, depends on the electron capture cross section, carrier thermal velocity, and on the trap concentration. High values for in a silicon -junction can reach .21, 22 Equation 20 assumes that the dark current arises from the silicon bulk and that surface generated dark current contribution to the total dark current is zero. This approximation is valid for well passivated mesa sidewalls. The formalism presented here describes the probabilities for dark electrons to be collected by the active mesa APD versus the guard diode and allows sufficient flexibility for the inclusion of the surface generated dark current contribution if needed.When the reverse bias voltage on the APD is increased to the level of avalanche gain onset, a fraction of will be multiplied and the instantaneous dark current in the APD will be described by the sum of multiplied dark electrons collected in the active silicon and unmultiplied dark electrons collected in the guard diode as given by Eq. 21. In Eq. 21, the quantities , , and are random variables with random values for each collection event. The quantity describes the probability of a thermally generated electron being collected in the active region and therefore being multiplied in the APD, while describes the probability of a thermally generated electron being collected by the guard diode and therefore not being multiplied. The values represent random variables of the gain imparted to a thermal electron generated during the time interval .The random variable is described by Eq. 22. The probability represents the probability of generation of a thermal electron during the time interval . It is assumed as before that the variables are independent and the variables are also independent. The mean dark current can now be expressed using Eqs. 21, 22 according to Eq. 23.Taking the limit , the sums in Eq. 23 become integrals according to Eq. 24.The mean dark current can now be expressed according to Eq. 25.To calculate the dark current variance, the mean square dark current is calculated to yield the result in Eq. 26.The variance of the dark current is given using Eqs. 25, 26 according to Eq. 27.Since and are random variables, a relationship between their mean square values and the square of the means is defined using excess noise factors. The excess noise factors for , , , and are given in Eq. 28.Equation 27 describes the dark current noise for the mesa APD pixel design in Fig. 1 with two terminals, a guard diode region and random gain. In Eq. 27, estimates for the mean values of and can be made by assuming that has a contribution only from the thermally generated bulk dark current. The thermally generated electron flux is calculated from Eqs. 20, 17 and depends on the temperature. For the case of zero surface dark current, represents the fractional volume of active silicon to the total device mesa silicon volume and represents the fractional volume of guard diode silicon to the total mesa silicon volume. The pyramidal mesa frustum volume is given by Eq. 29. In Eq. 29, is the pyramidal frustum base area and is the pyramidal frustum top area. For the mesa of height , with a guard diode boundary defined by parameter with the same slope as the mesa sidewall as shown in Fig. 1, the mean fractional volume of the active and guard silicon is given by Eq. 30.It is assumed for the calculation of the mesa APD pixel that the boundaries between the guard diode region and active region in the device are defined by the fixed boundary in Fig. 1, for which case the excess noise factors and in Eq. 28 are reduced to unity. 3.Pixel Circuit DesignThe pixel sensitivity depends on the electronic readout circuit design which must be capable of controlling linear and Geiger-mode operation of the detector as shown in Fig. 7, to image successfully over the full range of natural illumination conditions. The pixel readout circuit in Fig. 8 , based on a voltage source follower amplifier, may be used to control linear mode operation of the mesa APD and the circuit in Fig. 9 to control Geiger-mode operation. The linear mode electronic readout circuit for the APD in Fig. 8 combines a voltage source follower amplifier consisting of transistor M1 and load resistor , with pass transistor M3 to store the charge integrated at the gate capacitance of M1. The circuit is shown in the figure to be capable of performing correlated double sampling (CDS) using a second matching voltage source follower amplifier made up of transistor M2 with its own identical load resistor and pass transistor M4. The readout circuit in Fig. 9 implements active quenching for the Geiger-mode APD to shorten the detector reset time. During a normal image acquisition cycle, the APD and readout circuit are first reset by turning on then off reset transistor M2. Transistor M1 is then used to precharge the APD and the gate capacitance of M3 to the dc supply voltage. In response to just a single detected photon, the large avalanche gain in the APD produces a photocurrent capable of quickly discharging the self-capacitance of the APD. The capacitance at the gate of M3 is also discharged by the APD photocurrent, thereby turning on M3. The transistor M3 senses the onset of the avalanche in the APD and turns on M4 which diverts a significant fraction of the charge to ground instead of it flowing through the APD junction where it can fill traps and lengthen the quench time of the detector following a detection event. An advantage of the circuits in Figs. 8 and 9 includes being able to combine them into a single compact pixel readout circuit that is capable of controlling both linear and Geiger-mode operation for the wide dynamic range mesa APD detector as shown in Fig. 10 . The circuit in Fig. 10 combines the analog APD readout circuit shown on the right with a digital readout circuit shown on the left. When the APD is biased in linear mode, transistor M7 is turned off and the APD is controlled by the analog circuit on the right. When the APD is biased in Geiger-mode, transistors M3 and M4 are turned off and M7 is turned on, allowing the APD to be controlled by the digital circuit on the left. 3.1.Linear Mode Circuit CharacteristicsDuring a normal image frame acquisition by the linear circuit, both pass transistors M3 and M4 are initially conducting or turned on. Transistor M6 resets the APD as indicated in the pixel control clocking waveforms. Transistor M5 is used to precharge the APD self-capacitance and the gate capacitance of M1 and M2 to the supply voltage. Transistor M4 is allowed to stay turned on briefly after M5 is turned off to allow the voltages at the gates of M1 and M2 to equalize. Following the voltage equalization, M4 is turned off, marking the start of the frame integration time , indicated in Fig. 8. During the integration time a photocurrent , produced by incident illumination on the APD, discharges the precharged self-capacitance of the APD as well as the integration capacitance at the gate of transistor M1. At the end of the charge integration, pass transistor M3 is turned off and the charge is stored on the integration capacitor at the gate terminal of M1 which also serves as the input to the voltage source follower amplifier. Following the image frame integration time, the analog voltage output from the pixel circuit can be digitized using an A/D converter located on the periphery of the pixel array. The linear mode control circuit in Fig. 8 can provide low noise using CDS. In contrast to the digital readout circuit that contributes zero electronic readout noise, the analog circuit composed of the voltage source follower and pass transistor with integration capacitor at the input of the voltage source follower, contributes to the electronic noise. Correlated double sampling shown in Fig. 10 can be used to reduce the noise contribution from the integration capacitor, however, such CDS circuitry may increase the area footprint required for the dual mode readout circuit in Fig. 10, beyond the limit. Therefore, it is more appropriate in the noise analysis to consider the readout circuit for linear mode operation without CDS. Such an analysis is also more conservative because it assumes the worst case noise contribution from the circuit. Figure 11 shows the linear mode readout circuit without CDS used in the signal-to-noise ratio analysis of the APD pixel. The circuit in Fig. 11 is designed for implementation using the TSMC RF-CMOS process, and Table 3 lists important parameter values for the circuit. Table 3Linear mode electronic readout circuit parameters.
In Table 3, the gate length of M1 is given as because the foundry process offers a transistor with a thicker gate oxide and longer gate length. Using the higher voltage transistor in the pixel circuit allows biasing M1 into the active region when the gate capacitance is precharged at the start of signal integration. The RF-CMOS fabrication process supports a metal capacitor with capacitance per unit area of . Such a capacitor is linear and with a total area of , provides of capacitance at the gate of transistor M1. The total capacitance at the gate of M1 is given by Eq. 31. Transistor M1 of the voltage source follower in Fig. 11 can be biased into the active region with , when is precharged by M5 to . The total charge stored on for will therefore be given by Eq. 32.The input voltage at the gate of M1 is given by the initial charge deposited on the integration capacitance corresponding to the precharge voltage , minus the charge that was sourced by the photocurrent flowing through the APD during the integration time interval given by Eq. 33. The voltage transfer characteristic of the source follower amplifier in Fig. 11 can be calculated by knowing that M1 operates in the active region. If the integration capacitor is precharged to a voltage as described in Eq. 32, the voltage transfer characteristic for the amplifier output voltage as a function of can be calculated from Eq. 34 where is calculated iteratively taking into consideration the body effect. The voltage transfer characteristic of the voltage source follower amplifier is plotted in Fig. 12 .3.2.Linear Mode Circuit Noise CharacteristicsThe electronic noise from the linear mode circuit in Fig. 11 is obtained by analyzing the small signal noise models. The resistor in Fig. 13 represents the finite output impedance of M1 including channel length modulation effects. The possible difference in potential between the source and body terminals of M1, gives rise to a voltage-controlled current source with current given by . As a result, it can be modeled as a resistor with value equal to , as shown in Fig. 13. Since the three resistances , , and are in parallel with and since and also , then the noise due to the load resistor will dominate as shown in Fig. 13. The noise from MOSFET M1 is also indicated in Fig. 13. The finite resistance of MOSFET channel M3 contributes to the noise of the sample and hold output shown in Fig. 14 . The noise voltage from the sample and hold circuit in Fig. 14 is contributed at the gate of MOSFET M1 and this can be converted into a noise current variance at the output of the source follower amplifier according to Eq. 35.23 The total noise from the linear APD readout circuit of Fig. 11 is therefore expressed according to Eq. 36.The electronic noise result from Eq. 36 is used in the sensitivity analysis of Sec. 4 for the wide dynamic range mesa APD optical receiver based on the mesa APD of Fig. 1.4.Sensitivity of the Linear and Geiger-Mode Avalanche Photodiode Optical ReceiverThe sensitivity or signal-to-noise of the wide dynamic range mesa APD optical receiver is calculated using the mean photocurrent given in Eq. 4, together with the variance of the photocurrent in Eq. 15, the dark current in Eq. 27 and the electronic circuit readout noise in Eq. 36 according to Eq. 37. Other sources of noise such as the background noise which is a type of photon noise associated with light from extraneous optical sources other than the signal of interest, is not considered in the analysis. Background noise is mainly a problem when imaging at midwave IR and longwave IR wavelengths where objects in a scene emit radiation according to the Planck black body law.17 Leakage currents are also not considered in the noise analysis because the surface leakage currents usually occur at poorly passivated surfaces of semiconductor -junctions, for example, when ions on the semiconductor surface induce mirror charges in the semiconductor that form surface channels, enabling a surface leakage current to flow.20 Bulk leakage currents may occur as a result of defects in the crystal structure of the bulk semiconductor material. In the mesa structure APD design presented here, the surfaces of the mesa are passivated with thermally grown silicon dioxide on the (100) and (111) planes. Therefore, it is expected that surface leakage currents at the APD -junction will be minimized through good quality thermally grown oxide passivation. Surface leakage currents at the (111) mesa sidewalls of the APD are expected to be negligible compared to the bulk dark current in the device. The APD guard diode at the device periphery will also collect and limit the multiplication of any surface leakage currents.19 The signal-to-noise ratio is calculated for the linear mode mesa APD pixel according to Eq. 37 for the normal range of illumination conditions with the following assumptions listed in Table 4 . Table 4Sensitivity calculation parameters.
Using the parameter values listed in Table 4, the signal-to-noise ratio for the wide dynamic range APD optical receiver in linear mode is calculated according to Eq. 37 for different natural illumination conditions. The results are shown in Figs. 15, 16 17, 18 19, 20, 21 for AM 0 to late dusk. In Fig. 15 the APD stares into the sun and unity gain provides SNR above the unity threshold. The maximum integration time is also indicated due to the large APD photocurrent that quickly discharges the integration capacitor of the linear readout circuit. In Fig. 16 the maximum integration is indicated because of the large APD photocurrent that quickly discharges the integration capacitor of the linear readout circuit. In Fig. 17, gain becomes necessary to provide larger than unity signal-to-noise ratio. In Fig. 18, gain becomes necessary to provide larger than unity signal-to-noise ratio. In Fig. 19, gain becomes necessary to provide larger than unity signal-to-noise ratio. In Fig. 20, gain becomes necessary to provide larger than unity signal-to-noise ratio. In Fig. 21, gain becomes necessary to provide larger than unity signal-to-noise ratio. At late dusk, the capability of the linear mode imager to provide high sensitivity for imaging becomes limited, especially at larger bandwidths. After late dusk, however, the Geiger-mode operation of the mesa APD detector without a microlens is supported as calculated in Fig. 7. The signal-to-noise ratio for the Geiger-mode mesa APD optical receiver that uses the APD in Fig. 1 and the electronic readout circuit in Fig. 10 is calculated using the expression given by Eq. 38. In Eq. 38, the digital pixel control circuit does not add to the optical receiver noise.24, 25 This enabling characteristic of the Geiger-mode pixel control circuit allows the optical receiver to operate with high sensitivity at the low light levels indicated in Fig. 7. Using the parameter values listed in Table 4, without the entries for the linear circuit, the signal-to-noise ratio for the wide dynamic range APD optical receiver in Geiger-mode is calculated according to Eq. 38. The results are shown in Figs. 22, 23, 24, 25, 26 for late dusk to a cloudy night with no moon.In Fig. 22, the mesa APD pixel is shown to provide high sensitivity in the Geiger-mode. The SNR for and with resolution is with . In Fig. 23, the SNR for and with resolution is with . In Fig. 24, the SNR for and with resolution is with . In Fig. 25, the SNR for and with resolution is with . In Fig. 26, the SNR for and with resolution is with . The calculations of signal-to-noise ratio for the wide dynamic range mesa APD optical receiver operating in Geiger-mode show that the imager will be capable of recording an image even under the lowest natural illumination conditions in Fig. 26. This result is important because it confirms that the detector array design should be capable of imaging successfully under the broadest range of natural illuminations when cooled only to using a compact two-stage thermoelectric cooler. 5.Concluding RemarksIt is shown in this detailed design study that a novel back-illuminated silicon APD array architecture using (100) epitaxial silicon grown on -plane sapphire having a lattice-matched AlN antireflective layer provides a unique means for creating high sensitivity, wide dynamic range APD device arrays for large area, high resolution, solid-state APD-FPA imagers with zero optical crosstalk and 100% sensitive-area-fill-factor. Such APD-FPA imagers could be used for both passive and active 2-D and 3-D imaging applications. The detailed calculations and analysis presented here show that the APD-FPA imager design will be capable of operating with high signal-to-noise ratio over the full range of natural illumination conditions near room temperature. Such a camera represents a novel, enabling technology for scientific and industrial applications. referencesA. G. Stern and D. C. Cole,
“High quantum efficiency, back-illuminated, crystallographically etched, silicon-on-sapphire avalanche photodiode with very wide dynamic range, for manufacturable high resolution imaging arrays,”
Proc. SPIE, 7249 72490R
(2009). https://doi.org/10.1117/12.806272 Google Scholar
A. G. Stern and D. C. Cole,
“Design of a back-illuminated, crystallographically etched, silicon-on-sapphire avalanche photodiode with monolithically integrated microlens, for dual-mode passive and active arrays,”
Proc. SPIE, 7153 71530Z
(2008). https://doi.org/10.1117/12.804846 Google Scholar
I. Rech, A. Ingargiola, R. Spinelli, I. Labanca, S. Marangoni, M. Ghioni, and S. Cova,
“A new approach to optical crosstalk modelling in single-photon avalanche diodes,”
IEEE Photonics Technol. Lett., 20
(50), 330
–332
(2008). https://doi.org/10.1109/LPT.2007.915654 Google Scholar
S. Cihangir and S. Kwan,
“Characterization of indium and solder bump bonding for pixel detectors,”
Nucl. Instrum. Methods Phys. Res. A, 476
(3), 670
–675
(2002). https://doi.org/10.1016/S0168-9002(01)01771-5 Google Scholar
R. J. McIntyre,
“Multiplication noise in uniform avalanche diodes,”
IEEE Trans. Electron Devices, 13
(1), 164
–168
(1966). https://doi.org/10.1109/T-ED.1966.15651 Google Scholar
R. J. McIntyre,
“A new look at impact ionization—part I: a theory of gain, noise, breakdown probability, and frequency response,”
IEEE Trans. Electron Devices, 46
(8), 1623
–1631
(1999). https://doi.org/10.1109/16.777150 Google Scholar
P. Yuan, K. A. Anselm, C. Hu, H. Nie, C. Lenox, A. L. Holmes, B. G. Streetman, J. C. Campbell, and R. J. McIntyre,
“A new look at impact ionization—part II: gain and noise in short avalanche photodiodes,”
IEEE Trans. Electron Devices, 46
(8), 1632
–1639
(1999). https://doi.org/10.1109/16.777151 Google Scholar
T. Shibata, K. Asai, Y. Nakamura, M. Tanaka, K. Kaigawa, J. Shibata, and H. Sakai,
“AlN epitaxial growth on off-angle R-plane sapphire substrates by MOCVD,”
J. Cryst. Growth, 229
(1), 63
–68
(2001). https://doi.org/10.1016/S0022-0248(01)01051-X Google Scholar
W. Luo, X. Wang, L. Guo, H. Xiao, C. Wang, J. Ran, J. Li, and J. Li,
“Influence of AlN buffer layer thickness on the properties of GaN epilayer on Si(111) by MOCVD,”
Microelectron. J., 39
(12), 1710
–1713
(2008). https://doi.org/10.1016/j.mejo.2008.01.042 Google Scholar
Ü. Ozgür, G. Webb-Wood, H. O. Everitt, F. Yun, and H. Morkoc,
“Systematic measurement of
refractive indices,”
Appl. Phys. Lett., 79
(25), 4103
–4105
(2001). https://doi.org/10.1063/1.1426270 Google Scholar
C. K. Kim,
“The physics of charge-coupled devices,”
Charge-Coupled Devices and Systems, 1
–80 Wiley, New York
(1979). Google Scholar
S. D. Donati, Photodetectors, Devices, Circuits and Applications, 1
–36 Prentice Hall, Upper Saddle River, NJ
(2000). Google Scholar
C. Wehrli,
“Extraterrestrial solar spectrum,”
Physikalisch-Meteorologisches Observatorium-World Radiation Center (PMO/WRC),
(1985) Google Scholar
L. W. Thomason, B. M. Herman, and J. A. Reagan,
“The effect of atmospheric attenuators with structured vertical distributions on air mass determinations and Langley plot analyses,”
J. Atmos. Sci., 40 1851
–1854
(1983). https://doi.org/10.1175/1520-0469(1983)040<1851:TEOAAW>2.0.CO;2 Google Scholar
B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics, 246755789
–253 John Wiley and Sons, New York
(2007). Google Scholar
P. Kumar, P. Kwiat, A. Migdall, S. W. Nam, J. Vuckovic, and N. C. Wong,
“Photonic technologies for quantum information processing,”
Quantum Inf. Process., 3
(1–5), 215
–231
(2004). https://doi.org/10.1007/s11128-004-3102-4 Google Scholar
G. E. Stillman and C. M. Wolfe,
“Avalanche photodiodes,”
Semiconductors and Semimetals, 291
–391 Academic Press, New York
(1977). Google Scholar
R. J. McIntyre,
“The distribution of gains in uniformly multiplying avalanche photodiodes: theory,”
IEEE Trans. Electron Devices, 19
(6), 703
–713
(1972). https://doi.org/10.1109/T-ED.1972.17485 Google Scholar
J. C. Jackson, P. K. Hurley, B. Lane, and A. Mathewson,
“Comparing leakage currents and dark count rates in Geiger-mode avalanche photodiodes,”
Appl. Phys. Lett., 80
(22), 4100
(2002). https://doi.org/10.1063/1.1483119 Google Scholar
S. M. Sze, Physics of Semiconductor Devices, 495412
–427 2nd ed.John Wiley and Sons, New York
(1981). Google Scholar
A. Rohatgi and P. R. Choudhury,
“Process induced effects on carrier lifetime and defects in float-zone silicon,”
J. Electrochem. Soc., 127
(5), 1136
–1139
(1980). https://doi.org/10.1149/1.2129833 Google Scholar
W. Zimmerman,
“Experimental verification of the Shockley-Read Hall recombination theory in silicon,”
Electron. Lett., 9
(16), 378
–379
(1973). https://doi.org/10.1049/el:19730279 Google Scholar
D. A. Johns and K. Martin, Analog Integrated Circuit Design, 181
–217 John Wiley and Sons, New York
(1996). Google Scholar
M. Ghioni, S. Cova, F. Zappa, and C. Samori,
“Compact active quenching circuit for fast counting with avalanche photodiodes,”
Rev. Sci. Instrum., 67
(10), 3440
–3448
(1996). https://doi.org/10.1063/1.1147156 Google Scholar
F. Zappa, M. Ghioni, S. Cova, C. Samori, and A. C. Giudice,
“An integrated active-quenching circuit for single-photon avalanche diodes,”
IEEE Trans. Instrum. Meas., 49
(6), 1167
–1175
(2000). https://doi.org/10.1109/19.893251 Google Scholar
BiographyAlvin G. Stern recently completed his doctorate in electrical engineering at Boston University. He completed his BS degree in electrical engineering (BSEE) also at Boston University in May of 2000. He was a member of the technical staff in the Advanced Imaging Technology Group at MIT Lincoln Laboratory for four years prior to returning to Boston University in 2004 to pursue PhD studies. He completed an MSEE in 2007 in the area of superconducting single photon detector (SSPD) characterization. Research interests include the design through analysis and simulation means of solid-state silicon avalanche photodiode (APD) devices and their arrays for imaging applications. Research interests further include the design of high-speed electronic readout circuits to control APD detector array operation. Daniel C. Cole is an associate professor in the Department of Mechanical Engineering at Boston University (BU). He obtained his doctorate in theoretical physics from the Graduate Center of the City University of New York in 1985. He worked for IBM Microelectronics in the area of physical simulation of semiconductor processes and microelectronic devices until joining BU in 1998. During this time he also wrote a number of articles involving the thermodynamics of charged particles and electromagnetic radiation. During the last three years at IBM, he managed a department of physicists, engineers, and mathematicians involved with developing and applying large computational programs to semiconductor problems. His present research work is split between the simulation of physical processes to aid the development of nanoelectronic devices, and the further investigation of an area of physics often called stochastic electrodynamics. |