Fully integrated free-running InGaAs/InP single-photon detector for accurate lidar applications

TL;DR

This paper introduces a fully integrated free-running InGaAs/InP single-photon detector optimized for lidar, demonstrating comparable performance to superconducting detectors after correction algorithms, with advantages in cost and size.

Contribution

The paper presents a novel integrated InGaAs/InP NFAD-based detector with optimized parameters and correction algorithms, suitable for practical, accurate lidar applications.

Findings

Achieved 1.6 Mcps maximum count rate with 10% PDE and 950 cps dark count rate.

Performance after correction closely matches that of superconducting detectors, with only ~2% error.

Demonstrated practical lidar application with comparable results to advanced detectors.

Abstract

We present a fully integrated InGaAs/InP negative feedback avalanche diode (NFAD) based free-running single-photon detector (SPD) designed for accurate lidar applications. A free-piston Stirling cooler is used to cool down the NFAD with a large temperature range, and an active hold-off circuit implemented in a field programmable gate array is applied to further suppress the afterpulsing contribution. The key parameters of the free-running SPD including photon detection efficiency (PDE), dark count rate (DCR), afterpulse probability, and maximum count rate (MCR) are dedicatedly optimized for lidar application in practice. We then perform a field experiment using a Mie lidar system with 20 kHz pulse repetition frequency to compare the performance between the free-running InGaAs/InP SPD and a commercial superconducting nanowire single-photon detector (SNSPD). Our detector exhibits good performance with 1.6 Mcps MCR (0.6 {\mu}s hold-off time), 10% PDE, 950 cps DCR, and 18% afterpulse probability over 50 {\mu}s period. Such performance is worse than the SNSPD with 60% PDE and 300 cps DCR. However, after performing a specific algorithm that we have developed for afterpulse and count rate corrections, the lidar system performance in terms of range-corrected signal (Pr2) distribution using our SPD agrees very well with the result using the SNSPD, with only a relative error of ~2%. Due to the advantages of low-cost and small size of InGaAs/InP NFADs, such detector provides a practical solution for accurate lidar applications.

Tables

Table 1. Table 1: The calibration results of dark count rate and afterpulse probability at different temperatures (T) with 10% PDE.

T (K)	DCR (cps)	$P_{ap}$ (%)
243	6743	13.2
223	952	18.0
203	135	38.5

Equations

$$I\eta =\frac{R_s}{1-R_s\tau }-DCR.$$

$$Er{r}_{os}=\frac{R_o-R_s}{R_s}.$$

$$r_d=\frac{R}{1-R\tau }.$$

$$PDE=-\frac{1}{\mu }\ln (1-\frac{S_1}{f}),$$

$$P_{ap}=\frac{S_2-S_1-r_d(1-S_2\tau )}{S_1}.$$

$$P_{ap}(t)=A_{1}Exp(-t/{\tau }_1)+A_{2}Exp(-t/{\tau }_2)+A_{3}Exp(-t/{\tau }_3)+c,$$

$$R_{ap}(j)=\sum _{i}A(i,j),$$

$$A(i,j)=R(i)P_{nc}(i,j)P_{nap}(i,j)P_{ap}(j;i),$$

$$P_{nc}(i,j)=Exp[-\sum _{k=i}^{j}R(k)bi{n}_w],$$

$$P_{nap}(i,j)=Exp[-\sum _{k=0}^{j-i-1}{P}_{ap}(k)],$$

$$R_1(i)=R(i)-R_{ap}(i),$$

$$R_2(i)=\frac{R_1(i)}{1-R(i)\tau }-DCR(i).$$

Full text

OCIS codes: (010.3640) Lidar; (010.0280) Remote sensing and sensors; (040.1345) Avalanche photodiodes (APDs); (040.5570) Quantum detectors.

References and links

• [1]

S. A. Lisenko, M. M. Kugeiko, and V. V. Khomich, “Multifrequency lidar sounding of air pollution by particulate matter with separation into respirable fractions,” Atmos. Ocean. Opt. 29, 288 (2016).

