Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Spaceborne, low-noise, single-photon detection for satellite-based quantum communications

Open Access Open Access

Abstract

Single-photon detectors (SPDs) play important roles in highly sensitive detection applications, such as fluorescence spectroscopy, remote sensing and ranging, deep space optical communications, elementary particle detection, and quantum communications. However, the adverse conditions in space, such as the increased radiation flux and thermal vacuum, severely limit their noise performances, reliability, and lifetime. Herein, we present the example of spaceborne, low-noise, high reliability SPDs, based on commercial off-the-shelf (COTS) silicon avalanche photodiodes (APD). Based on the high noise-radiation sensitivity of silicon APD, we have developed special shielding structures, multistage cooling technologies, and configurable driver electronics that significantly improved the COTS APD reliability and mitigated the SPD noise-radiation sensitivity. This led to a reduction of the expected in-orbit radiation-induced dark count rate (DCR) increment rate from ∼219 counts per second (cps) per day to ∼0.76 cps/day. During a continuous period of continuous operations in orbit which spanned of 1029 days, the SPD DCR was maintained below 1000 cps, i.e., the actual in-orbit radiation-induced DCR increment rate was ∼0.54 cps/day, i.e., two orders of magnitude lower than those evoked by previous technologies, while its photon detection efficiency was > 45%. Our spaceborne, low-noise SPDs established a feasible satellite-based up-link quantum communication that was validated on the quantum experiment science satellite platform. Moreover, our SPDs open new windows of opportunities for space research and applications in deep-space optical communications, single-photon laser ranging, as well as for testing the fundamental principles of physics in space.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Single-photon detectors (SPDs) have been extensively used as one of the key components in extremely weak light detectionn applications over the past two decades. These include light detection and ranging [13], laser time transfer [4,5], scintillation detection of elementary particles [6], deep-space optical communications [7,8], and quantum communications [917]. Among these applications, quantum communications constitute the most rigorous application that demands increased SPD detection efficiency (DE) and low-dark count rate (DCR), given that the optical signal cannot be amplified like classical light signals according to the quantum noncloning principle [18] to guarantee the high fidelity of the single photon during its transmission. In addition to the establishment of a feasible large-scale quantum communication, increased DE, and low DCR, SPDs can effectively improve and extend existing application performances. For example, use of the SPD to replace the Lidar’s light power detector can significantly reduce the dependence of high-power lasers and optical amplifiers [19,20], and can efficiently simplify the system design and reduce power consumption.

At present, the available SPDs mainly include photomultiplier tubes (PMTs), avalanche photodiodes (APDs), and superconducting nanowire single-photon detectors (SNSPDs). In satellite-based quantum communication which is the primary interest of this work, an extensively used wavelength is near 800 nm. In this wavelength, PMTs suffer from low-detection efficiency (below 10%) [21,22], and SNSPDs require heavy and bulky cryocoolers. By contrast, low-cost and compact silicon-based APDs have an ∼65% detection efficiency near the 800 nm wavelength, and are more suited for our space-based missions [23]. Given the use of state-of-the-art, silicon-based APDs, the Execelitas’s commercial off-the-shelf (COTS) super-low K factor (SLiK) silicon APD can achieve a higher DE with a special electric field distribution [24] compared with other silicon APD structures, such as reach-through structures. However, this also led to higher radiation-noise sensitivity [25]. However, given the adverse conditions in space, the silicon-APD-based SPD performances, especially the noise performance and reliability would face unprecedented challenges, as evidenced in previous pioneering work and space applications, such as in the National Aeronautics and Space Administration’s (NASA’s) ice, cloud, and land elevation satellite (ICESat) [13,26] and in Chinese navigation satellites [5]. For example, in NASA’s ICESat, eight lab-assembled, silicon-based APD SPDs were intended for cloud and aerosol backscattering measurements in the geoscience laser altimeter system (GLAS). Even though they were carefully designed and implemented, three of them failed during the GLAS system level tests prior to the ICESat launch [27]. Owing to space radiation, their radiation-induced DCR increased ∼60 cps per day per SPD during orbit [3], which is unacceptable for quantum communications.

The quantum experiment science satellite (QUESS), i.e., Micius, is the satellite designed to carry out quantum experiments [28,29] at the space scale. One of its targets is to conduct ground-satellite quantum teleportation [16] that constitutes the primary interest of our work. The ground–satellite quantum teleportation is the first application of the spaceborne, low-noise SPD, receives uplink configurations, and places the complicated multiphoton setup on the ground. Compared to downlink configurations, two substantial challenges for uplink channels are the atmospheric turbulence at the beginning of the transmission path, which would cause a) increased link attenuation owing to beam wandering and broadening, and b) in-orbit receiver SPD radiation-induced DCR increases, thereby demanding SPDs with increased DEs and low DCR.

To evaluate the SLiK silicon APD sensitivity to radiation damage and in-orbit performance, we have exposed the APDs in the proton flux produced in the laboratory. According to the experimental results, the in-orbit radiation-induced DCR increment rate of the SLiK silicon APD would be ∼219 cps/day. To mitigate the APD radiation-induced DCR increment rate and guarantee the SPD reliability, we have developed special radiation shield structures, multistage cooling technologies, and special driver electronics for the spaceborne, low-noise SPD. Based on these measures, the expected SPD in-orbit DCR increment rate would be ∼0.76 cps/day, thus satisfying the satellite-based quantum communication applications.

Before launch, a series of space environmental tests were carried out to verify the SPD performances and reliabilities, and the results were proven satisfactory. The SPDs were then launched into orbit onboard the QUESS on 16 August 2016. During the past two years of in-orbit operations, we observed an ultralow DCR increase of ∼0.54 cps/day by cooling the APDs to ∼ −50°C, and by adjusting their excess voltages. This DCR increase was two orders of magnitude lower than previous results [3]. As the example of spaceborne low-noise SPDs, they were used in the demonstration of a ground–satellite quantum teleportation experiment [16], during which the DCRs of the SPDs remained below 150 cps. Furthermore, we evaluated numerically the importance of spaceborne, low-noise, Si APD SPDs in ground-to-satellite quantum key distributions (QKDs) [3035]. The results showed that our SPDs could support ground-to-satellite QKDs for 4.6 years compared to conventional detectors which can only last a month.

2. Challenges and solutions

2.1 Radiation damage and simulations

