# Patent

## Title
Characterization of single-photon detectors using a continuous wave laser source

## Abstract

Operational parameters of a single-photon detector are determined with a continuous wave laser source. At a fixed trigger, a dark count probability and a series of count probabilities at different optical powers are determined. A particular optical power is selected by using a wide-range variable attenuator to attenuate the optical power of the continuous wave laser. The dark count probability and the count probabilities are determined for different trigger rates. The operational parameters include efficiency, afterpulsing constant, and detrap time. The operational parameters are computed by fitting the computed dark count probabilities and count probabilities to a user-defined relationship.

## Background
This application is a continuation of prior U.S. patent application Ser. No. 12/834,956, filed Jul. 13, 2010, which is herein incorporated by reference in its entirety. This application is related to U.S. patent application Ser. No. 12/834,940, filed Jul. 13, 2010, which is herein incorporated by reference in its entirety.


BACKGROUND
The present disclosure relates generally to optical detectors, and more particularly to characterization of single-photon detectors using a continuous wave laser source.
An optical detector that has a sufficiently high detection sensitivity that it can detect a single photon is referred to as a single-photon detector (SPD). Single-photon detectors operating over various ranges of wavelengths have been used for various applications, including ultrahigh resolution optical time domain reflectometry, photon-correlation spectroscopy, and device characterization. At present, there is particularly strong interest in the application of SPDs for optical telecommunications. One method for high-security transmission over an optical telecommunications network is quantum cryptography, which utilizes quantum key distribution via single photons. Fiberoptic transport networks typically operate over a wavelength range of  ̃1.3-1.6 μm. Over this wavelength range, single-photon detectors are typically based on avalanche photodiodes (APDs) fabricated from InP/InGaAs/InGaAsP heterostructures.
A critical performance parameter characterizing an SPD is its detector efficiency. For optical telecommunications, SPDs typically have a low detector efficiency ( ̃10-30%); furthermore, the detector efficiency drifts with time. Current methods for measuring the detector efficiency of SPDs use complicated timing control schemes and ultralow input power levels. Current measurement equipment is typically bulky and expensive, including electronic time delays, pulsed lasers, high dynamic range attenuators, and high sensitivity power meters. What is needed is a simple, low-cost, compact system for measuring the detector efficiency, and other characteristics, of a single-photon detector.
BRIEF SUMMARY
Operational parameters of a single-photon detector are determined with a continuous wave laser source. For each trigger rate selected from a plurality of trigger rates, a dark count probability is computed from recorded dark counts, and a plurality of count probabilities are computed from recorded counts. Each particular count probability is determined from a particular optical power selected from a plurality of optical powers. The operational parameters include efficiency, afterpulsing constant, and detrap time. The operational parameters are computed by fitting the computed dark count probabilities and count probabilities to a user-defined relationship.
These and other advantages of the disclosure will be apparent to those of ordinary skill in the art by reference to the following detailed description and the accompanying drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows a schematic of a generic optical telecommunications system;
FIG. 2A shows a schematic of a measurement system, including a reference pulse laser source and a power meter, for measuring the operational parameters of a single-photon detector;
FIG. 2B shows a schematic of a field measurement system for measuring the operational parameters of a single-photon detector;
FIG. 3 shows bias pulses applied to a single-photon detector;
FIG. 4 shows plots of probability count as a function of the average number of photons per detector gate;
FIG. 5 shows plots of probability count for a continuous wave laser source and a pulsed laser source as a function of the average number of photons per detector gate at a fixed trigger rate;
FIG. 6A shows plots of detector efficiency for a continuous wave laser source as a function of function of the average number of photons per detector gate for different values of trigger rate;
FIG. 6B shows plots of afterpulse parameter as a function of time;
FIG. 7 shows a flowchart of a method for measuring the operational parameters of a single-photon detector; and
FIG. 8 shows a schematic of a computational system for implementing a measurement control system.