• [2]

M. D. Eisaman, J. Fan, A. Migdall, and S. V. Polyakov, “Invited review article: Single-photon sources and detectors,” Rev. Sci. Instrum. 82, 071101 (2011).

• [3]

F. Marsili, V. B. Verma, J. A. Stern, S. Harrington, A. E. Lita, T. Gerrits, I. Vayshenker, B. Baek, M. D. Shaw, R. P. Mirin, and S. W. Nam, “Detecting single infrared photons with 93% system efficiency,” Nat. Photonics 7, 210-214 (2013).

• [4]

X. Yang, H. Li, W. Zhang, L. You, L. Zhang, X. Liu, Z. Wang, W. Peng, X. Xie, and M. Jiang, “Superconducting nanowire single photon detector with on-chip bandpass filter,” Opt. Express 22, 16267-16272 (2014).

• [5]

J. Wu, L. You, S. Chen, H. Li, Y. He, C. Lv, Z. Wang, and X. Xie, “Improving the timing jitter of a superconducting nanowire single-photon detection system,” Appl.  Opt. 56, 2195-2200 (2017).

• [6]

H. Xia, G. Shentu, M. Shangguan, X. Xia, X. Jia, C. Wang, J. Zhang, J. S. Pelc, M. M. Fejer, Q. Zhang, X. Dou, and J. W. Pan, “Long-range micro-pulse aerosol lidar at 1.5 $\mu$m with an upconversion single-photon detector,” Opt. Lett. 40, 1579-1582 (2015).

• [7]

L. Hogstedt, A. Fix, M. Wirth, C. Pedersen, and P. Tidemand-Lichtenberg, “Upconversion-based lidar measurements of atmospheric $CO_{2}$,” Opt. Express 24, 5152-5161 (2016).

• [8]

M. Shangguan, H. Xia, C. Wang, J. Qiu, G. Shentu, Q. Zhang, X. Dou, and J. W. Pan, “All-fiber upconversion high spectral resolution wind lidar using a Fabry-Perot interferometer,” Opt. Express 24, 19322-19336 (2016)

• [9]

H. Xia, M. Shangguan, C. Wang, G. Shentu, J. Qiu, Q. Zhang, X. Dou, and J. W. Pan, “Micro-pulse upconversion Doppler lidar for wind and visibility detection in the atmospheric boundary layer,” Opt. Lett. 41, 5218-5221 (2016).

• [10]

M. A. Itzler, X. Jiang, M. Entwistle, K. Slomkowski, A. Tosi, F. Acerbi, F. Zappa, and S. Cova, “Advances in InGaAsP-based avalanche diode single photon detectors,” J. Mod. Opt. 58, 174-200 (2011).

• [11]

J. Zhang, M. A. Itzler, H. Zbinden, and J. W. Pan, “Advances in InGaAs/InP single-photon detector systems for quantum communication,” Light Sci. Appl. 4, e286 (2015).

• [12]

N. Namekata, S. Sasamori, and S. Inoue, “800 MHz Single-photon detection at 1550-nm using an InGaAs/InP avalanche photodiode operated with a sine wave gating,” Opt. Express 14, 10043-10049 (2006).

• [13]

J. Zhang, R. Thew, C. Barreiro, and H. Zbinden, “Practical fast gate rate InGaAs/InP single-photon avalanche photodiodes,” Appl. Phys. Lett. 95, 091103 (2009).

• [14]

N. Namekata, S. Adachi, and S. Inoue, “1.5 GHz single-photon detection at telecommunication wavelengths using sinusoidally gated InGaAs/InP avalanche photodiode,” Opt. Express 17, 6275-6282 (2009).

• [15]

J. Zhang, P. Eraerds, N. Walenta, C. Barreiro, R. Thew, and H. Zbinden, “2.23 GHz gating InGaAs/InP single-photon avalanche diode for quantum key distribution,” Proc. SPIE 7681, 76810Z (2010).

• [16]