Space radiation comprises mainly solar energetic particle events, galactic cosmic radiation, and trapped radiation belts. For a low Earth orbit (LEO), the major radiation sources are trapped radiation belts which consist of trapped electrons and protons in the South Atlantic Anomaly, and solar protons at the poles of the Earth’s magnetosphere [3638]. Trapped electrons can be shielded effectively to negligible levels with aluminum with a thickness of ∼3 mm [25]. However, the radiation of protons has two major effects on the silicon lattice in silicon-based APDs, namely, ionizing dose damage and accumulated displacement damage [39]. For silicon APDs, the accumulated displacement damage is the major factor that increases the bulk leakage current appreciably or the DCR when the APD is operated in the Geiger mode.

The displacement damage creates different types of silicon lattice defects in the depletion region of APD, whereby the primary initial defects include vacancies, and interstitials. Accordingly, complex defects, or defect centers, are then formed with various numbers of vacancies and interstitials, such as Frenkel pairs, divacancies, vacancy–phosphorus pairs, and defect clusters [40]. These defect centers constitute new electron–hole pair generation centers, influence the properties of existing generation centers [41], and subsequently increase the DCR of the APD. The concentration of the defects depends only on the radiation particle nonionizing energy loss (NIEL, total energy that goes into displacements) and not on the type and initial energy of the particle [42,43]. Based on this fact, it is practical to use mono-energetic particles, such as 50 MeV protons, to scale and simulate the displacement damage of different particles. Therefore, we used 50 MeV protons, which can penetrate the front glass window (the window in front of the APD was 1.0 mm thick and was made of Borosilicate glass which had a density of $\textrm{2}\textrm{.6 g/c}{\textrm{m}^3}$) of APD package with slight spectra spreading to estimate and simulate the displacement damage to the SLiK silicon APD produced by the radiation particles in the in-orbit environment. In the case of the QUESS, its orbit features a circular Sun-synchronous orbit of approximately 500 km and an inclination of 97.4°. Based on the space environment information system (SPENVIS) [44], we used the AP8-min and AE8-min models to evaluate the total displacement damage dose, which was equal to ${\sim} 5.92 \times {10^6}{\textrm{MeV/g(Si)}}$ (behind the 3 mm thick aluminum of the satellite case) for durations of two years, or ${\sim} 8.10 \times {10^3}{\textrm{MeV/g(Si)}}$ per day. In the case of the APD operating in Geiger mode, its DCR can be expressed as [45],

$$\textrm{DCR} = \textrm{S} \cdot \int_{{\textrm{W}_1}}^{{\textrm{W}_1} + {\textrm{W}_2}} {{\textrm{P}_{\textrm{pair}}}(\textrm{x}) \cdot {\textrm{G}_{\textrm{TOT}}}\textrm{dx}}$$
where ${\textrm{W}_1}$ and ${\textrm{W}_2}$ are the position and length of the depletion in the APD, $\textrm{P}\textrm{pair}$ is the total probability of induction of an avalanche breakdown for an electron or a hole, and ${\textrm{G}_{\textrm{TOT}}}$ is the total carrier generation rate. According to the Shockley–Read–Hall theory, and following the counting of all carrier generation centers induced by the material’s natural or radiation-induced defects, the total carrier generation rate, $\textrm{G}_{\textrm{TOT}}$, can be expressed as,
$${\textrm{G}_{TOT}} = \sum\limits_{d = 1}^\textrm{D} {\frac{{{n_i}}}{{{\tau _n}\exp \left[ {\frac{{ - ({E_{T, d}} - {E_i})}}{{KT}}} \right] + {\tau _p}\exp \left[ {\frac{{ - ({E_{T, d}} - {E_i})}}{{KT}}} \right]}}}$$
where $\textrm{D}$ is the number of defects in the depletion area, which increases almost linearly with the displacement damage dose or the in-orbit time, ${E_{T, d}}$ is the $d\textrm{t}h$ carrier generation center energy level, ${E_i}$ is the intrinsic Fermi energy level, and ${n_i}$ is the intrinsic carrier concentration. Additionally, ${\tau _n}$ and $\tau p$ are the lifetimes of the electron and hole, respectively, K is the Boltzmann constant, and T is the absolute temperature.

As the radiation dose or in-orbit time increase, the number of defects in the APD depletion area increases. Correspondingly, the APD DCR increases linearly as a function of the displacement damage dose or the in-orbit time. However, owing to the different electron–hole pair generation rate for different carrier generation centers, the variations of the DCR increment rate as a function of the displacement damage dose or the in-orbit time would temporarily deviate from linearity at small dose spans or short in-orbit time spans (e.g., days or months). However, it would be a very difficult task to derive the APD macroscopic properties, such as the noise performance, from the characteristics of the individual microscopic carrier generation center, which needs to know the carrier generation center energy level in the bandgap, the carrier generation center concentrations associated with each level, the capture and emission probabilities for electrons and holes for each level, etc. In the practical analysis and calculation of the APD noise performance, a more direct and easier mode is used in the APD radiation damage [4648], namely the DCR increment $\Delta \textrm{DCR}$. Owing to radiation damages, this variable can be expressed as,

$$\frac{{\Delta DCR}}{V} = \frac{{{n_i}\phi }}{{2{K_{gn}}}}$$
where $\textrm{V}$ is the APD depletion region volume, ${n_i}$ is the intrinsic carrier density, $\phi$ is the radiation dose, and $K\textrm{gn}$ is the damage coefficient for the material type in the depletion region. In the case of the SLiK APD, $\textrm{V}$, ${n_i}$, and $K\textrm{gn}$, are constants. The DCR increment is linearly proportional to the displacement damage dose.

To evaluate the APD performance degradation in the presence of a radiation environment, we have performed a series of radiation tests on the ground to evaluate how the displacement damage affects the macroscopic properties of the APDs. Seven samples based on Excelitas SLiK APDs and labeled #1–7 were radiated with a proton flux of 50 MeV. Figure 1 shows the DCRs during the radiation test, which increase linearly as a function of radiation fluence in agreement with previous findings [1]. For an operating temperature of −12°C and for an excess voltage of 17 V, the DCR sensitivity is ${\sim} 0.027{\textrm{cps}} \cdot {\textrm{g(Si)/MeV}}$. Combined with the estimated in-orbit radiation dose, the DCR increase is ∼219 cps/day, which is too high for several space applications, including quantum communications [912,1417]. In addition to DCRs, we also measured the APD breakdown voltage, responsivity, and dark current. No significant changes were observed for the responsivity and dark current, whereas the breakdown voltage presented a measurable increment up to ∼3 V when the displacement damage dose reached $3.70 \times {10^9}{\textrm{MeV/g(Si)}}$. However, this fluence value is far beyond the total displacement damage dose during two years of QUESS operations.

 figure: Fig. 1.