DETAILED DESCRIPTION
FIG. 1 shows a schematic of a generic optical telecommunications system. Multiple optical transceivers (XCVRs) send and receive lightwave signals via optical transport network 102 (such as a fiberoptic network). Shown are four representative transceivers, referenced as XCVR 1 104, XCVR 2 106, XCVR 3 108, and XCVR 4 110, respectively. The receivers in the transceivers include single-photon detectors (SPDs). Each transceiver has a corresponding transmit wavelength (λn,T) and a corresponding receive wavelength (λn,R), where n=1-4.
FIG. 2A shows a dual-mode measurement system for measuring the operational parameters of a single-photon detector. The measurement system can compare measurements using a short-pulse laser source with measurements using a continuous wave (CW) laser source. Continuous wave laser source 202 transmits optical beam 231, and short-pulse laser source 204 transmits optical beam 233. Optical beam 231 and optical beam 233 are transmitted to coupler 206 [for example, an equal-intensity (50:50) coupler]. The output of coupler 206 is optical beam 235, which is transmitted through variable attenuator 208. The output of variable attenuator 208 is transmitted as optical beam 237 to single-photon detector 210, or, alternatively, to power meter 212. The electronic charge pulses from the output of single-photon detector 210 are detected and counted by counter 214. The operations of CW laser source 202, short-pulse laser source 204, variable attenuator 208, single-photon detector 210, counter 214, and power meter 212 are controlled by measurement control system 220. Various elements of the measurement system can also be manually controlled by a user.
In the example shown in FIG. 2A, the transport medium is optical fiber, and the wavelength of all the optical beams (optical beam 231—optical beam 237) is λ=1.55 μm. As discussed above, SPDs for wavelengths used in optical telecommunications are typically based on an avalanche photodiode (APD). Operation of an APD is a function of the reverse bias voltage applied across the APD. When the bias voltage exceeds the breakdown voltage, an incident photon can initiate a carrier avalanche, resulting in a charge pulse that can be electronically detected. A carrier avalanche can also be initiated by a source other than an incident photon (such as a trapped charge or a phonon). These other sources can generate background counts, resulting in background noise.
To reduce background noise, an APD can be operated in a gated mode, as shown in FIG. 3. The horizontal axis 302 represents time, and the vertical axis 304 represents bias voltage. Shown are three reference voltages: reference voltage 312 is the breakdown voltage; reference voltage 314, below the breakdown voltage, is the pulse baseline voltage; and reference voltage 316, above the breakdown voltage, is the pulse peak voltage. The bias voltage is modulated by a periodic sequence of pulses (shown as ideal rectangular pulses), referred to as bias pulses (also referred to as gate pulses). Shown are three representative bias pulses: bias pulse 332, bias pulse 334, and bias pulse 336. The amplitude of the bias pulses ranges from the pulse baseline voltage 314 to the pulse peak voltage 316. The pulse width (also referred to as the gate window) is referenced as pulse width 350. The period between pulses is referenced as period 360. The inverse of the period between pulses is referred to as the trigger rate. The interval delimited by the end of one pulse and the start of another pulse is referenced as afterpulse interval 370.
During a bias pulse, the bias voltage is greater than the breakdown voltage, and an incident photon can initiate a carrier avalanche that results in a charge pulse that can be electronically detected. During an afterpulse interval, the bias voltage is less than the bias voltage. The charge state of the APD, however, does not relax instantaneously. During an afterpulse interval, an incident photon can still trigger a carrier avalanche, resulting in a detected signal. In addition, as discussed above, even in the absence of incident photons (dark conditions), other sources can trigger carrier avalanches, resulting in dark counts. Dark counts are dependent on the temperature of the SPD. If the incident photons are generated by a pulsed laser source that is synchronized with the bias pulses, the effects of non-photon induced carrier avalanches can be reduced. For characterization of commercial SPDs, a pulsed laser source at ultralow powers is typically used. A CW laser source at relatively higher powers would be desirable, however, since the measurement system would be simpler and more economical.
Return to FIG. 2A. Short-pulse laser source 204, used for reference measurements, emits optical pulses with an optical pulse width of  ̃30 ps; the short-pulse laser source 204 is synchronized with the single-photon detector 210. Variable attenuator 208 is a double-stage variable attenuator with a total range of 85 dB. Each stage is calibrated independently at relatively high input power levels of −50 dBm by using the power meter 212. The calibration enables precise control of the power at levels below the sensitivity of the power meter 212. By adjusting variable attenuator 208, the average number of photons per detector gate (average number of photons per measurement interval), μ, is adjusted down to 0.02 for both the CW laser source 202 and the short-pulse laser source 204. For each detection event, single-photon detector 210 generates an electrical Nuclear Instrumentation Module Standard (NIM) pulse. The NIM pulses are integrated by electronic counter 214.
To determine the average number of photons per detector gate in the CW regime, the effective gate window of the detector τg is measured. During each counting cycle, a 1-ns-wide electrical bias pulse brings the APD above the breakdown voltage. Conventionally, the maximum efficiency η0 achieved at the peak of the bias pulse is measured with a properly aligned narrow optical pulse and serves as the major specification for commercial SPDs. With a CW laser source, however, the efficiency changes during the duration of the bias pulse in a non-linear fashion. To simplify calculations, an SPD is modelled in the CW regime as operating at its maximum efficiency value η0 for an effective gate window τg. To determine τg experimentally, the probability of a count Pc is measured over a wide power range for both CW and pulsed lasers at a fixed trigger rate R. In the discussion below, R=100 kHz; however, other trigger rates can be used as well. The values of the CW laser power Pcw and the average pulsed laser power Pp are then selected such that the corresponding count probabilities are equal: Pc(Pcw)=Pc(Pp). The effective gate window is then calculated as τg=10(Pp−Pcw−10 log10(R))/10, where the power values are measured in dBm. When averaged over the entire power range, τg is observed to fall in a range of 0.49 ns-0.59 ns for various commercial SPDs.
FIG. 4 shows plots of the measured count probability as a function of the average number of photons per detector gate. The horizontal axis 402 represents the average number of photons per detector gate, μ. The vertical axis 404 represents the measured count probability, Pc. The open squares 406 show the results for CW laser source 202, and the filled dots 408 show the results for short-pulse laser source 204. More specifically, the plots show the measured count probabilities Pc(μ) as a function of the average number of photons per detector gate μcw for CW laser source 202 and as a function of the average number of photons per detector gate μp for short-pulse laser source 204, where




  
    μ
    cw
  
  =
  
    
      
        
          10
          
            (
            
              
                P
                cw
              
              /
              10
            
            )
          
        
        
          1000
          ×
          h
          ×
          v
        
      
      ×
      
        τ
        g
      
      ⁢
      
        
      
      ⁢
      and
      ⁢
      
        
      
      ⁢
      
        μ
        p
      
    
    =
    
      
        
          10
          
            (
            
              
                P
                p
              
              /
              10
            
            )
          
        
        
          1000
          ×
          h
          ×
          v
          ×
          R
        
      
      .
    
  