Y. Nambu, S. Takahashi, K. Yoshino, A. Tanaka, M. Fujiwara, M. Sasaki, A. Tajima, S. Yorozu, and A. Tomita, “Efficient and low-noise single-photon avalanche photodiode for 1.244-GHz clocked quantum key distribution,” Opt. Express 19, 20531-20541 (2011).

• [17]

X. L. Liang, J. H. Liu, Q. Wang, D. B. Du, J. Ma, G. Jin, Z. B. Chen, J. Zhang, and J. W. Pan, “Fully integrated InGaAs/InP single-photon detector module with gigahertz sine wave gating,” Rev. Sci. Instrum. 83, 083111 (2012).

• [18]

N. Walenta, T. Lunghi, O. Guinnard, R. Houlmann, H. Zbinden, and N. Gisin, “Sine gating detector with simple filtering for low-noise infra-red single photon detection at room temperature,” J. Appl. Phys. 112, 063106 (2012).

• [19]

A. Restelli, J. C. Bienfang, and A. L. Migdall, “Single-photon detection efficiency up to 50% at 1310nm with an InGaAs/InP avalanche diode gated at 1.25 GHz,” Appl. Phys. Lett. 102, 141104 (2013).

• [20]

Z. L. Yuan, B. E. Kardynal, A. W. Sharpe, and A. J. Shields, “High speed single photon detection in the near-infrared,” Appl. Phys. Lett. 91, 041114 (2007).

• [21]

X. Chen, E. Wu, G. Wu, and H. Zeng, “Low-noise high-speed InGaAs/InP-based single-photon detector,” Opt. Express 18, 7010-7018 (2010).

• [22]

P. Eraerds , M. Legre, J. Zhang, H. Zbinden, and N. Gisin, “Photon counting OTDR: advantages and limitations,” J. Lightwave Technol. 28, 952-964 (2010).

• [23]

M. Ren, X. Gu, Y. Liang, W. Kong, E. Wu, G. Wu, and H. Zeng, “Laser ranging at 1550 nm with 1-GHz sine-wave gated InGaAs/InP APD single-photon detector,” Opt. Express 19, 13497-13502 (2011).

• [24]

H. Xia, D. Sun, Y. Yang, F. Shen, J. Dong, and T. Kobayashi, “Fabry-Perot interferometer based Mie Doppler lidar for low tropospheric wind observation,” Appl.  Opt. 46, 7120-7131 (2007).

• [25]

H. Xia, X. Dou, D. Sun, Z. Shu, X. Xue, Y. Han, D. Hu, Y. Han, and T. Cheng, “Mid-altitude wind measurements with mobile Rayleigh Doppler lidar incorporating system-level optical frequency control method,” Opt. Express 20, 15286-15300 (2012).

• [26]

H. Xia, X. Dou, M. Shangguan, R. Zhao, D. Sun, C. Wang, J. Qiu, Z. Shu, X. Xue, Y. Han, and Y. Han, “Stratospheric temperature measurement with scanning Fabry-Perot interferometer for wind retrieval from mobile Rayleigh Doppler lidar,” Opt. Express 22, 21775-21789 (2014).

• [27]

H. Xia, M. Shangguan, G. Shentu, C. Wang, J. Qiu, M. Zheng, X. Xie, X. Dou, Q. Zhang, and J. W. Pan, “Brillouin optical time-domain reflectometry using up-conversion single-photon detector,” Opt. Commun. 381, 37-42 (2016).

• [28]

J. Rarity, T. Wall, K. Ridley, P. Owens, and P. Tapster, “Single-photon counting for the 1300-1600-nm range by use of Peltier-cooled and passively quenched InGaAs avalanche photodiodes,” Appl.  Opt. 39, 6746-6753 (2000).

• [29]

R. E. Warburton, M. Itzler, and G. S. Buller, “Free-running, room temperature operation of an InGaAs/InP single-photon avalanche diode,” Appl. Phys. Lett. 94, 071116 (2009).