Fig. 1. Dark count rate (DCR) versus displacement damage dose. During radiation, sample #2 was cooled to −35°C and was biased with an excess voltage of 17 V, while samples #3 and #7 were unbiased at 25°C, and the other samples were cooled to −12°C and biased with an excess voltage of 17 V. For DCR measurements, all samples were cooled to −12°C and biased with an excess voltage of 17 V.

Download Full Size | PDF

2.2 Shield and thermal design

Plainly, an effective direct way to reduce the radiation-induced APD DCR increment rate is to reduce the radiation dose to APD, such as the selection of low orbits, small-orbit inclination angles, and the addition of protection shield layers. For a given experimental mission, the satellite orbit parameters are basically unmodifiable. Theoretically, radiation can be fully shielded with a shield layer with a sufficient thickness, but it needs extremely thick shields owing to secondary particle effects [49,50]. Based on the use of the SPENVIS tools, we had evaluated the shielding effects of the shield layers with different aluminum thicknesses, as shown as Fig. 2. As the shield thickness increases, the secondary particles dominate in the APD equivalent radiation dose, and the equivalent 50.0 MeV proton fluence remains almost unchanged after the AI shield thickness reaches 200 mm. According to the previous estimation, the displacement damage dose should be lower than $2.70 \times {10^4}{\textrm{MeV/g(Si)}}$, while the in-orbit radiation of the SPD induced a DCR increment which was less than 1 cps/day. This required an AI shield thickness > 400 mm. Correspondingly, the quality of the shielding layer was ∼723 kg, which is unaffordable for most satellites, especially for small satellites.

 figure: Fig. 2.

Fig. 2. Displacement damage dose versus spherical Al shield radius.

Download Full Size | PDF

Taking into account the weight and volume limitation of the satellite, we have traded off the shield performance and the volume, weight, and thermal capacity. Furthermore, an aluminum shielding layer with a thickness of 12 mm—whose specific heat capacity was ∼4 times of tantalum—and a tantalum shielding layer with a thickness of 4 mm were added around the APD, the total quality was ∼5 kg, and the equivalent aluminum thickness was ∼22 mm. This reduced the radiation dose damage by ∼2.5 times. In addition, to achieve the APD deep cooling, the radiation shield layer was thermally isolated and cooled to approximately −15°C by radiation cooling. The shielding and installation of the spaceborne low-noise Si APD single photon detectors, namely, SPD #1 and SPD #2, as shown in Fig. 3.

 figure: Fig. 3.

Fig. 3. Shielding and installation of the spaceborne, low-noise, Si APD, single photon detectors.

Download Full Size | PDF

According to Eq. (2), the carrier generation rate of the carrier generation center decreased exponentially as a function of temperature. Consequently, the APD DCR can be effectively reduced by cooling the APD to lower temperatures, which is consistent with previous findings [51,52]. Furthermore, we performed a series of tests on our APDs at different operating temperatures, and excess voltages to validate the effectiveness of cooling to mitigate the radiation-induced DCR.

Figure 4 shows the test results for one irradiated sample, illustrating that the DCR increases exponentially as a function of temperature and linearly as a function of the excess voltage. In addition, the slope S for DCR versus the excess voltage also increases exponentially as a function of the operating temperature T. For sample #4, the slope can be expressed as,

$$S = \alpha \exp (\beta T)$$
We obtained $\alpha = 6807.53$ and $\beta ={-} 0.09757$ from the testing data. According to Eq. (4), at the same excess voltage, the DCR decreased by ∼2.65 times for a temperature drop of 10°C. This demonstrated the effectiveness of deep cooling to mitigate the radiation damage. By further cooling the APD operating temperature from −12°C to −60°C, the DCR decreased by ∼108 times. Combination of the shield and deep cooling yielded the expected in-orbit DCR increment of ∼0.76 cps/day.

 figure: Fig. 4.

Fig. 4. Deep cooling effect on DCR of sample #4. The results for the other six samples are the same. S represents the slope for the DCR versus excess voltage plots at different temperatures.

Download Full Size | PDF

2.3 Driving electronics and tests

In our implementation, each SPD contained an Excelitas SLiK APD in a TO-can, with an active diameter of 180 µm and ∼60% photon detection efficiency at a wavelength of 800 nm, a two-stage thermal electrical cooler (TEC), and a thermistor. The two-stage TEC and the thermistor were used to realize closed-loop control of the temperature. The package was gas-sealed to reduce the damage from other external factors, such as water vapor and extreme temperature. The driving circuit of each SPD, as shown in Fig. 5, included a robust passive quenching circuit, a high-voltage (HV) manager, an operating temperature controller, and a signal discrimination circuit.

 figure: Fig. 5.