Here h is Planck's constant and ν=1/λ is the optical frequency. When τg=0.49 ns, the measured count probabilities in the CW and the pulsed case are equal.

The dependence shown in FIG. 4 saturates because avalanche SPDs are not sensitive to the number of registered photons (a registered photon is a photon that sets off an avalanche). For Poissonian photon statistics, this saturation was first described by B. F. Levine et al. (“Near room temperature 1.3 μm single photon counting with a InGaAs avalanche photodiode,”Electron Lett., vol. 20, no. 14, pp. 596-598, July 1984) by the expression:

Pc=1−(1−Pdc)×exp(−μη),  (E1)

where Pdc is the probability of a dark count, μ is the average number of photons per detector gate, and η is the detector efficiency. The measured data, however, deviates from the dependence in (E1) in two ways. First, the measured count probability is found to be somewhat higher at lower μ values; this result suggests the presence of afterpulses, as described by M. A. Itzler et al. (“Single photon avalanche diodes (SPADs) for 1.5 μm photon counting applications,” J. Mod. Opt., vol. 54, no. 2-3, pp. 283-304, February 2007) and X. Jiang et al. (“InGaAsP—InP avalanche photodiodes for single photon detection,” IEEE Journal of Selected Topics in Quantum Electronics, vol. 13, no. 4, pp. 895-905, August 2007). Second, the entire curve seems to be shifted toward higher μ; this result indicates a slight decrease in efficiency. These two effects become more pronounced with higher trigger rates. These deviations due to the afterpulsing effect are accounted for by modifying the expression for the count probability (see below).