• [30]

M. Liu, C. Hu, J. C. Campbell, Z. Pan, and M. M. Tashima, “Reduce afterpulsing of single photon avalanche diodes using passive quenching with active reset,” IEEE J. Quantum Electron. 44, 430-434 (2008).

• [31]

C. Hu , M. Liu , X. Zheng and J. Campbell, “Dynamic range of passive quenching active reset circuit for single photon avalanche diodes,” IEEE J. Quantum Electron. 46, 35-39 (2010).

• [32]

R. T. Thew, D. Stucki, J.-D. Gautier, H. Zbinden, and A. Rochas, “Free-running InGaAs/InP avalanche photodiode with active quenching for single photon counting at telecom wavelengths,” Appl. Phys. Lett. 91, 201114 (2007).

• [33]

J. Zhang, R. Thew, J.-D. Gautier, N. Gisin, and H. Zbinden, “Comprehensive characterization of InGaAsP-InP avalanche photodiodes at 1550 nm with an active quenching ASIC,” IEEE J. Quantum Electron. 45, 792-799 (2009).

• [34]

M. A. Itzler, X. Jiang, B. Nymann, and K. Slomkowski, “InP-based negative feedback avalanche diodes,” Proc. SPIE, 7222, 72221K (2009).

• [35]

M. A. Itzler, X. Jiang, B. M. Onat, and K. Slomkowski, “Progress in self-quenching InP-based single photon detectors,” Proc. SPIE , 7608, 760829 (2010).

• [36]

X. Jiang, M. A. Itzler, K. O’Donnell, M. Entwistle, and K. Slomkowski, “InGaAs/InP negative feedback avalanche diodes,” Proc. SPIE , 8033, 8033K (2011).

• [37]

Z. Yan, D. R. Hamel, A. K. Heinrichs, X. Jiang, M. A. Itzler, and T. Jennewein, “An ultra low noise telecom wavelength free running single photon detector using negative feedback avalanche diode,” Rev. Sci. Instrum. 83, 073105 (2012).

• [38]

B. Korzh, N. Walenta, T. Lunghi, N. Gisin, and H. Zbinden, “Free-running InGaAs single photon detector with 1 dark count per second at 10% efficiency,” Appl. Phys. Lett. 104, 081108 (2014).

• [39]

http://www.idquantique.com/photon-counting/photon-counting-modules/id220/

1 Introduction

Aerosol lidars play an important role in many fields, such as atmosphere environment management, agricultural meteorology, hydrological cycle and radiation budget. Compared with ultraviolet (UV) and Visible systems, 1.5 $\mu$m aerosol lidars offer several advantages, including higher maximum permissible exposure to human eyes, lower atmospheric attenuation, minor disturbance from Rayleigh scattering, and weaker sky radiance. For instance, when using a multifrequency lidar for retrieving the spatial distribution of respirable fractions of aerosol in the lower atmosphere, 1.5 $\mu$m aerosol lidars provide more reliable concentration estimates of coarse aerosol particles [1].

Using single-photon detectors (SPDs) at 1.5 $\mu$m can greatly improve the performance of aerosol lidars. Currently the main techniques for single-photon detection at the telecom wavelength [2] include superconducting nanowire single-photon detectors (SNSPDs), up-conversion SPDs and InGaAs/InP single-photon avalanche diodes (SPADs). SNSPDs exhibit excellent performance of high photon detection efficiency (PDE), low dark count rate (DCR) and low timing jitter [3, 4, 5], however, the requirement of cryogenic conditions limits the use for practical applications. Up-conversion SPDs exhibit moderate performance, and recently have been used for lidar experiments [6, 7, 8, 9]. Compared with the above detectors, InGaAs/InP SPADs have relatively poor performance, however, such detectors are widely used in practical applications due to their advantages of small size and low cost [10, 11]. The key parameters of InGaAs/InP SPADs, i.e., PDE, DCR, afterpulse probability and maximum count rate (MCR), influence each other, so that these parameters should be compromised and optimized according to the requirements of applications. Apart from SPAD device, the quenching electronics is the crucial part for a SPAD [11]. InGaAs/InP SPADs are often operated either in gating mode or in free-running mode, which are suited for synchronous and asynchronous single-photon detections, respectively.