Fig. 5. Schematic of the driving circuit for single-photon detectors (SPDs) (ADC, analog-to-digital-converter; DAC, digital-to-analog-converter; COMP, comparator; DFF, D Flip-flop; PWM, pulse-width modulation; TEC, thermal electric cooler, HV, high voltage. The ADCs, DACs, and PWMs, are connected to a field-programmable gate array (FPGA) in which the temperature and the high-voltage (HV) closed loop control were carried out.

Download Full Size | PDF

According to Eq. (1), the APD DCR was decided by the carrier generation rate and the total probability, $\textrm{P}\textrm{pair}$. The avalanche breakdown increased linearly as a function of the excess voltage (the voltage above the APD breakdown voltage) because the electric field in the APD depletion region was basically uniform based on SLiK structure, as shown in Fig. 4. The APD itself DCR can be expressed as,

$${d_s} = \textrm{SVe}$$
where $\textrm{S}$ is the slope of the plot of the DCR versus the excess voltage which depends on the operating temperature and radiation dose, and $Ve$ is the excess voltage. Conversely, the APD detection efficiency was decided by the quantum efficiency and the photon-induced carrier breakdown probability. The APD detection efficiency is a constant after the APD was produced. By contrast, the photon-induced carrier breakdown probability increased exponentially with the excess voltage [24]. Thus, the APD detection efficiency can be approximately expressed as,
$${\eta _d} = \eta (1 - \exp ( - Ve/Vc))$$
in which $\eta$ is the APD quantum efficiency and equals ∼0.8, $V\textrm{c}$ is the APD characteristic voltage which depends on the depletion layer thickness and the weighted average of the ratio of the ionization coefficient of holes to that of electrons. In the case of the SLiK APD, $V\textrm{c}$ was ∼9.1 V, according to our test results. Based on Eqs. (5) and (6), we can further optimize the system performance through the adjustment of the excess APD voltage. For example, in the uplink ground–satellite experiment [16], the teleportation fidelity can be expressed as,
$$F = \frac{{{\eta _\textrm{l}}{\eta _d}{F_s}{F_l} + 1.795{R_d}{t_w}}}{{{\eta _\textrm{l}}{\eta _d} + 3.59{R_d}{t_w}}}$$
where ${\eta _\textrm{l}}$ is the is the link efficiency from the ground photon source to the input of APD, and is approximately equal to −45 dB, and ${\eta _d}$ is the APD detection efficiency. ${\textrm{F}_s}$ and ${F_l}$ are the photon source and link fidelities, which are equal to ∼0.8 and ∼0.99, respectively. ${R_\textrm{d}}$ is the APD dark count rate, including the link background photon count rate which is ∼150 cps, the thermal noise dark count and the radiation-induced dark count. Additionally, ${t_w}$ is the time window, which is ∼1 ns. According to Eqs. (5) and (7), we can obtain the teleportation fidelity at different excess voltages and radiation damage levels, as shown in Fig. 6.

 figure: Fig. 6.

Fig. 6. Teleportation fidelity at different excess voltages and excess DCR voltages slopes.

Download Full Size | PDF

Lowering the APD excess voltage can improve the teleportation fidelity, especially when the excess DCR voltage slope S is very large owing to radiation damage. Of course, lowering the APD excess voltage would reduce the APD detection efficiency. Accordingly, it will take more time to obtain an adequate number of effective detection events. Based on this fact, we have designed an adjusted HV manager with a broad adjustable range. In the HV manager, a coarse voltage of ∼300 V is generated using voltage regulator diodes. The high voltage is further tuned through a digital-to-analog converter (DAC) for APD at different operating temperatures and detection efficiencies, as shown in Fig. 5. In addition, the temperature drift of the voltage regulator diodes was compensated with the use of a closed-loop controller.

For deep-cooling APD, multistage cooling technologies were introduced, including a stage radiation cooling and a two-stage TEC. The former was implemented with a heat pipe, the shield layer was cooled to ∼15°C, and the APD was installed on the shield layer. To achieve additional cooling, the APD and fine control of the APD operating temperature, a digital proportional, integral, derivative (PID) closed-loop controller based on a two-stage TEC and a thermistor, were implemented. The value of the operating temperature was adjusted according the APD’s DCR in a radiation environment to improve SPD reliability and save satellite energy.

In addition, thermal annealing also partially reduced the radiation-induced DCR [25,53,54]. We annealed radiated samples #2 and #3 in an oven at 50°C and 80°C for 15 h and 265 h, respectively. The results illustrated that thermal annealing could indeed reduce the radiation-induced DCR. To extend the orbital lifetimes of the APDs, it is valuable to reserve their thermal annealing abilities. Herein, we used an H-bridge pulse-width modulation (PWM) driver to drive the two-stage TEC, which could heat or cool the APD. Accordingly, the APD could operate between approximately + 80°C and −60°C, as required. Note that laser annealing is another way to mitigate radiation-induced damages, as recently reported in [55], but it needs high-power laser sources.

For the Excelitas SLiK APD, the after-pulsing probability decreased exponentially with time and its tail attenuated to very low values after a few hundred nanoseconds [56]. A passive quenching circuit [57,58] with a dead time of a few microseconds was used to drive the APD, which can significantly increase its after-pulsing resistance [52]. The after-pulsing probability before launching was less than 0.1%. Based on these studies [52,54] and our satellite orbit parameters, the expected after-pulsing probability was approximately 0.5% at the end of the two-year mission.

Lastly, we performed a series of space environmental performance and functional tests of the spaceborne, low-noise, Si APD SPDs on the ground, including random vibration, vacuum thermal cycling, and life tests. The results were satisfactory. Typical results are listed in Table 1.

Tables Icon

Table 1. Characteristics of the spaceborne low-noise, Si avalanche photodiode (APD), single-photon detectors (SPDs) in ground tests

3. In-orbit operations and performance

On 16 August 2016, the spaceborne, low-noise, Si APD SPDs were placed into orbit, onboard the QUESS. Figure 7(a)–7(c) show the DCRs and operating parameters of the two SPDs, namely SPD #1 and SPD #2, during the 1029 days of in-orbit operations. Operating parameters of temperature (−40°C) and excess voltages (21 V, 18 V) were first applied to the SPDs to achieve high-detection efficiency and low DCR. As the in-orbit cumulative radiation dose increased, the SPD DCR increased significantly. On the 84th day, the DCR of SPD #1 increased from ∼30 cps to ∼200 cps, and the SPD operating temperature was then lowered to −50°C. The excess voltage was slightly adjusted to mitigate the radiation-induced DCR increment. In the first 262 days, the radiation-induced DCRs were kept below 150 cps, and the detection efficiency remained ∼60% at the wavelength of 780 nm. During this period, ground–satellite quantum teleportation—the practical application of the spaceborne, low-noise Si APD SPDs—was demonstrated [16].

 figure: Fig. 7.

Fig. 7. In-orbit characteristics of SPDs for 1029 consecutive operational days on the QUESS. (a) SPD DCR versus in-orbit time. (b) SPD excess voltage. (c) SPD operating temperature.

Download Full Size | PDF

On the 289th day, the DCR of the SPD #2 increased abruptly from ∼300 cps to ∼1600 cps, and deviated from linearity. This phenomenon can be explained according to Eq. (2). When the energy level of the new radiation-induced carrier generation center was very close the intrinsic Fermi energy level, the carrier generation rate of the new radiation-induced carrier generation center would be very high, and would subsequently lead to the sudden DCR increment. Note that the DCR increase was observed in the ground tests when the radiation dose was low, as shown in Fig. 1. To mitigate the sudden radiation-induced DCR increment of SPD #2, a trade-off was established between the DCR and the DE. Optimization of the receiving end signal-to-noise (SNR) rate further lowered the SPD #2 excess voltage. By adjusting the SPD #2, the excess voltage at the cost of a DE of ∼15%, we increased the receiving end SNR by approximately four times.

In addition, because the QUESS telescope points at a fixed ground station, the SPD DCRs were unavoidably affected by background noise, especially by the moonlight. This result is shown in Fig. 8. The DCRs fluctuated as a function of the moon phase of the ground-station, and the maximum DCR owing to the moonlight was ∼200 cps. It is an effective way to obtain low-background noise by avoiding the effects attributed to the full moon.

 figure: Fig. 8.

Fig. 8. DCRs of SPD #1 and SPD #2 during the first 262 days in orbit. Black and red dots show the observed data. Green and blue curves represent the radiation-induced DCR and moonlight-induced DCR, respectively. SPD #1 and SPD #2 feature a DCR increase of ∼0.56 cps/day during the first 262 days in orbit on the QUESS.

Download Full Size | PDF

The in-orbit average DCR of SPD #1–2 in the past 1029 days was showed in Fig. 9. By deep cooling and SNR optimization, the spaceborne, low-noise Si APD SPDs have achieved an ultralow, in-orbit DCR increment rate of ∼0.54 cps/day, which is two orders of magnitude lower than that of previous results [1,2], while the DE was maintained above 45%.

 figure: Fig. 9.

Fig. 9. In-orbit average DCR of SPD #1 and SPD #2 over a 1029-day period.

Download Full Size | PDF

4. Ground–satellite QKD simulations

Spaceborne, low-noise, Si APD SPDs, play an important role in ground-to-satellite or uplink QKD and its dark count mainly affects the bit error rate of QKD. To generate a secret key, QKD requires the DCR to be kept at a low level to produce a low-quantum bit error rate. Spaceborne, low-noise, Si APD SPDs naturally satisfies this requirement, thus making them suitable for ground-to-satellite QKD. We performed a numerical simulation to study how the DCR affects ground-to-satellite QKD. This involved the adoption of the parameters of an actual satellite-based QKD experiment [59], namely a channel loss of ∼40 dB, a pulse generation rate of 100 MHz, a detection efficiency of 60%, and an optical misalignment of 0.5%. The DCRs were the observed data of average of SPD #1-2 (see Fig. 9), the conventional increase of ∼30 cps/day, and our increase of ∼0.54 cps/day. We considered a decoy-state QKD with security against general attacks in the finite-key setting, in which we chose a total number of pulses $\textrm{N} = {10^{11}}$ and a security parameter $\varepsilon = {10^{ - 7}}$. The simulation results are shown in Fig. 10.

 figure: Fig. 10.

Fig. 10. Expected secret-key rate for ground-to-satellite QKD. The DCRs are (i) those observed in a previous study [16] (blue curve), (ii) those of the spaceborne, low-noise, Si APD SPDs in the present study (green curve), and (iii) those of the previous in-orbit SPDs (red curve). Other simulation parameters are adopted from the satellite-based QKD experiment [16].

Download Full Size | PDF

We found that our SPDs could enable ground-to-satellite QKD for more than 4.6 years (∼1,700 days), whereas a conventional SPD could be used for only one month (∼30 days). Besides QKD, the spaceborne, low-noise, Si APD SPD has already enabled the first demonstration of ground–satellite quantum teleportation [16].

5. Conclusion

In this study, we have presented the results from in-orbit tests of spaceborne, low-noise, Si APD SPDs, based on a special shielding structure, multistage cooling technologies, and special driver electronics. By controlling the close-loop DCR, we achieved a low DCR increase of ∼0.54 cps/day for 1029 days of in-orbit operations, while we detected efficiencies > 45% at 780 nm. As a potential application in addition to ground–satellite quantum teleportation, ground–satellite QKD was numerically simulated, and results indicated that the spaceborne, low-noise, Si APD SPDs could substantially extend the lifetime of the satellite for QKD.

In future work, we will test the method of annealing [54,55] in orbit. Moreover, the active quenching method could be applied to the next generation spaceborne, low-noise, Si APD SPDs, and their operation temperature could be further cooled to a lower temperature, such as −100°C. This may achieve a higher maximum count rate and a lower in-orbit DCR increment. In addition, for the medium and high-orbit satellite, its optoelectronic devices will face high-radiation doses. According to Eq. (4), the APD DCR increment is linearly proportional to the depletion region volume for the same radiation dose. Thus, it is also an effective APD antiradiation measure to use smaller active APD diameters. Furthermore, we have radiated the SPDs with small active diameters, such as the ID Quantique’s ID100-SMF20 with a silicon APD of $20\mu \textrm{m}$, and the $PD050-CBT from Micro Photon Devices with a silicon APD of $50\mu \textrm{m}$. The radiation test results indicated that the small active diameter of the APD led to an extremely low DCR radiation sensitivity.

Our SPDs provided a practical choice for weak optical signal reception in space applications, such as satellite-based quantum applications, deep-space optical communications, and laser time transfers.

Funding

Chinese Academy of Sciences Key Technology Talent Program; National Natural Science Foundation of China (11705191, U1738202); Youth Innovation Promotion Association of the Chinese Academy of Sciences (2017500); Strategic Pioneer Research Projects of Defense Science and Technology.

Disclosures

The authors declare no conflicts of interest.

References

1. X. Sun, M. A. Krainak, J. B. Abshire, J. D. Spinhirne, C. Trottier, M. Davies, H. Dautet, G. R. Allan, A. T. Lukemire, and J. C. Vandiver, “Space-qualified silicon avalanche-photodiode single-photon-counting modules,” J. Mod. Opt. 51(9-10), 1333–1350 (2004). [CrossRef]  

2. X. Sun, J. D. Spinhirne, and M. A. Krainak, “In orbit performance of Si avalanche photodiode single photon counting modules in the Geoscience Laser Altimeter System on ICESat,” Proc. SPIE 6372, 63720P (2006). [CrossRef]  

3. X. Sun, P. L. Jester, J. B. Abshire, and E. S. Chang, “Performance of the GLAS space lidar receiver through its seven-year space mission,” in Lasers and Electro-Optics (CLEO), 2011 Conference on, (IEEE, 2011), 1–2.

4. I. Prochazka and F. Yang, “Photon counting module for laser time transfer via Earth orbiting satellite,” J. Mod. Opt. 56(2-3), 253–260 (2009). [CrossRef]  

5. W. Meng, H. Zhang, Z. Zhang, and I. Prochazka, “The application of single photon detector technique in Laser Time Transfer for Chinese Navigation Satellites,” in SPIE Optics+ Optoelectronics, (International Society for Optics and Photonics, 2013), 87730E.

6. C. Fabjan and H. Schopper, “Detectors for Particles and Radiation. Part 1: Principles and Methods,” Landolt Börnstein (2011).

7. D. M. Boroson, J. J. Scozzafava, D. V. Murphy, B. S. Robinson, and M. Lincoln, “The lunar laser communications demonstration (LLCD),” in Space Mission Challenges for Information Technology, 2009. SMC-IT 2009. Third IEEE International Conference on, (IEEE, 2009), 23–28.

8. D. M. Boroson, B. S. Robinson, D. V. Murphy, D. A. Burianek, F. Khatri, J. M. Kovalik, Z. Sodnik, and D. M. Cornwell, “Overview and results of the lunar laser communication demonstration,” in Free-Space Laser Communication and Atmospheric Propagation XXVI, (International Society for Optics and Photonics, 2014), 89710S.

9. J. Rarity, P. Tapster, P. Gorman, and P. Knight, “Ground to satellite secure key exchange using quantum cryptography,” New J. Phys. 4, 82 (2002). [CrossRef]  

10. C. Bonato, A. Tomaello, V. Da Deppo, G. Naletto, and P. Villoresi, “Feasibility of satellite quantum key distribution,” New J. Phys. 11(4), 045017 (2009). [CrossRef]  

11. J. Bourgoin, E. Meyer-Scott, B. L. Higgins, B. Helou, C. Erven, H. Huebel, B. Kumar, D. Hudson, I. D’Souza, and R. Girard, “A comprehensive design and performance analysis of low Earth orbit satellite quantum communication,” New J. Phys. 15(2), 023006 (2013). [CrossRef]  

12. Z. Tang, R. Chandrasekara, Y. C. Tan, C. Cheng, L. Sha, G. C. Hiang, D. K. Oi, and A. Ling, “Generation and analysis of correlated pairs of photons aboard a nanosatellite,” Phys. Rev. Appl. 5(5), 054022 (2016). [CrossRef]  

13. R. Chandrasekara, Z. Tang, Y. C. Tan, C. Cheng, A. Ling, and D. K. Oi, “Quantum Physics with Cubesats: In-Orbit Observation of Photon Pair Correlations on Board the Galassia Spacecraft,” (2016).

14. S. K. Joshi, J. Pienaar, T. C. Ralph, L. Cacciapuoti, W. McCutcheon, J. Rarity, D. Giggenbach, V. Makarov, I. Fuentes, and T. Scheidl, “Space QUEST mission proposal: experimentally testing decoherence due to gravity,” arXiv preprint arXiv:1703.08036 (2017).

15. R. Bedington, J. M. Arrazola, and A. Ling, “Progress in satellite quantum key distribution,” npj Quantum Inf 3(1), 30 (2017). [CrossRef]  

16. J. G. Ren, P. Xu, H. L. Yong, L. Zhang, S. K. Liao, J. Yin, W. Y. Liu, W. Q. Cai, M. Yang, L. Li, K. X. Yang, X. Han, Y. Q. Yao, J. Li, H. Y. Wu, S. Wan, L. Liu, D. Q. Liu, Y. W. Kuang, Z. P. He, P. Shang, C. Guo, R. H. Zheng, K. Tian, Z. C. Zhu, N. L. Liu, C. Y. Lu, R. Shu, Y. A. Chen, C. Z. Peng, J. Y. Wang, and J. W. Pan, “Ground-to-satellite quantum teleportation,” Nature 549(7670), 70–73 (2017). [CrossRef]  

17. J. Yin, Y. Cao, Y. H. Li, S. K. Liao, L. Zhang, J. G. Ren, W. Q. Cai, W. Y. Liu, B. Li, H. Dai, G. B. Li, Q. M. Lu, Y. H. Gong, Y. Xu, S. L. Li, F. Z. Li, Y. Y. Yin, Z. Q. Jiang, M. Li, J. J. Jia, G. Ren, D. He, Y. L. Zhou, X. X. Zhang, N. Wang, X. Chang, Z. C. Zhu, N. L. Liu, Y. A. Chen, C. Y. Lu, R. Shu, C. Z. Peng, J. Y. Wang, and J. W. Pan, “Satellite-based entanglement distribution over 1200 kilometers,” Science 356(6343), 1140–1144 (2017). [CrossRef]  

18. W. Wootters, “Wootters, WK, and WH Zurek, 1982, Nature (London) 299, 802,” Nature (London) 299(5886), 802–803 (1982). [CrossRef]  

19. D. Shin, A. Kirmani, V. K. Goyal, and J. H. Shapiro, “Photon-efficient computational 3-D and reflectivity imaging with single-photon detectors,” IEEE Trans. Comput. Imaging 1(2), 112–125 (2015). [CrossRef]  

20. A. Swatantran, H. Tang, T. Barrett, P. DeCola, and R. Dubayah, “Rapid, high-resolution forest structure and terrain mapping over large areas using single photon lidar,” Sci. Rep. 6(1), 28277 (2016). [CrossRef]  

21. W. G. Lawrencea, G. Varadia, G. Entinea, E. Podniesinskib, and P. K. Wallaceb, “A comparison of avalanche photodiode and photomultiplier tube detectors for flow cytometry,” in Proc. of SPIEVol, (2008), 68590M–68591.

22. G. Kell, A. Bülter, M. Wahl, and R. Erdmann, “τ-SPAD: a new red sensitive single-photon counting module,” in Advanced Photon Counting Techniques V, (International Society for Optics and Photonics, 2011), 803303.

23. R. H. Hadfield, “Single-photon detectors for optical quantum information applications,” Nat. Photonics 3(12), 696–705 (2009). [CrossRef]  

24. H. Dautet, P. Deschamps, B. Dion, A. D. MacGregor, D. MacSween, R. J. McIntyre, C. Trottier, and P. P. Webb, “Photon counting techniques with silicon avalanche photodiodes,” Appl. Opt. 32(21), 3894–3900 (1993). [CrossRef]  

25. X. Sun and H. Dautet, “Proton radiation damage of Si APD single photon counters,” in Radiation Effects Data Workshop, 2001 IEEE, (IEEE, 2001), 146–150.

26. B. Schutz, H. Zwally, C. Shuman, D. Hancock, and J. DiMarzio, “Overview of the ICESat mission,” Geophys. Res. Lett. 32(21), L21S01 (2005). [CrossRef]  

27. J. Kataoka, T. Toizumi, T. Nakamori, Y. Yatsu, Y. Tsubuku, Y. Kuramoto, T. Enomoto, R. Usui, N. Kawai, and H. Ashida, “In-orbit performance of avalanche photodiode as radiation detector on board the picosatellite Cute1.7 + APD II,” J. Geophys. Res. Atmos. 115(A5), 1292–1300 (2010). [CrossRef]  

28. P. Jianwei, “Progress of the quantum experiment science satellite (QUESS) Micius project,” 空间科学学报 38, 604–609 (2018).

29. M. Raska, “China's Quantum Satellite Experiments: Strategic and Military Implications,” (2016).

30. C. H. Bennett and G. Brassard, “Experimental quantum cryptography: the dawn of a new era for quantum cryptography: the experimental prototype is working,” ACM Sigact News 20(4), 78–80 (1989). [CrossRef]  

31. C.-Z. Peng, T. Yang, X.-H. Bao, J. Zhang, X.-M. Jin, F.-Y. Feng, B. Yang, J. Yang, J. Yin, and Q. Zhang, “Experimental free-space distribution of entangled photon pairs over 13 km: towards satellite-based global quantum communication,” Phys. Rev. Lett. 94(15), 150501 (2005). [CrossRef]  

32. R. Ursin, F. Tiefenbacher, T. Schmitt-Manderbach, H. Weier, T. Scheidl, M. Lindenthal, B. Blauensteiner, T. Jennewein, J. Perdigues, and P. Trojek, “Entanglement-based quantum communication over 144 km,” Nat. Phys. 3(7), 481–486 (2007). [CrossRef]  

33. J. Yin, J.-G. Ren, H. Lu, Y. Cao, H.-L. Yong, Y.-P. Wu, C. Liu, S.-K. Liao, F. Zhou, and Y. Jiang, “Quantum teleportation and entanglement distribution over 100-kilometre free-space channels,” Nature 488(7410), 185–188 (2012). [CrossRef]  

34. X.-S. Ma, T. Herbst, T. Scheidl, D. Wang, S. Kropatschek, W. Naylor, B. Wittmann, A. Mech, J. Kofler, and E. Anisimova, “Quantum teleportation over 143 kilometres using active feed-forward,” Nature 489(7415), 269–273 (2012). [CrossRef]  

35. S. K. Liao, W. Q. Cai, W. Y. Liu, L. Zhang, Y. Li, J. G. Ren, J. Yin, Q. Shen, Y. Cao, Z. P. Li, F. Z. Li, X. W. Chen, L. H. Sun, J. J. Jia, J. C. Wu, X. J. Jiang, J. F. Wang, Y. M. Huang, Q. Wang, Y. L. Zhou, L. Deng, T. Xi, L. Ma, T. Hu, Q. Zhang, Y. A. Chen, N. L. Liu, X. B. Wang, Z. C. Zhu, C. Y. Lu, R. Shu, C. Z. Peng, J. Y. Wang, and J. W. Pan, “Satellite-to-ground quantum key distribution,” Nature 549(7670), 43–47 (2017). [CrossRef]  

36. E. Stassinopoulos and J. P. Raymond, “The space radiation environment for electronics,” Proc. IEEE 76(11), 1423–1442 (1988). [CrossRef]  

37. W. Schimmerling, “The Space Radiation Environment: An Introduction,” Science@NASA. (2012).

38. A. Globus and J. Strout, “Orbital Space Settlement Radiation Shielding,” preprint, issued July 2015 available on-line at space. alglobus. net (2015).

39. J. Kodet, I. Prochazka, J. Blazej, X. Sun, and J. Cavanaugh, “Single photon avalanche diode radiation tests,” Nucl. Instrum. Methods Phys. Res., Sect. A 695, 309–312 (2012). [CrossRef]  

40. G. Lindström, “Radiation damage in silicon detectors,” Nucl. Instrum. Methods Phys. Res., Sect. A 512(1-2), 30–43 (2003). [CrossRef]  

41. J. Srour, C. J. Marshall, and P. W. Marshall, “Review of displacement damage effects in silicon devices,” IEEE Trans. Nucl. Sci. 50(3), 653–670 (2003). [CrossRef]  

42. A. Akkerman, J. Barak, M. Chadwick, J. Levinson, M. Murat, and Y. Lifshitz, “Updated NIEL calculations for estimating the damage induced by particles and γ-rays in Si and GaAs,” Radiat. Phys. Chem. 62(4), 301–310 (2001). [CrossRef]  

43. I. Jun, M. A. Xapsos, S. R. Messenger, E. A. Burke, R. J. Walters, G. P. Summers, and T. Jordan, “Proton nonionizing energy loss (NIEL) for device applications,” IEEE Trans. Nucl. Sci. 50(6), 1924–1928 (2003). [CrossRef]  

44. M. Kruglanski, N. Messios, E. De Donder, E. Gamby, S. Calders, L. Hetey, and H. Evans, “Space environment information system (SPENVIS),” in EGU General Assembly Conference Abstracts, 7457 (2009).

45. W. Kindt and H. Van Zeijl, “Modelling and fabrication of Geiger mode avalanche photodiodes,” IEEE Trans. Nucl. Sci. 45(3), 715–719 (1998). [CrossRef]  

46. H. N. Becker and A. H. Johnston, “Dark current degradation of near infrared avalanche photodiodes from proton irradiation,” (2004).

47. H. N. Becker, T. F. Miyahira, and A. H. Johnston, “The influence of structural characteristics on the response of silicon avalanche photodiodes to proton irradiation,” IEEE Trans. Nucl. Sci. 50(6), 1974–1981 (2003). [CrossRef]  

48. J. Srour, S. Chen, S. Othmer, and R. Hartmann, “Radiation damage coefficients for silicon depletion regions,” IEEE Trans. Nucl. Sci. 26(6), 4783–4791 (1979). [CrossRef]  

49. O. Zeynali, D. Masti, A. Ebrahimi, M. Orvatinia, and N. Maryam, “The design and simulation of the shield reduce ionizing radiation effects on electronic circuits in satellites,” Electr. Electron. Eng. 1, 112–116 (2011).

50. J. Wilson, F. Cucinotta, J. Miller, J. Shinn, S. Thibeault, R. Singleterry, L. Simonsen, and M. Kim, “Materials for shielding astronauts from the hazards of space radiations,” Mater. Res. Soc. Symp. Proc. 551, 3 (1998). [CrossRef]  

51. N. Nightingale, “A new silicon avalanche photodiode photon counting detector module for astronomy,” Exp. Astron. 1(6), 407–422 (1990). [CrossRef]  

52. Y.-S. Kim, Y.-C. Jeong, S. Sauge, V. Makarov, and Y.-H. Kim, “Ultra-low noise single-photon detector based on Si avalanche photodiode,” Rev. Sci. Instrum. 82(9), 093110 (2011). [CrossRef]  

53. Y. C. Tan, R. Chandrasekara, C. Cheng, and A. Ling, “Silicon avalanche photodiode operation and lifetime analysis for small satellites,” Opt. Express 21(14), 16946 (2013). [CrossRef]  

54. E. Anisimova, B. L. Higgins, J.-P. Bourgoin, M. Cranmer, E. Choi, D. Hudson, L. P. Piche, A. Scott, V. Makarov, and T. Jennewein, “Mitigating radiation damage of single photon detectors for space applications,” EPJ Quantum Technology 4(1), 10 (2017). [CrossRef]  

55. J. G. Lim, E. Anisimova, B. L. Higgins, J.-P. Bourgoin, T. Jennewein, and V. Makarov, “Laser annealing heals radiation damage in avalanche photodiodes,” EPJ Quantum Technology 4(1), 11 (2017). [CrossRef]  

56. A. W. Ziarkash, S. K. Joshi, M. Stipčević, and R. Ursin, “Comparative study of afterpulsing behavior and models in single photon counting avalanche photo diode detectors,” Sci. Rep. 8(1), 5076 (2018). [CrossRef]  

57. S. Cova, M. Ghioni, A. Lacaita, C. Samori, and F. Zappa, “Avalanche photodiodes and quenching circuits for single-photon detection,” Appl. Opt. 35(12), 1956–1976 (1996). [CrossRef]  

58. S. Cova, M. Ghioni, A. Lotito, I. Rech, and F. Zappa, “Evolution and prospects for single-photon avalanche diodes and quenching circuits,” J. Mod. Opt. 51(9-10), 1267–1288 (2004). [CrossRef]  

59. Y. Li, S.-K. Liao, X.-L. Chen, W. Chen, K. Cheng, Y. Cao, H.-L. Yong, T. Wang, H.-Q. Yang, and W.-Y. Liu, “Space-bound optical source for satellite-ground decoy-state quantum key distribution,” Opt. Express 22(22), 27281–27289 (2014). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (10)

Fig. 1.
Fig. 1. Dark count rate (DCR) versus displacement damage dose. During radiation, sample #2 was cooled to −35°C and was biased with an excess voltage of 17 V, while samples #3 and #7 were unbiased at 25°C, and the other samples were cooled to −12°C and biased with an excess voltage of 17 V. For DCR measurements, all samples were cooled to −12°C and biased with an excess voltage of 17 V.
Fig. 2.
Fig. 2. Displacement damage dose versus spherical Al shield radius.
Fig. 3.
Fig. 3. Shielding and installation of the spaceborne, low-noise, Si APD, single photon detectors.
Fig. 4.
Fig. 4. Deep cooling effect on DCR of sample #4. The results for the other six samples are the same. S represents the slope for the DCR versus excess voltage plots at different temperatures.
Fig. 5.
Fig. 5. Schematic of the driving circuit for single-photon detectors (SPDs) (ADC, analog-to-digital-converter; DAC, digital-to-analog-converter; COMP, comparator; DFF, D Flip-flop; PWM, pulse-width modulation; TEC, thermal electric cooler, HV, high voltage. The ADCs, DACs, and PWMs, are connected to a field-programmable gate array (FPGA) in which the temperature and the high-voltage (HV) closed loop control were carried out.
Fig. 6.
Fig. 6. Teleportation fidelity at different excess voltages and excess DCR voltages slopes.
Fig. 7.
Fig. 7. In-orbit characteristics of SPDs for 1029 consecutive operational days on the QUESS. (a) SPD DCR versus in-orbit time. (b) SPD excess voltage. (c) SPD operating temperature.
Fig. 8.
Fig. 8. DCRs of SPD #1 and SPD #2 during the first 262 days in orbit. Black and red dots show the observed data. Green and blue curves represent the radiation-induced DCR and moonlight-induced DCR, respectively. SPD #1 and SPD #2 feature a DCR increase of ∼0.56 cps/day during the first 262 days in orbit on the QUESS.
Fig. 9.
Fig. 9. In-orbit average DCR of SPD #1 and SPD #2 over a 1029-day period.
Fig. 10.
Fig. 10. Expected secret-key rate for ground-to-satellite QKD. The DCRs are (i) those observed in a previous study [16] (blue curve), (ii) those of the spaceborne, low-noise, Si APD SPDs in the present study (green curve), and (iii) those of the previous in-orbit SPDs (red curve). Other simulation parameters are adopted from the satellite-based QKD experiment [16].

Tables (1)

Tables Icon

Table 1. Characteristics of the spaceborne low-noise, Si avalanche photodiode (APD), single-photon detectors (SPDs) in ground tests

Equations (7)

Equations on this page are rendered with MathJax. Learn more.

DCR = S W 1 W 1 + W 2 P pair ( x ) G TOT dx
G T O T = d = 1 D n i τ n exp [ ( E T , d E i ) K T ] + τ p exp [ ( E T , d E i ) K T ]
Δ D C R V = n i ϕ 2 K g n
S = α exp ( β T )
d s = SVe
η d = η ( 1 exp ( V e / V c ) )
F = η l η d F s F l + 1.795 R d t w η l η d + 3.59 R d t w
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.