The afterpulsing effect is characterized by the conditional probability Q(t) of having an avalanche for the detector gate being open after a time interval t since a successful detection event. It is assumed that this conditional probability decays exponentially with a time constant τ:Q(t)=Q0exp(−t/τ), where Q0 is referred to as the afterpulsing constant, and τ is referred to as the detrap time. Then Qn is the conditional probability of the avalanche occurrence after exactly n clock cycles in the detector that is gated with the trigger rate R:

Qn=Q0exp(−n/Rτ).  (E2)

Within this model, the overall probability of a registered event (count probability) Pc is found to be determined by the relationship:




  
    
      
        
          
            
              
                P
                c
              
              =
              
                
                  P
                  ⁢
                  
                    
                  
                  ⁢
                  1
                  ⁢
                  
                    {
                    •
                    }
                  
                
                +
                
                  
                    1
                    2
                  
                  ⁢
                  
                    
                      
                        P
                        ⁢
                        
                          
                        
                        ⁢
                        2
                        ⁢
                        
                          {
                          •
                          }
                        
                      
                      +
                      
                        P
                        ⁢
                        
                          
                        
                        ⁢
                        3
                        ⁢
                        
                          {
                          •
                          }
                        
                      
                    
                  
                
              
            
            ,
            where
          
          ⁢
          
            

          
          ⁢
          
            
              P
              ⁢
              
                
              
              ⁢
              1
              ⁢
              
                {
                •
                }
              
            
            =
            
              
                2
                -
                
                  e
                  
                    
                      1
                      /
                      R
                    
                    ⁢
                    
                      
                    
                    ⁢
                    τ
                  
                
                -
                
                  
                    (
                    
                      1
                      -
                      
                        P
                        dc
                      
                    
                    )
                  
                  ⁢
                  
                    (
                    
                      1
                      -
                      
                        P
                        ph
                      
                    
                    )
                  
                  ⁢
                  
                    (
                    
                      1
                      -
                      
                        Q
                        0
                      
                    
                    )
                  
                
              
              2
            
          
          ,
          
            

          
          ⁢
          
            
              P
              ⁢
              
                
              
              ⁢
              2
              ⁢
              
                {
                •
                }
              
            
            =
            
              
                [
                
                  2
                  -
                  
                    e
                    
                      
                        1
                        /
                        R
                      
                      ⁢
                      
                        
                      
                      ⁢
                      τ
                    
                  
                  -
                  
                    
                      (
                      
                        1
                        -
                        
                          P
                          dc
                        
                      
                      )
                    
                    ⁢
                    
                      (
                      
                        1
                        -
                        
                          P
                          ph
                        
                      
                      )
                    
                    ⁢
                    
                      (
                      
                        1
                        -
                        
                          Q
                          0
                        
                      
                      )
                    
                  
                
                ]
              
              2
            
          
          ,
          and
        
        ⁢
        
          

        
        ⁢
        
          
            P
            ⁢
            
              
            
            ⁢
            3
            ⁢
            
              {
              •
              }
            
          
          =
          
            4
            ⁢
            
              
                
                  (
                  
                    
                      e
                      
                        
                          1
                          /
                          R
                        
                        ⁢
                        
                          
                        
                        ⁢
                        τ
                      
                    
                    -
                    1
                  
                  )
                
                ⁡
                
                  [
                  
                    1
                    -
                    
                      
                        (
                        
                          1
                          -
                          
                            P
                            dc
                          
                        
                        )
                      
                      ⁢
                      
                        (
                        
                          1
                          -
                          
                            P
                            ph
                          
                        
                        )
                      
                    
                  
                  ]
                
              
              .
            
          
        
      
    
    
      
        (
        E3
        )
      
    
  




Here, Pph is the probability of detecting photons generated by the specific source in the absence of both dark counts and afterpulses. For a Poissonian source of single photons, such as attenuated coherent laser light corresponding to optical beam 237,