Gating mode can effectively reduce DCR of InGaAs/InP SPADs, and gating frequencies have been increased from a few MHz in the early stage to the regime of GHz [11]. High-frequency gating techniques, including sine wave gating [12, 13, 14, 15, 16, 17, 18, 19] and self-differencing [20, 21], can highly shorten avalanche duration time and thus significantly suppress the afterpulsing effect, which greatly helps to increase count rates of InGaAs/InP SPADs. However, gating mode is not well suited for remote sensing applications due to the limit of small duty cycle for asynchronous single-photon detection [22, 23]. Even using high-frequency gating technique, the duty cycle can reach around 20%, which still results in pretty low effective detection efficiency [11].

For lidar and remote sensing applications [24, 25, 26, 27], the most practical solution is using free-running InGaAs/InP SPADs. So far several techniques have been used to implement free-running operations. First, passive quenching is certainly a common approach for free-running mode. Rarity et al. demonstrated a passively quenched InGaAs/InP SPD in 2000 [28]. The important technical issue to be settled in such detectors is to suppress afterpulsing effect [11]. Using long recovery time can reduce the afterpulse probability. However, the performance of count rate is severely limited. Warburton et al. presented a solution for free-running operation through lowering excess bias and increasing temperature [29]. The Virginia group demonstrated a scheme called passive quenching and active reset (PQAR) [30, 31], combining the advantages of parasitic capacitance minimization by chip-to-chip wire bonding and fast electronic switch to swiftly reset the bias voltage after the hold-off time. Second, active quenching can also be used for free-running mode given a short quenching time. The Geneva group demonstrated a practical free-running detector using integrated circuit of active quenching, in which both the parasitic capacitance and the propagation delay of the quenching circuit were minimized [32, 33]. Third, new SPAD devices called negative feedback avalanche diodes (NFADs) provide a simple and practical solution for free-running mode [34, 35, 36]. NFADs monolithically integrate thin film resisters ($\sim$ 900 K$\Omega$) inside the semiconductor structure of devices for passive quenching operation, which can fundamentally solve the problem of parasitic capacitance compared with the PQAR scheme. Yan et al. implemented a NFAD based free-running detector with 10% PDE and 100 cps DCR at 193 K [37]. Korzh et al. demonstrated a similar detector and optimized the operation conditions to further reduce DCR down to 1 cps at 10% PDE [38]. However, these two free-running SPDs used pretty long hold-off time to further suppress the afterpulse probability, which results in low MCR and even signal distortion in lidar applications.

In this paper, we present a fully integrated free-running SPD based on InGaAs/InP NFAD (Princeton Lightwave) with 1.6 Mcps MCR, 10% PDE and 950 cps DCR. We then develop a specific algorithm for afterpulse correction and count rate correction for lidar applications. After the corrections, the lidar system performance in terms of range-corrected signal ($Pr^{2}$) distribution using such detector is considerably comparable to that of using a commercial SNSPD with only a relative error of $\sim$2%.

2 Mie lidar at 1.5 $\mu$m

The experimental setup of 1.5 $\mu$m Mie lidar using free-running InGaAs/InP SPD is shown in Fig. 1. The whole laser system of lidar utilizes master-oscillator power-amplifier architecture. A continuous wave (CW) laser from a distributed feedback diode (DFB, 1548.1 nm) is chopped into a pulse train using an electro-optic modulator (EOM, Photline, MXER-LN-10) with high extinction ratio (35 dB). The EOM is driven by a pulse generator, which controls the shape and pulse repetition frequency (PRF) of the laser. In the experiment, the PRF of the pulses is set to 20 kHz, which indicates the maximum unambiguous detection range is $\sim$7.5 km. The weak laser pulses are fed into an erbium-doped fiber amplifier (EDFA, Keyopsys, PEFA-EOLA), which emits a pulse train with 110 $\mu$J pulse energy and 200 ns pulse width. A large-mode-area fiber with numerical aperture of 0.08 is used to increase the threshold of stimulated Brillouin scattering and to avoid self-saturation of amplified spontaneous emission (ASE). The laser is then collimated and sent to the atmosphere.