Pph=1−exp(−μη).  (E4)

For Q0→0 or τ→0, (E3) reduces to the more simple form (E1).

The derivation of (E3) is accomplished in two steps. The first step consists of evaluating the marginal probability of an afterpulse Pap:

Pap(n)=Pc(n-1)Q1+(1−Pcn-1))Pc(n-2)Q2+(1−Pc(n-1))(1−Pc(n-2))Pc(n-3)Q3+ . . . ,  (E5)

where the superscript n denotes that the corresponding probability is evaluated at the n-th time interval. Note that Qk is a conditional probability of an afterpulse given that no counts occurs during the time t=k/R.

Substitution of (E2) into (E5) yields:

Pap(n)=Pc(n-1)Q1+Pap(n-1)(1−Pc(n-1))Q1/Q0.  (E6)

After a sufficiently long measurement time (sufficiently large number of trigger cycles), the probabilities Pap(n) and Pc(n) reach their stationary values and can be replaced by Pap and Pc, respectively, in (E6), which is then solved for Pap:





  
    
      
        
          P
          ap
        
        =
        
          
            
              
                P
                c
              
              ⁢
              
                Q
                1
              
            
            
              1
              -
              
                
                  (
                  
                    1
                    -
                    
                      P
                      c
                    
                  
                  )
                
                ⁢
                
                  
                    Q
                    1
                  
                  /
                  
                    Q
                    0
                  
                
              
            
          
          .
        
      
    
    
      
        (
        E7
        )
      
    
  




The second step of the derivation consists of extracting Pc from Pc=1−(1−Pph)(1−Pap), where Pap is given by (E7). For Q0→0 or τ→0 (that is, Q1→0), (E7) yields Pap=0. Therefore, Pc, as evaluated from Pc=1−(1−Pdc)(1−Pph), reduces to the simple form given in (E1).
FIG. 5 shows plots of the measured count probability for a CW laser as a function of average number of photons per detector gate at different trigger rates, over a range of approximately 100 kHz to 1 MHz. A different range can be used, depending on the detector. To show sufficient detail, two sets of plots are shown: set 506A in the principal field, and set 506B in the inset. The horizontal axes, horizontal axis 502A and horizontal axis 502B, represent the average number of photons per detector gate,




  
    
      μ
      cw
    
    =
    
      
        
          10
          
            (
            
              
                P
                cw
              
              /
              10
            
            )
          
        
        
          1000
          ×
          h
          ×
          v
        
      
      ×
      
        τ
        g
      
    
  
  ,




over the range from 0.1 to 20. The data is taken at four different trigger rates (R): data points 520, data points 522, data points 524, and data points 526 correspond to measured values for R=100 kHz, 250 kHz, 500 kHz, and 1 MHz, respectively. As discussed below, the plotted lines correspond to fitted curves.