The backscattering signal from the atmosphere is collected into a single-mode fiber using a pigtailed coupler. The background noise is filtered out by inserting a fiber Bragg grating (FBG, Advanced Optics Solution GmbH) with a ultra-narrow bandwidth of 6 pm and an insertion loss less than 3 dB. The backscattering signal is detected by an InGaAs/InP NFAD based free-running SPD. The detection signals of the received photons are recorded on a multi-channel scaler (FAST ComTec, MCS6A) and then transmitted to a computer. In order to optimize the performance of free-running InGaAs/InP SPD for lidar applications, a commercial superconductor nanowire single photon detector (SNSPD, Single Quantum, Eos 210 CS) is used to detect the atmospheric backscattering signal simultaneously for comparison. The SNSPD exhibits pretty good calibration performance with 60% PDE, 300 cps DCR, low timing jitter and high MCR. Therefore, such detector is well suited to be used as a reference in the experiment.

3 Free-running InGaAs/InP single-photon detector

Backscattering signal from a lidar decays very fast along detection range so that the count rate of signal has a large dynamic range. Therefore, a SPD with high MCR is required to avoid detector saturation, particularly in the near range. The main design considerations of our free-running InGaAs/InP detector system include moderate DCR and afterpulse performance with high MCR and integrated module for practical use. The design diagram and photo of the detector are shown in Figs. 2(a) and 2(b), respectively. To achieve low DCR performance, cooling down the NFAD device is crucial. In the detector, a free-piston Stirling cooler (FPSC, SC-UE15R) is used, which has a powerful cooling capacity and can easily cool down the NFAD to 163 K. The NFAD is fixed on the cold side surface of FPSC and then the whole cold side is encapsulated inside a cavity. The cavity is connected with the detector controller via a micro rectangular connector for temperature control and micro coaxial (MCX) connectors for electronic signal transmission. A target temperature is set by rotating a potentiometer and displayed on the detector panel. Since the cooling power of FPSC is tuned by a direct current (DC) voltage in its driver circuit, a proportional-integral-derivative (PID) program is developed in a field-programmable gate array (FPGA) to regulate the voltage in real-time. In such a way, the cooling temperature that is displayed on the detector panel in real-time can be finally stabilized at the target temperature with an error less than 1 K.

The reverse bias voltage of NFAD is tuned by another potentiometer on the detector panel to guarantee that the NFAD is operated in Geiger mode. Once an avalanche occurs, the original avalanche signal is alternating current (AC) coupled to a low-noise amplifier (LNA) with a gain of 40 dB. The amplified signal is then discriminated. The output signals of the discriminator are sent to the FPGA for further processing. The FPGA counts the discriminated signals and drivers the displayer on the panel, and outputs reset pulses as well. The reset LVTTL pulses are amplified to 5 V via an operational amplifier (OPA) and then coupled to the anode of NFAD. When the reset pulses are at HIGH level, the bias is below breakdown voltage and hence the width of the reset pulses are exactly the hold-off time. Once the status of the reset pulses is switched from HIGH level to LOW level, the NFAD is directly reset to the initial state. Due to the capacitance response of NFAD, the falling edges of reset pulses produce negative derivative signals that are superposed with avalanche signals. In order to avoid wrong detections, the output of the discriminator is latched until 20 ns after the falling edges of reset pulses so that negative derivative signals are ignored.

4 Optimization of parameters

Afterpulse probability is an additional parameter to be considered in lidar applications compared with SNSPD. This parameter is related to the operation conditions of NFAD including excess bias voltage, temperature, and hold-off time. Setting a long hold-off time can effectively reduce aftepulse probability, which, however, considerably limits the parameter of MCR and results in a ringing effect to the detected backscattering signal in lidar applications. Figures 3(a) and 3(b) show the backscattering signals detected by the free-running InGaAs/InP SPD with different settings of hold-off time. In the case of 2 $\mu$s hold-off time one can clearly see the ringing effect of detector count rate in the peak area while the ringing disappears in the case of 0.6 $\mu$s hold-off time.

In order to illustrate the ringing effect, we develop a MATLAB program to simulate the detector count rate of free-running InGaAs/InP SPD with different settings of hold-off time. The input intensity is assumed to be a step signal from 0 to $10^{7}$ photons with a rising time of 2 $\mu$s, which are close to the values in the field experiment. Given a setting of 2 $\mu$s hold-off time, the ringing effect distinctly appears, as shown in Fig. 3(c), from which one can conclude that such effect highly relies on the input photon flux and the hold-off time setting of SPD. Further, let us define the maximum count rate as $R_{o}$ and the stable count rate as $R_{s}$, the relationship between $R_{s}$ and the input intensity is given by

[TABLE]

where $\eta$ is the detection efficiency, $\tau$ is the hold-off time, and $I$ is the input photon flux. The maximum relative overshoot error is calculated as

[TABLE]

Figure 3(d) plots $Err_{os}$ as a function of hold-off time, which indicates that such overshoot error exponentially increases with the increase of $\tau$. Therefore, for lidar applications using InGaAs/InP SPD the hold-off time has to be chosen as small as possible to avoid the count rate ringing. In our experiment, the hold-off time is set as 0.6 $\mu$s, corresponding to a MCR of $\sim$1.6 Mcps, to guarantee that the maximum relative overshoot error is less than 1.5%.

We then calibrate the key parameters of the free-running InGaAs/InP SPD with 0.6 $\mu$s hold-off time including PDE, DCR, and afterpulse probability distribution. The setup for calibration is shown in Fig. 4. A laser diode (LD, picoquant PDL 800-D) is driven by a pulse generator (PG) with a frequency of 20 kHz that is the same as the PRF of our lidar system. The laser pulses are split by a beam splitter (BS). One output port of the BS is monitored by a power meter (PM, Exfo IQS-1600), and the other one is connected with a precise attenuator (ATT) to further attenuate the intensity down to single-photon level. The detector outputs are recorded by a time-to-digital converter (TDC, picoquant PicoHarp300), which is triggered by synchronized signals from the PG.

For DCR measurement, the LD is switched off. Given a count rate $R$, the inherent dark count rate $r$ is calculated as

[TABLE]

When the LD is switched on, the peak photon count rate $S_{1}$ and total count rate $S_{2}$ are recorded, respectively. Then PDE is calculated as

[TABLE]

where $\mu$ is the mean photon number per pulse, $f$ is the frequency of laser pulses. The distribution of afterpulse counts can be obtained from the TDC, and the total afterpulse probability over the whole 50 $\mu$s period is evaluated by

[TABLE]

Table 1 shows the calibration results of DCR and afterpulse probability with 10% PDE at three typical temperatures with hold-off time of 0.6 $\mu$s. At low temperatures, the SPD exhibits low DCR performance but high afterpulse probability. Therefore, these two parameters have to be compromised for lidar applications. In the lidar experiment, the temperature of the SPD is set as 223 K with performance of 10% PDE, 950 cps DCR and 18% afterpulse probability. Such performance is moderately better than the commercial product of free-running InGaAs/InP SPD (ID220-FR-SMF, IDQ), which has 1 kcps DCR and 10% PDE with hold-off time of 10 $\mu$s [39]. The afterpulse probability distribution of our free-running SPD is further measured using the TDC, as shown in Fig. 5. This distribution is fitted by the following formula [33]

[TABLE]

where $A_{1},\tau_{1},A_{2},\tau_{2},A_{3},\tau_{3},c$ are fitting parameters. Acquiring such distribution is useful for afterpulse correction in the lidar experiment.