As discussed above, for each trigger rate R, the dependence of Pc on μcw can be described by (E3), computed with the corresponding values of R. The following operational parameters of the detector can then be extracted by fitting (E3) to the experimental data: the efficiency for each trigger rate η(R), the afterpulsing constant Q0, and the detrap time τ.
The rate dependent dark count probability Pdc is first measured. Then a simultaneous fit is performed for all four experimentally obtained dependencies of the count probability on the average photon number per detector gate. By inserting the resulting parameters back into (E3), fitted plots can be made. Two full-range representative fitted plots (corresponding to the low end and the high end of the range of trigger rates) are shown in set 506A of FIG. 5: fitted plot 530 corresponding to R=100 kHz and fitted plot 536 corresponding to R=1 MHz. The high quality of the fits can be clearly seen. By using only a limited range of μcw=3-20 (corresponding to an easily measurable interval of −61 dBm≦Pcw≦−53 dBm) a nearly perfect match is obtained for the data within a rms deviation of 3×10−4. Therefore, the operational parameters of the detector can be extracted from measured photon counts at a relatively high number of incident photons.
FIG. 6A shows plots of the extracted detector efficiency as a function of trigger rate. The horizontal axis 602 represents the trigger rate R (kHz). The vertical axis 604 represents the extracted detector efficiency η(%). Points 610 (open squares) show the values for a CW laser source at relatively high powers; points 612 (filled dots) show the values for a pulsed laser source at relatively high powers. For low trigger rates R=100 kHz and R=250 kHz, the CW laser source and the pulsed laser source yield nearly identical results. At R=250 kHz, the measured efficiency is 21%, which is in close agreement with the manufacturer's specification of 20%. The data indicates that the efficiency drops with trigger rate. The drop is more pronounced for a CW laser source, resulting in a 2% deviation from the value of 19.7% measured with a pulsed laser source at R=1 MHz.
The observed trends in the efficiency are also reflected in the afterpulsing probabilities, as shown in the plots in FIG. 6B. The horizontal axis 632 represents time (μs). The vertical axis 634 represents Q(t)=Q0exp(−t/τ), with the afterpulsing constant Q0 and the detrap time τ extracted from the fits. Plot 640 shows the results for a CW laser source; plot 642 shows the results for a pulsed laser source. The CW laser source and the pulsed laser source yield similar results for times greater than 4-5 μs, that is, corresponding to trigger rates less than 200-250 kHz. For small time values, the afterpulsing effect for a CW laser source grows significantly. The CW detrap time τcw≈2.5-3.0 μs appears to be slightly smaller than the values obtained for the pulsed laser source, τpulsed≈4.0-5.0 μs, which are more consistent with the values reported by the device manufacturers. One possible reason for this discrepancy is that the relaxation in the CW regime might be a combination of several processes beyond afterpulsing; the other processes, for example, might be the charge persistence or twilight effect as being dominant for small time values for a CW laser source (see J. Zhang et al., “Comprehensive characterization of InGaAs/InP avalanche photodiodes at 1550 nm with an active quenching ASIC,” IEEE J. Quantum Electron., vol. 45, Issue 7, pp. 792-799, July 2009, and S. V. Polyakov and A. L. Migdall, “High accuracy verification of a correlated photon-based method for determining Photon counting detection efficiency”, Optics Express, vol. 15, No. 4, pp. 1390-1407, February 2007).
An embodiment of a measurement system for characterizing the operational parameters of a single-photon detector is shown in FIG. 2B. It is similar to the measurement system in FIG. 2A, except the components used in FIG. 2A for initial calibration of the variable attenuator, comparison with a reference pulsed laser source, and determination of the effective gate window, are no longer present. Note that the effective gate window is a parameter that can also be reported by the detector manufacturer. In one embodiment, the SPD under test is decoupled from its normal operating network connection (see FIG. 1) and coupled to the measurement system. As described above with reference to FIG. 2A, continuous wave laser source 202 transmits optical beam 231 through variable attenuator 208. The output of variable attenuator 208, optical beam 237, is transmitted to single-photon detector 210. The electronic charge pulses from the output of single-photon detector 210 are detected and counted by counter 214. Measurement control system 220 controls the operation of CW laser source 202, variable attenuator 208, single-photon detector 210, and counter 214.
FIG. 7 shows a flowchart of steps for a method for measuring the operational parameters of single-photon detector 210. The SPD measurement system is configured as shown in FIG. 2B. For a specific value of R, counts are measured under dark conditions and under multiple values of input CW laser power. The measurements are then repeated for multiple values of R. After all the measurements have been completed, the operational parameters are computed. Since the operational parameters are a function of temperature, the temperature of single-photon detector 210 is stabilized at a user-specified temperature prior to the set of measurements. Measurements can be taken over a user-specified set of temperatures, if desired.
In step 702, the trigger rate R is set at a first value for single-photon detector 210. The process then passes to step 704, in which the optical power from CW laser source 202 is turned OFF. The process then passes to step 706, in which dark counts Ddark over a specified measurement time are recorded by counter 214. The dark count probability is then computed from Pdc=Ddark/R. The process then passes to step 708, in which the optical power from CW laser source 202 is turned ON. The process then passes to step 710, in which variable attenuator 208 is adjusted to set the input optical power to single-photon detector 210 to a first value. The process then passes to step 712, in which the counts D are recorded over a specified measurement time. The count probability is then computed from Pc=D/R.
The process then passes to decision step 714. If more measurements at the same trigger rate are to be taken, then the process returns to step 710, in which variable attenuator 208 is adjusted to set the input optical power to single-photon detector 210 to a different value. For a fixed value of trigger rate, step 710-step 714 are iterated for a user-specified number of cycles. The process then passes to decision step 716. If measurements are to be taken at a different trigger rate, then the process returns to step 702, in which the trigger rate is set to a different value. Step 702-step 716 are iterated for a user-specified number of cycles. When the set of measurements has been completed, then the process passes from step 716 to step 718.
In step 718, the following operational parameters of the detector are computed by fitting the measured data to (E3): the efficiency as a function of trigger rate η(R), the afterpulsing constant Q0, and the detrap time τ.
FIG. 8 shows a schematic of an embodiment of a computational system for implementing a measurement control system 220 (see FIG. 2B). One skilled in the art can construct the computational system 802 from various combinations of hardware, firmware, and software. One skilled in the art can construct the computational system 802 from various combinations of electronic components, such as general purpose microprocessors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), and field-programmable gate arrays (FPGAs).
Computational system 802 comprises computer 804, which includes a central processing unit (CPU) 806, memory 808, and data storage device 810. Data storage device 810 comprises at least one non-transitory, persistent, tangible computer readable medium, such as non-volatile semiconductor memory, a magnetic hard drive, and a compact disc read only memory.
Computational system 802 can further comprise one or more user input/output interfaces that interface computer 804 with user input/output devices. For example, user input/output interface 820 interfaces computer 804 with user input/output device 840. Examples of user input/output device 840 include a keyboard, a mouse, and a local access terminal. Data, including computer executable code, can be transferred to and from computer 804 via user input/output interface 820.
Computational system 802 can further comprise video display interface 822, which drives video display 842 and can receive user input from video display 842.
Computational system 802 can further comprise one or more instrumentation interfaces. For example, instrumentation interface A 824 interfaces computer 804 with CW laser source 202; instrumentation interface B 826 interfaces computer 804 with variable attenuator 208; instrumentation interface C 828 interfaces computer 804 with single-photon detector 210; and instrumentation interface D 830 interfaces computer 804 with counter 214 (see FIG. 2B).
Computational system 802 can further comprise one or more communications network interfaces that interface computer 804 with communications networks, such as local area networks and wide area networks. Data, including computer executable code, can be transferred to and from computer 804 via communications network interfaces. For example, communications network interface 832 interfaces computer 804 with communications network 852. Computational system 802 can be remotely accessed and controlled via communications network 852.
As is well known, a computer operates under control of computer software, which defines the overall operation of the computer and applications. CPU 806 controls the overall operation of the computer and applications by executing computer program instructions that define the overall operation and applications. The computer program instructions can be stored in data storage device 810 and loaded into memory 808 when execution of the program instructions is desired. The method steps shown in the flowchart in FIG. 7 can be defined by computer program instructions stored in memory 808 or in data storage device 810 (or in a combination of memory 808 and data storage device 810) and controlled by the CPU 806 executing the computer program instructions. For example, the computer program instructions can be implemented as computer executable code programmed by one skilled in the art to perform algorithms implementing the method steps shown in the flowchart in FIG. 7. Accordingly, by executing the computer program instructions, the CPU 806 executes algorithms implementing the method steps shown in the flowchart in FIG. 7.
The foregoing Detailed Description is to be understood as being in every respect illustrative and exemplary, but not restrictive, and the scope of the inventive concept disclosed herein is not to be determined from the Detailed Description, but rather from the claims as interpreted according to the full breadth permitted by the patent laws. It is to be understood that the embodiments shown and described herein are only illustrative of the principles of the present disclosure and that various modifications may be implemented by those skilled in the art without departing from the scope and spirit of the disclosure. Those skilled in the art could implement various other feature combinations without departing from the scope and spirit of the disclosure.