5 Afterpulse and count rate corrections for lidar experiment

In order to apply our free-running InGaAs/InP SPD to practical lidar system, we develop a specific algorithm for afterpulse and count rate corrections, after which the lidar performance has been significantly improved. The flow diagram of the correction algorithm is described in Fig. 6. Given a measured afterpulse probability distribution $P_{ap}(i)$ and backscattering signal count rate distribution $R(i)$, the afterpulse count rate $R_{ap}$ in bin $j$ is

[TABLE]

where $A(i,j)$ represents the afterpulse count rate probability in bin $j$ due to a detection click in bin $i$. $A(i,j)$ can be calculated as

[TABLE]

where $P_{nc}(i,j)$ represents the probability that no photon count occurs between bin $i$ and bin $j$, $P_{nap}(i,j)$ represents the probability that no afterpulse occurs between bin $i$ and bin $j$, and $P_{ap}(j;i)$ represents the probability that a photon count in bin $i$ results in an afterpulse count in bin $j$. $P_{nc}(i,j)$ and $P_{nap}(i,j)$ are further calculated by

[TABLE]

[TABLE]

where $bin_{w}$ is the duration time of a bin. Therefore, the count rate in bin $i$ without the afterpulsing contribution is given by

[TABLE]

and further considering the limit of hold-off time the corrected photon count rate in bin $i$ is

[TABLE]

In the field lidar experiment, the measured results of detector count rate and range-corrected signal (normalized $Pr^{2}$) versus range using both SNSPD and InGaAs/InP NFAD are shown in Figs. 7(a) and 7(b), respectively. After correcting the results using InGaAs/InP NFAD by the afterpulse and count rate correction algorithm, one can find out that the lidar performances in terms of normalized $Pr^{2}$ versus range in two cases of using SNSPD and InGaAs/InP NFAD with the corrections agree very well with each other, which clearly shows the usefulness of the correction algorithm. Further, we use the normalized $Pr^{2}$ results using SNSPD as a reference and calculate the relative error of the results using InGaAs/InP NFAD, as shown in Fig. 7(c). In the case without correction, the relative error between the two scenarios increases and then is stable at around 30% with the increase of range. However, with correction such relative error is drastically decreased down to a pretty low level. For instance, in the range between 1 km to 3 km, the maximum error is only 2.2% while the average error is 1.3%. The large errors in the range close to the limit, i.e., $\sim$6 km, are due to the fact that detector count rate significantly decreases in that range and thus the calculation of relative errors is not accurate because of statistical fluctuations. Due to the geometrical overlap factor in the biaxial lidar, the raw signal increases rapidly from zero to maximum in the near range. Since the scaler receives raw signals from two detectors, a slight delay between two channels results in very large errors due to the steep slopes in the very near range, i.e., < 0.5 km. In practice, the results of normalized $Pr^{2}$ in very near range are not used. Such comparison indicates that for lidar applications our free-running InGaAs/InP SPD is completely comparable to SNSPD with the help of afterpulse and count rate corrections.

6 Conclusion

In conclusion, we have presented a fully integrated free-running SPD based on an InGaAs/InP NFAD for 1.5 $\mu$m Mie lidar system. We have optimized the key parameters of InGaAs/InP SPD and investigated the ringing effect of backscattering signal for lidar applications. After afterpulse and count rate corrections, the field lidar performance in terms of normalized $Pr^{2}$ versus range using the free-running InGaAs/InP SPD is considerably comparable to that of using a commercial SNSPD, with only $\sim$2% relative error. Our technique provides a practical solution for single-photon detection dedicated to accurate lidar applications.

Funding

This work has been financially supported by the National Basic Research Program of China Grant No. 2013CB336800, the National Natural Science Foundation of China Grant No. 11674307, and the Chinese Academy of Sciences.

Acknowledgments

We would like to thank Wen-Hao Jiang, Chong Wang and Jiawei Qiu for useful discussions.