## Claims

1. A method comprising:
recording dark counts with a single-photon detector at a trigger rate;
computing a dark count probability for the trigger rate based on the dark counts;
recording counts at the trigger rate with the single-photon detector, the counts generated by a continuous wave laser at an optical power;
computing a count probability based on the counts; and
calculating an operating parameter of the single-photon detector by fitting the dark count probability and the count probability to a user-defined relationship defining a probability of detecting photons generated by a specific source in an absence of the dark counts and afterpulses, wherein the probability of detecting the photons generated by the specific source in the absence of the dark counts is exponentially related to an average number of the photons per detector gate.

2. The method of claim 1, wherein the calculating the operating parameter comprises calculating one of:
an efficiency η;
an afterpulsing constant Q0; and
a detrap time τ.

3. The method of claim 1, further comprising:
generating the optical power by attenuating a constant optical power generated by the continuous wave laser.

4. The method of claim 1, wherein a wavelength of the continuous wave laser ranges from 1.3 microns to 1.6 microns.

5. The method of claim 1, wherein the count probability is derived from evaluating a probability of an afterpulse occurring at a predetermined time interval after a photon count.

6. The method of claim 1, wherein the trigger rate is fixed.

7. The method of claim 3, wherein the user-defined relationship comprises:



  
    
      
        P
        c
      
      =
      
        
          P
          ⁢
          
            
          
          ⁢
          1
          ⁢
          
            {
            •
            }
          
        
        +
        
          
            1
            2
          
          ⁢
          
            
              
                P
                ⁢
                
                  
                
                ⁢
                2
                ⁢
                
                  {
                  •
                  }
                
              
              +
              
                P
                ⁢
                
                  
                
                ⁢
                3
                ⁢
                
                  {
                  •
                  }
                
              
            
          
        
      
    
    ,
    
      wherein
      ⁢
      
        :
      
    
  




  
    
      
        P
        ⁢
        
          
        
        ⁢
        1
        ⁢
        
          {
          •
          }
        
      
      =
      
        
          2
          -
          
            e
            
              
                1
                /
                R
              
              ⁢
              
                
              
              ⁢
              τ
            
          
          -
          
            
              (
              
                1
                -
                
                  P
                  dc
                
              
              )
            
            ⁢
            
              (
              
                1
                -
                
                  P
                  ph
                
              
              )
            
            ⁢
            
              (
              
                1
                -
                
                  Q
                  0
                
              
              )
            
          
        
        2
      
    
    ,
    
      

    
    ⁢
    
      
        P
        ⁢
        
          
        
        ⁢
        2
        ⁢
        
          {
          •
          }
        
      
      =
      
        
          [
          
            2
            -
            
              e
              
                
                  1
                  /
                  R
                
                ⁢
                
                  
                
                ⁢
                τ
              
            
            -
            
              
                (
                
                  1
                  -
                  
                    P
                    dc
                  
                
                )
              
              ⁢
              
                (
                
                  1
                  -
                  
                    P
                    ph
                  
                
                )
              
              ⁢
              
                (
                
                  1
                  -
                  
                    Q
                    0
                  
                
                )
              
            
          
          ]
        
        2
      
    
    ,
    
      

    
    ⁢
    
      
        P
        ⁢
        
          
        
        ⁢
        3
        ⁢
        
          {
          •
          }
        
      
      =
      
        4
        ⁢
        
          
            (
            
              
                e
                
                  
                    1
                    /
                    R
                  
                  ⁢
                  
                    
                  
                  ⁢
                  τ
                
              
              -
              1
            
            )
          
          ⁡
          
            [
            
              1
              -
              
                
                  (
                  
                    1
                    -
                    
                      P
                      dc
                    
                  
                  )
                
                ⁢
                
                  (
                  
                    1
                    -
                    
                      P
                      ph
                    
                  
                  )
                
              
            
            ]
          
        
      
    
    ,
    
      

    
    ⁢
    
      
        P
        ph
      
      =
      
        1
        -
        
          exp
          ⁡
          
            (
            
              -
              μη
            
            )
          
        
      
    
    ,
  


Pph=1−exp(−μη),
Pc is the count probability,
Pdc is the dark count probability,
μ is the average number of the photons per detector gate, and
R is the trigger rate.

8. The method of claim 3, wherein the optical power is attenuated using a wide-range variable attenuator.