## Optimal photon pairs for quantum communication protocols

Mikołaj Lasota 1,* and Piotr Kolenderski

1

1 Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toru´ n, Poland

* miklas@umk.pl

## ABSTRACT

We theoretically investigate the problem of finding optimal characteristics of photon pairs, produced in the spontaneous parametric down-conversion (SPDC) process, for fiber-based quantum communication (QC) protocols. By using the accessible setup parameters, the pump pulse duration and the extended phase-matching function width, we minimize the temporal width of SPDC photons within the general scenario. This allows one to perform more effectively the temporal filtering procedure, which aims at reducing the noise acquired by the measurement devices. Moreover, we compare the obtained results with the achievable parameter values for SPDC sources based on b -Barium Borate (BBO) crystal. We also investigate the influence of non-zero detection timing jitter. Finally, we apply our optimization strategy to a simple quantum key distribution scheme to show that the full optimization of an SPDC source can potentially extend the maximal security distance by several tens of kilometres, which is around 30% more as compared to previous approaches.

## Introduction

Quantum communication is a vast field of physical science focused on improving the process of information distribution among spatially separated entities with the use of quantum mechanics. The exploration of various types of quantum correlation and application of fundamental quantum laws has lead to a plethora of proposals for novel communication protocols, including quantum key distribution 1,2 , secret sharing 3 , quantum teleportation 4 and dense coding 5 . However, the initially proposed theoretical versions of these protocols have typically assumed ideal performance of the setup elements required for their physical realization, which is unreachable in practice. As a consequence, the performance of real-life implementations of QC protocols has been severely limited.

Realization of many such schemes requires using sources of single photons or entangled photon pairs. One of the most popular types of them are the devices utilizing the phenomenon of spontaneous parametric down-conversion 6, 7 . They have many advantages, including high quality of the emitted photons 8, 9 , high generation and collection efficiency 10-12 and relatively low cost of their construction. Therefore, they have been extensively used in practical implementations of many QC protocols 13-17 . However, photons born in the SPDC process are not monochromatic. Thus, they propagate through dispersive media ( e.g. the standard telecommunication fibers) with wavelength-dependent velocity. As a consequence, their temporal width, defined as the standard deviation of the probability distribution function for the time of their arrival at the destination point, grows with the length of the utilized dispersive quantum channel. It forces the experimenter to define longer detection windows for the photon measurement system in order not to lose considerable amount of real signals. However, the longer detection windows are, the more noise is registered during the realization of a given QC protocol, negatively affecting its performance.

Taking into account the above consideration, the minimization of the temporal width of the emitted photons after their propagation through telecommunication fibers is very important. Nevertheless, to the best of our knowledge such a general optimization has not yet been done. Only recently it was shown that changing the properties of pairs of photons can indeed significantly influence the performance of quantum protocols 18-20 . A preliminary optimization of SPDC source for a specific setup configuration was also performed 20 . In this manuscript we generalize this optimization to an arbitrary QC scheme. We also investigate the influence of non-zero timing jitter of realistic detectors on the obtained results and discuss whether the theoretically optimal values can be implemented in practice. Finally, we estimate the advantage stemming from the optimization of SPDC source on the maximal security distance of a basic quantum key distribution (QKD) scheme.

In our work we consider SPDC source of photon pairs based on a nonlinear crystal parametrized by the effective phasematching function width s , pumped by laser pulses of temporal width t p . Throughout this manuscript we assume that the information on the time moments at which the laser sends pump pulses to the crystal is available to everyone interested. Photons generated by the SPDC source are subsequently transfered to single-photon detectors through quantum channels of length LA and LB , characterized by the group velocity dispersion (GVD) equal to 2 b A and 2 b B , respectively. This scheme

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-dc3bc231-58c9-4603-8263-9ababc64dcdc.jpg)

Figure 1. Detection scheme. The generic setup configuration for the detection of SPDC photons.

Figure 2. The QKD scheme. A simple scheme of discrete-variable QKD with the source of entangled photon pairs placed outside of Alice's and Bob's laboratories. R denotes polarization rotators.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-fee90a39-488f-48dc-bf92-6c0b1bf2ad28.jpg)

Figure 3. The dependence of t Ah on the pump laser. The temporal width of the heralded photon A, t Ah , plotted as a function of the duration time of a pump laser pulse for the case when two standard single-mode fibers (SMFs) of length LA = 1km and (a) LB = 1km, (b) LB = 100km are placed between the crystal and single-photon detectors. The legend corresponds to both panels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-12865441-f14c-4e79-9b33-a82395b8a543.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-c9092bb6-6a5a-4f21-89f8-95acc004db43.jpg)

is illustrated in Fig. 1. We use the following notation to simplify the analytical calculations: DX ≡ b X LX (for X = A B , ). To show the potential to improve the performance of QC schemes by optimizing the SPDC source we also consider a basic setup configuration for the realization of BB84 protocol in the entanglement-based variant, presented in Fig. 2. All the subsequent figures shown in this manuscript are made for b A = b B = -1 15 . × 10 -26 s 2 / m. For the standard wavelength of 1550nm it corresponds 21 to the dispersion value of 18ps / ( nm km , that is typical for single-mode fibers (SMFs) · ) 22 , which are the most common telecommunication channels utilized in practical QC schemes. During the QKD security analysis we also assume the typical value of their attenuation coefficient: a A = a B = 0 2dB km. . /

## Results

## Optimization of temporal widths

As has already been stated in the Introduction, in the context of QC applications it is desirable for the temporal widths of SPDC photons to be as narrow as possible. Therefore, a natural question is: what are the optimal values of the source parameters, t p and s , for which the temporal widths of SPDC photons, written explicitly in the Methods section (formulas (3) and (4)), are the lowest? In practice it is much easier to calibrate the temporal width of pump laser pulses than to modify the effective phase-matching function for the nonlinear crystal, since the latter usually requires replacing the crystal itself. Therefore, we first consider the situation in which the experimenter can only change the pump laser utilized by the SPDC source, while the crystal is fixed.

In this case the temporal width of the unheralded photon A, t A , reaches its lowest value, equal to t low A =( + 2 | DA | s 2 ) / ( 2 s ) , for t p = √ 2 | DA | . Since t A dos not depend on DB , the above result is identical for the symmetric and asymmetric setup configurations. In the symmetric case also the temporal width of the heralded photon A, t Ah , reaches its minimum for the same value of t p . It reads: t low sym , Ah =[ 2 | DA | ( D A 2 s 4 + ) ( 4 / | DA | s 2 + ) ] 2 2 1 / 2 . On the othe hand, the optimization of t Ah over t p for

Figure 4. The dependence of t Ah on the pump laser and the nonlinear crystal. Logarithm of the temporal width of the heralded photon A, t Ah , at the entrance to the detector, plotted as a function of the duration time of the pump laser pulse, t p , and the effective phase-matching function width of the nonlinear crystal, s , for the case when the source is connected with the detectors by SMFs of length LA = 1km and (a) LB = 1km or (b) LB = 100km. The contours shown in both plots represent values from log 10 t Ah = -11 2 to log . 10 t Ah = -9 2 with constant 0 2 spacing. . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-7ecded7a-1b48-4199-a04b-a9d369dcdd0d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-881a241d-afbf-4f35-ab50-a2dabc17e5b8.jpg)

the asymmetric setup configuration is much more complicated. In this case the function t Ah ( t p ) does not always have a global minimum and the conditions for its existence heavily depend on the relationship between DA , DB and s . The details of this dependence can be found in the Methods section, along with the derivation of the above formulas for t low A and t low sym , Ah .

The examples of the relationship t Ah ( t p ) for the symmetric and highly asymmetric setup configurations can be seen in Fig. 3. In the symmetric case, presented in panel (a), the function t Ah ( t p ) has a well-defined minimum for any s . Its value depends on the effective phase-matching function width relatively weekly, while the value of t p for which this minimum is reached is independent of s . The situation is much different in the highly asymmetric case, illustrated in panel (b). Here both the optimal value of t p and the minimal value of t Ah significantly depend on the effective phase-matching function width. Moreover, for s = 100GHz (corresponding to 0 13nm at 1550nm in terms of wavelength) none of the conditions for the . existence of the global minimum of t Ah is fulfilled. In this situation t Ah ( t p ) is monotonically increasing function (see the orange dot-dashed line). It can also be seen that in the asymmetric scenario the comparison between the functions of t Ah ( t p ) plotted for different s heavily depends on t p . For example, while for very short pump pulses the value of t Ah calculated for s = 1THz is much smaller than for s = 10GHz, it is the opposite for large t p . The situation like this cannot be seen in the symmetric case.

Contrary to the scenario when the nonlinear crystal is fixed, full optimization of a SPDC source over the parameters t p and s cannot be done analytically in the general case. Nevertheless, it can be performed for the symmetric setup configuration, when DA = DB ≡ D . This task has already been done in our previous paper 18 , where it was shown that in the symmetric case the optimal values of the SPDC source parameters are t sym p = √ 2 | D | and s sym = √ 2 / | D | . For these numbers the function t Ah reaches its absolute minimum, equal to t sym Ah = √ 2 | D | . In the symmetric case t Ah exhibits high symmetry both as a function of t p and s . It can be seen in Fig. 4(a), where the temporal width of the photon A is plotted for LA = LB = 1km. For comparison, in Fig. 4(b) we plot t Ah for the highly asymmetric case of LA = 1km and LB = 100km.

## Dependence on the length of the heralding arm

For the asymmetric QC scheme it is possible to reduce the temporal width of SPDC photons propagated through one of its arms by introducing a proper amount of dispersion to the other arm ( e.g. by adjusting its length). This can have positive effect on the performance of QC protocols in some setup configurations, as has already been shown in the context of the asymmetric QKD scenario 19 . However, the framework used in the aforementioned work was based on the analysis of spectral correlation and generated photon widths. It gave general insight into the physical mechanisms yielding the optimal performance of the QKD scheme, but must be reformulated to be directly related to the typical experimental scenario. Here we accomplish this goal using the parameters s and t p .

Figure 5. The dependence of t Ah on nonlinear crystal and the length of the heralding arm. Logarithm of the temporal width of the heralded photon A at the entrance to the detector, t Ah , shown as a function of the length of the heralding SMF quantum channel, LB , and the effective phase-matching function width of the nonlinear crystal, s , plotted for the case when the source is connected with the detector A by another SMF quantum channel of length LA = 1km. For every pair of values ( LB , s ) the calculated temporal width has been optimized over the pump laser pulse duration t p . The overshadowed area near the left [right] edge of the figure corresponds to the range of ( LB , s ) , for which the optimal value of t Ah is reached for t p → ¥ [ t p → 0]. For other combinations of LB and s the optimal value of t p is given by the formula (7). The spacing between the neighboring contours becomes smaller than 0 2 for log . 10 t Ah &lt; -11 2 in order to better illustrate how this function behaves . near its minimum.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-82c1fa49-7798-4557-a936-7b5cf634097f.jpg)

In Fig. 5 it can be seen how the temporal width t Ah , optimized over the duration time of the pump laser pulses, changes with different values of LB and s . This picture shows that the conditions for the function t Ah ( t p ) to have a global minimum, derived in the Methods section, are always fulfilled when the channel parameters, DA and DB , are of the same order of magnitude. Only for highly asymmetric schemes this function can be minimized asymptotically for t p → ¥ or t p → 0. Furthermore, it is interesting to note that if the effective phase-matching function width is significantly larger than its optimal value, s opt , extending the length of the heralding arm always leads to the reduction of t Ah . This effect is more prominent for smaller distances, while for LB → ¥ the temporal width of the heralded photon asymptotically decreases to a fixed value. On the other hand, when s &lt; s opt , extending the heralding arm has the opposite effect on t Ah to the one described above. Therefore, one can conclude that for the asymmetric QKD scheme the maximal security distance in one arm can be extended by introducing more dispersion to the other arm, as long as s &gt; s opt . This conclusion is similar to the one stated in our previous work 19 . However, it is important to underline that contratry to the aforementioned paper, here the global time reference, i.e. the timing information on the pump pulses, is assumed to be known by Alice and Bob, as has already been stated in the Introduction. Consequently, it can potentially have much broader practical application than it was previously thought 19 .

Since the simplest way of introducing more dispersion to the heralding arm is to use longer telecommunication fiber, one should be aware of the fact that such action would always reduce the probability of registering the heralding photon. Therefore, one should avoid it as long as the length of the heralded arm is short enough to provide the QKD security even without introducing additional dispersion, as it would unnecessarily descrease the key generation rate. Only when the length of the heralded arm is indeed too long for the security of the traditional setup configuration, lenghtening the heralding fiber may provide positive results for the participants of the QKD protocol.

## Dependence on the detector timing jitter

It can be seen in Fig. 3-Fig. 5 that if the SMF connecting the source with the detector A is of the order of 1km, the temporal width t Ah can be reduced even below the level of 10ps. This value is comparable with the timing jitter of the best currently existing single-photon detectors 23-26 . In order to estimate the range of fiber lengths for which non-zero jitter can have significant

Figure 6. The influence of non-zero timing jitter on .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-cd4ca4c5-bf1d-41ff-b232-78c90e8bf5d3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-458b7bac-3465-4e1d-a2b9-1bf72945146b.jpg)

t Ah Temporal width of the heralded photon, t J Ah , given by the formula (8),optimized over the SPDC source parameters s and t p , plotted as a function of the length of SMFs separating the source and the photon detection systems in the case of symmetric setup configuration. The detectors' timing jitter is defined as the standard deviation of the detection time probability function.

Figure 7. Realistic values of s for BBO crystal. The effective phase-matching function width, s , calculated for 775nm → 1550nm + 1550nm type I SPDC process, plotted as a function of the angle a between the central propagation directions for the pump photons and the generated photons. The plots are made for the crystal length equal to L cryst = 1cm (solid lines) and L cryst = 1mm (dashed lines), and for the following widths of the transverse spatial modes collected by the SMFs, Wf : 10 m m(black lines), 100 m m (red lines), 1mm (yellow lines). Blue dotted (green dot-dashed) line correspond to the optimal value of s , calculated for the symmetric QC setup with SMF quantum channels of 100km (1km) length.

influence on the temporal widths of SPDC photons, in Fig. 6 we compare the temporal widths of the heralded photon A, t J Ah , optimized over the source parameters, s and t p , calculated as a function of the propagation distance for a few different values of the timing jitter t JA and t JB . The mathematical formula for t J Ah is derived in the Methods section. The plots in Fig. 6 are made for the symmetric QC scheme. As one can see there, if LA and LB are shorter than a few kilometers, the jitter significantly influences t J Ah even if it is much smaller than in the case of the state-of-the-art single-photon detectors. Therefore, to exploit the full potential of the optimization method presented in this paper further development of photon detection technology will be needed. At present, however, it is certainly possible to make the influence of detection jitter negligible if the propagation distance is of the order of tens of kilometers or more. To conclude, the results of our investigation, presented in Fig. 6, indicate that if the experimenter wants to fully optimize the short-distance QC scheme, the detection jitter of realistic single-photon detectors can become one of the most important factors. On the other hand, for long-distance communication schemes the jitter can be safely neglected.

## Realistic values of the effective phase-matching function width

In principle, in the case of any specific QC setup configuration, using the optimization rules presented in this paper allows the experimenter to easily find the most favourable values of a pump laser pulse duration and an effective phase-matching function width. However, one may wonder if these theoretically optimal values would be achievable for realistic SPDC sources. It is much easier to answer this question in the context of the pump laser pulse duration, owing to the variety of commercially available lasers, ranging from the CW to femtosecond ones. Since the optimal value of t p generally grows with the propagation distance and already for LA = LB = 1 m it is approximately equal to 150fs, one can safely say that the theoretically optimal pump laser pulse duration should be achievable for basically every realistic QC scheme.

Performing similar analysis in the context of the effective phase-matching function width associated with different kinds of nonlinear crystals is much more complexed. The value of s depends not only on the type of nonlinear material, but also on several other parameters such as the crystal length or its optical axis orientation 27 . However, in order to get some intuition in this matter, we analyzed here a specific case of BBO crystal cut for degenerate type I SPDC process, in which 775nm pump photons are converted to pairs of 1550nm photons. The results of our investigation are presented in Fig. 7, where the effective phase-matching function width was plotted as a function of the angle a between the central propagation directions of the pump photons and the generated photons. The calculations were made for several different values of the crystal length, L cryst , and the

Figure 8. Key generation rate. The lower bound for the key generation rate, K , plotted as a function of the length LA of the SMF used to connect the SPDC source with the laboratory of Alice, LA , for the BB84 protocol performed in (a) symmetric, (b) asymmetric version of the QKD scheme presented in Fig. 2. In panel (a) the plots are calculated for the following values of the source parameters: s = 1THz and t p = 1ns (dashed yellow line), s = 1THz and t p = √ 2 | b | L (dotted red line), s = √ 2 / ( | b | L ) and t p = √ 2 | b | L (solid black line). The three aforementioned curves are drawn assuming ideal single-photon detectors with no timing jitter. Additionally, dot-dashed black line illustrates the lower bound for the key generation rate calculated in the case when the jitter of all the detectors utilized by Alice and Bob is 100ps, while the source parameters are s = √ 2 / ( | b | L ) and t p = √ 2 | b | L . In panel (b) all of the plots are made for s = 1THz, while the values of t p are numerically optimized and the jitter is assumed to be zero.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-c9b3d0b5-223b-4996-a00f-147e6a323770.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00989v2-figure-5c21d676-2033-44e6-abb2-c0878c071e8e.jpg)

width of transverse spatial mode collected by the SMFs, Wf . Additionally, the optimal values of s for symmetric QC setup configuration using SMFs of length L = 1km and L = 100km were indicated in this figure for comparison.

The most important conclusion that can be drawn from Fig. 7 is that for the source based on BBO crystal analyzed here the theoretically optimal values of the effective phase-matching function width can be very difficult to obtain in most practical situations. This goal seems to be especially hard to achieve for the case of collinear source configuration, i.e. when a = 0, which is often the most convenient one in practice. In this situation, even when using exceptionally long BBO crystals, one may hope to obtain s opt width only for short-distance QC schemes. In principle, smaller widths of the effective phase-matching function can be get when the values of Wf are sufficiently large and the BBO crystal is cut to emit pairs of photons at broad angle from the direction of propagation of the pump laser pulses. However, this kind of SPDC source would be significantly more difficult to construct. Moreover, its pair production efficiency and heralding efficiency would most likely be much smaller than for the case of collinear configuration. This would negatively affect the performance of many QC protocols 28 .

## An example of application: quantum key distribution

The potential of the presented method for the optimization of a SPDC source for its use in QC applications can be seen in Fig. 8(a), where we plot the lower bound for the key generation rate that can be obtained from the realization of BB84 protocol in the symmetric version of the setup configuration schematically illustrated in Fig. 2. We analyzed the cases of (i) non-optimized source with s = 1THz and t p = 1ns, (ii) the source with the same s , but optimized over the value of t p and (iii) the fully optimized SPDC source. The value of the pump laser pulse duration used in the case (i) is roughly the same as in one of the rare experimental realizations of long-distance QKD with SPDC sources 29 . Since the authors of the aforementioned paper did not calculate the effective phase-matching function width for their crystal, we decided to use a typical value here. Technical details of the security analysis can be found in the Methods section.

It can be seen in Fig. 8(a) that in principle by fully optimizing the source the maximal security distance for the analyzed scheme can be extended by almost sixty kilometers for each of the two existing quantum channels, which is around 30% compared to the non-optimized case. Moreover, even partial optimization of the source, just over the pump laser pulse duration, can provide the legitimate participants of the BB84 protocol with about 20% of additional security distance. It is also important to notice, that the results plotted in Fig. 8(a) do not change considerably if we assume that Alice and Bob use single-photon detectors characterized by detection timing jitter of t JA = t JB = 100ps, which is well above the best achievable value for the modern devices 23-25 . In this situation the maximal security distance is shortened only by a few kilometers compared to the case with ideal single-photon detectors. This result is consistent with the jitter influence analysis presented earlier in this work.

The results shown in Fig. 8(a) were obtained for the security analysis of the symmetric version of the QKD setup. In Fig. 8(b) we present the results concerning more general situation, in which the two SMFs connecting the SPDC source with Alice and Bob are not of the same length. In this part of our work we specifically focus on checking how changing the length of

Bob's fiber can influence the maximal security distance between the source and Alice. This investigation is motivated by the possibility of decreasing the temporal width of the heralded SPDC photon by extending the distance between the source and the heralding detector, discussed before. While such possibility is available only when s &gt; s opt (see Fig. 5 (a)), the plots shown in Fig. 7 strongly suggest that this requirement can be fulfilled in most practical situations. As can be seen in Fig. 8(b), the maximal security distance between the source and Alice can be increased by several tens of kilometers just by optimizing the value of LB for any given LA , instead of fixing it on some short level.

## Discussion

In this work we performed theoretical optimization of SPDC photon pairs for QC schemes with two dispersive quantum channels of arbitrary lengths. It was done over the pump laser pulse duration and the effective phase-matching function width of nonlinear crystal. We derived an analytical formula for the best setting of the pump laser for a given crystal in the most general case. Moreover, we performed full numerical optimization of a SPDC source, demonstrating the possibility to further refine the performance of quantum protocols. We also showed that the temporal width of a SPDC photon can be minimized in one of two possible ways, depending on the exact value of the effective phase-matching function width: either by increasing the dispersion in the quantum channel or by decreasing it. The first (second) of these possibilities is available when the effective phase-matching function width is larger (smaller) than its optimal value.

To compare theoretical predictions of our work with capabilities of realistic SPDC sources we investigated the source based on BBO crystal, designed for type I SPDC process generating pairs of 1550nm photons. For such source we performed analytical estimation of the effective phase-matching function width. It should be noted here that precise calculation of this parameter can be done only numerically and is beyond the scope of this analysis. The obtained results suggests that for most QC schemes the achievable value of the effective phase-matching function width would be significantly larger than the theoretically optimal one. While in some situations the optimal value could be achieved, it would often require relatively large angles between the pump laser pulse direction and the propagation directions of the generated photons. However, such setup configuration would negatively affect the efficiency of SPDC source. The above consideration raises the question in what situations it would be more beneficial to abandon the full optimization of the SPDC source based on BBO crystal and use the collinear configuration to produce pairs of photons, and when it would be better to push for the full optimization at the expense of efficiency of the source. Further analysis of this problem would be required to reliably answer such question. Moreover, similar investigation performed for other types of nonlinear crystals would be necessary.

To demonstrate the potential for improving the performance of QC protocols by optimizing SPDC source, we analyzed simple entanglement-based QKD scheme. We showed that the maximal secure communication distance can be significantly increased just by properly adjusting the pump laser. Furthermore, if the full optimization of the source is possible, the improvement may even reach 30% in comparison with the practical non-optimized scheme. We also showed that in realistic cases the detection timing jitter reduces the maximal security distance by no more than a few kilometers.

The method for improving the QKD security presented in this paper is based on the reduction of the amount of noise registered by single-photon detectors during the protocol. Therefore, its effectiveness could be significantly smaller if the decrease of the key generation rate to zero at the maximal security distance was mainly caused by some other factors than the reduction of signal-to-noise ratio below the critical level. Specifically, if the SPDC source is used in a prepare-and-measure type of QKD setup configuration, the security of such scheme could be higly dependent on the probability for producing more than one pair of photons, which is always non-zero in realistic situations. However, the damaging influence of the multipair generation events on QKD security can be efficiently reduced by using decoy-pulse method 30 , which greatly limits the possibility to attack multiphoton pulses by a potential eavesdropper. While most of the recent record-breaking long-distance realizations of QKD protocols reported in the literature have been implemented using weak coherent pulses and decoy-pulse method 31-33 , many papers suggest that heralded single-photon sources could potentially be better for this task 34-36 . This notion can be supported by taking into account the recent advances in the field of heralding efficiency of the SPDC sources 10-12 . While the strong temporal broadening of the generated signals has always been a serious obstacle for using these sources in fiber-based long-distance communication, their optimization method presented here allows to overcome this important problem.

Taking into account the above considerations, the noise registered by the measurement systems during the key generation procedure appears to be much bigger issue for long-distance QC than the aforementioned imperfection of photon pair sources. Since in our work we considered dark counts as the only source of noise, one can expect that the SPDC source optimization method can provide even better results in more realistic cases. It seems to be especially promising for the QKD performed in commercial fibers populated by strong classical signals, where the level of channel noise caused by those signals is typically very high 37 . Our results can be particularly useful in the case of asymmetric QKD scheme in which the distance between one of the parties and the source is relatively small and the goal is to maximally extend the security length of the quantum channel connecting the source with the other party. A good example of such scenario can be found when considering a communication between a single individual user and a distant node in a multilevel quntum network with several access networks connected to

the central backbone 38,39 . Then, the maximal security distance between two separate access networks could be substantially increased by introducing more dispersion to the quantum channels connecting the individual users with their respective central nodes, as we also demonstrated in this work.

## Methods

## Temporal widths of SPDC photons

The spectral wavefunction of the pairs of photons produced by an SPDC source can be written in the following approximate form 27,40 :

$$\phi ( v _ { A }, v _ { B } ) \approx \exp \left ( - \frac { \left ( v _ { A } - v _ { B } \right ) ^ { 2 } } { \sigma ^ { 2 } } - \frac { \left ( v _ { A } + v _ { B } \right ) ^ { 2 } \tau _ { p } ^ { 2 } } { 4 } \right ),$$

where n A , n B are frequency detunings from the respective central frequencies. To calculate the temporal wavefunction of the pair of SPDC photons after their propagation through the dispersive media we utilize the following formula 18 :

$$\psi _ { D _ { A } D _ { B } } ( t _ { A }, t _ { B } ) = \frac { 1 } { 4 \pi i \sqrt { D _ { A } D _ { B } } } \int d t _ { A } ^ { \prime } \int d t _ { B } ^ { \prime } \exp \left ( \frac { i ( t _ { A } - t _ { A } ^ { \prime } ) ^ { 2 } } { 4 D _ { A } } + \frac { i ( t _ { B } - t _ { B } ^ { \prime } ) ^ { 2 } } { 4 D _ { B } } \right ) \psi \left ( t _ { A } ^ { \prime }, t _ { B } ^ { \prime } \right ),$$

where y ( t ′ , t ′ ) denotes the initial temporal wavefunction. It can be obtained from f n ( A , n B ) through Fourier transform.

A B Without any loss of generality we focus on calculating the temporal width of the photon entering the detector A (photon A) in Fig. 1. If an experimenter knows nothing about the detection time of the other photon (photon B), the probability distribution function for the detection time of photon A can be calculated as the marginal distribution pA ( t A ) = ∫ d t B | y DADB ( t A , t B ) | 2 . In this case the temporal width of photon A reads:

$$\tau _ { A } = \frac { \sqrt { \left ( \tau _ { p } ^ { 2 } + D _ { A } ^ { 2 } \sigma ^ { 2 } \right ) \left ( 4 + \sigma ^ { 2 } \tau _ { p } ^ { 2 } \right ) } } { 2 \sigma \tau _ { p } }$$

On the other hand, if the detection time of photon B is known to be TB , the probability distribution function for the detection time of photon A takes the form of pAh ( t A , t B = TB ) = | y DADB ( t A , t B = TB ) | 2 / [ ∫ d t A | y DADB ( t A , t B = TB ) | 2 ] . Its temporal width is then given by:

$$\tau _ { A h } = \sqrt { \frac { 1 6 \left ( \tau _ { p } ^ { 2 } - D _ { A } D _ { B } \sigma ^ { 2 } \right ) ^ { 2 } + \left ( D _ { A } + D _ { B } \right ) ^ { 2 } \left ( \sigma ^ { 2 } \tau _ { p } ^ { 2 } + 4 \right ) ^ { 2 } } { 4 \left ( \tau _ { p } ^ { 2 } + D _ { B } ^ { 2 } \sigma ^ { 2 } \right ) \left ( \sigma ^ { 2 } \tau _ { p } ^ { 2 } + 4 \right ) } }.$$

The temporal width of photon B in the non-heralded and heralded case can be obtained immediately from the expressions (3) and (4), respectively, by switching DA to DB and vice versa.

## Optimization of the pump laser in the asymmetric case

The conditions for the function t Ah ( t p ) to have a well-defined global minimum are very complicated in the general case. However, they can be considerably simplified if we assume that DA and DB have the same sign, which is certainly a justified assumption in realistic situations. To write them explicitly we first introduce the following notation:

$$\xi _ { i, j } = & \left [ i \frac { D _ { A } - D _ { B } + j \sqrt { ( D _ { A } - D _ { B } ) ^ { 2 } - 8 D _ { A } ( D _ { A } + D _ { B } ) } } { 2 D _ { A } D _ { B } } \right ] ^ { 1 / 2 },$$

$$\zeta _ { i, j } = & \left [ i \frac { D _ { A } - D _ { B } + j \sqrt { ( D _ { A } - D _ { B } ) ^ { 2 } - 8 D _ { B } ( D _ { A } + D _ { B } ) } } { D _ { B } ( D _ { A } + D _ { B } ) } \right ] ^ { 1 / 2 }.$$

Here we focus on the typical QC scheme with SMFs, in which case it is always DA DB , &lt; 0. Then the right-hand side of the expression (4) reaches its minimum for

$$\tau _ { p } ^ { ( - ) } = 2 \sqrt { - \frac { 2 ( D _ { A } + D _ { B } ) - \sigma ^ { 2 } D _ { B } ( D _ { A } - D _ { B } ) + \sigma ^ { 4 } D _ { A } D _ { B } ^ { 2 } } { 8 + 2 \sigma ^ { 2 } ( D _ { A } - D _ { B } ) + \sigma ^ { 4 } D _ { B } ( D _ { A } + D _ { B } ) } }$$

in the three following cases: (1) when 10 7 . DA ≈ ( 4 √ 2 + ) 5 DA &lt; DB &lt; [ ( 4 √ 2 -5 ) / 7 ] DA ≈ 0 094 . DA , (2) when DB ≤ ( 4 √ 2 + ) 5 DA and one of the inequalities s &lt; x + 1 , -1 or x + + 1 , 1 &lt; s is true, (3) when [ ( 4 √ 2 -5 ) / 7 ] DA ≤ DB and one of the inequelities s &lt; z -1 , + 1 or z - -1 , 1 &lt; s is true. If none of the above sets of conditions is fulfilled, then the function t Ah ( t p ) does not have a global minimum. If DB ≤ ( 4 √ 2 + ) 5 DA but x + 1 , -1 ≤ s ≤ x + + 1 , 1 it is monotonically increasing, meaning that the lowest temporal width of photon A is reached for t p → 0. On the other hand if [ ( 4 √ 2 -5 ) / 7 ] DA ≤ DB but z -1 , + 1 ≤ s ≤ z - -1 , 1 the function t Ah ( t p ) always decreases when t p grows. Therefore, in this situation the lowest temporal width of photon A is reached for t p → ¥ .

## Temporal widths of SPDC photons when the detection timing jitter is non-zero

When the timing jitter, t JA , is non-zero the difference between the detection time of photon A, t A , and the time of its arrival at the measurement system, t 0 A , can be described by the probability distribution function q tA ( , t 0 A , t JA ) = MA exp [ -( t A -t 0 A ) 2 / ( 2 t 2 JA )] , where MA is the normalization constant. Then, the probability distribution for the detection time of this photon in the case when the detection time of photon B is unknown can be calculated as p A ( t A ) = ∫ d t 0 A pA ( t 0 A ) q tA ( , t 0 A , t JA ) . The marginal distribution function pA ( x ) has already been defined in the text between the equations (2) and (3). It is straightforward to check that the standard deviation of p A ( t A ) is equal to t J A =( t 2 A + t 2 JA ) 1 / 2 . The above formula gives the temporal width of the non-heralded photon A for the case of non-zero jitter.

While the value of t J A depends only on the timing jitter of the detector A, the analogous temporal width of photon A found in the heralded case, t J Ah , would be influenced also by the timing jitter of the other detector, t JB . In order to calculate it, one has to take the joint probability formula for the detection of photon A at the time t A and the detection of photon B at the time t B , which can be derived from (2) as pAB ( t A , t B ) = | y DADB ( t A , t B ) | 2 , and modify it to the following form: p AB ( t A , t B ) = ∫ d t 0 A ∫ d t 0 B pAB ( t 0 A , t 0 B ) q tA ( , t 0 A , t JA ) q tB ( , t 0 B , t JB ) . In the above formula t 0 B is the arrival time of photon B to the heralding detector. The probability distribution of the detection time of photon A, conditioned on the detection of photon B at the time TB , is then given by p Ah ( t A , t B = TB ) = p AB ( t A , t B = TB ) / ∫ d t A p AB ( t A , t B = TB ) . The standard deviation of the resulting function is

$$\tau _ { A b } ^ { I } = \sqrt { \tau _ { A b } ^ { 2 } + \tau _ { J A } ^ { 2 } + X \tau _ { J B } ^ { 2 } },$$

where t Ah is the temporal width of the photon A calculated for zero jitter case and X =( t 2 p -DADB s 2 ) 2 ( s t 2 2 p -4 ) 2 / [( t 2 p + DB 2 s 2 ) 2 ( s t 2 2 p + ) ] 4 2 .

## Effective phase-matching function width for the BBO crystal

The approximate value of the effective phase-matching function width for a particular nonlinear crystal can be calculated by using the following formula: s =[( d 2 k / Wf 2 + 5 / L 2 ) / d w 2 ] 1 / 2 , where Wf is the width of transverse spatial mode collected by the SMF and L is the length of the crystal 27 . Furthermore, by d k and dw we defined the partial derivatives of the phase mismatch D kz over the transverse component of wave vector of the produced signal photons and their angular frequencies, respectively. In this work we are interested in type I SPDC process taking place in BBO crystal. It is a negative uniaxial crystal, which means that the pump photons have to be extraordinarily polarized, while the polarizations of signal and idler photons are always ordinary 41 . Assuming that the pump pulses propagate along the z direction, the phase mismatch for the investigated process is given by D kz =( w s + w i ) n e ( w s + w q i , ) / c -[ w s no ( w s ) / c -k sx 2 ] 1 / 2 -[ w i no ( w i ) / c -k 2 ix ] 1 / 2 , where w s ( w i ) is the angular frequency of the signal (idler) photon and ksx ( kix ) is its transverse wave vector component. The refractive index for pump photons depends on the angle between the Z axis and the optic axis, q , as follows 41 :

$$n ^ { c } ( \omega, \theta ) = n _ { o } ( \omega ) \sqrt { \frac { 1 + \tan ^ { 2 } \theta } { 1 + [ n _ { o } ( \omega ) / n _ { e } ( \omega ) ] ^ { 2 } \tan ^ { 2 } \theta } }.$$

The approximate formula for the dependence of the refractive index of the ordinarily [extraordinarily] polarized photons propagating in the BBO crystal, no ( w ) [ ne ( w ) ], on their angular frequency can be found in Ref. 41 . The expression for q can be obtained by solving the equation D kz = 0.

## QKD security analysis

For the BB84 protocol realized in the setup configuration illustrated in Fig. 2 the lower bound on the key generation rate is given by K = p exp [ 1 -2 H Q ( )] , where H x ( ) = -x log 2 x -( 1 -x ) log 2 ( 1 -x ) is the Shannon entropy and Q denotes the quantum bit error rate (QBER) in the raw key generated by the legitimate participants of the protocol 42 . In the above formula p exp is the probability of accepting a given event by Alice and Bob for the process of key generation. Obviously, both Q and p exp depend

on the duration time of the detection windows chosen by the participants of the protocol. For a single such window of width xt X the probability for registering a photon of temporal width t X is given by

$$\eta ( \xi ) = ( 2 \pi ) ^ { - 1 / 2 } \int _ { - \xi / 2 } ^ { \tilde { \xi } / 2 } d y \exp ( - y ^ { 2 } / 2 ) = \text{erif} ( \xi / 2 \sqrt { 2 } ).$$

On the other hand, the probability for registering a dark count in one of the two single-photon detectors can be calculated as

$$P _ { X } ( \xi ) & = 2 d \xi \tau _ { X },$$

where d is the dark count rate for a given single-photon detector. For the calculations performed in this work we assume that d = 1kHz both in the case of Alice's and Bob's detectors. Also, for every pair of the investigated values of LA and LB we separately optimize the parameters x A and x B in order to get the best possible outcomes for the legitimate parties. We consider the situation in which the dark counts are the only source of errors in the raw key. Since narrowing the detection windows reduces all the possible errors that are uncorrelated with the real signals in exactly the same way, adding such errors to the model can be easily made just by appropriate increase of d . On the other hand, the errors that are connected to the real signals, e.g. polarization rotation, would only slightly change the obtained results and not in qualitative way.

We consider the case of perfect SPDC source, always emitting a single pair of photons when the pump pulse propagates through the crystal. With this assumption the probability p exp for the scheme illustrated in Fig. 2 can be approximated by

$$P _ { \exp } & \quad \approx \quad T _ { A } \eta ( \xi _ { A } ) T _ { B } \eta ( \xi _ { B } ) + T _ { A } \eta ( \xi _ { A } ) [ 1 - T _ { B } \eta ( \xi _ { B } ) ] P _ { B h } ( \xi _ { B } ) + [ 1 - T _ { A } \eta ( \xi _ { A } ) ] T _ { B } \eta ( \xi _ { B } ) P _ { A h } ( \xi _ { A } ) + \\ & \quad + \quad [ 1 - T _ { A } \eta ( \xi _ { A } ) ] [ 1 - T _ { B } \eta ( \xi _ { B } ) ] P _ { A } ( \xi _ { A } ) P _ { B h } ( \xi _ { B } ),$$

where TA ( TB ) is the transmittance of quantum channel connecting the SPDC source with Alice (Bob), given by TA = 10 -a A L A / 10 ( TB = 10 -a BLB / 10 ). The probabilities for a dark count to be registered by Alice or Bob in a particular detection window, denoted by PAh and PBh respectively, can be calculated by inserting the expression (4), and the analogous expression for the temporal width of the heralded photon B, into the formula (11). On the other hand, in order to obtain PA one should use the equation (3) instead of (4). This probability, appearing only in the last term on the right-hand side of Eq. (13), is needed to properly account for the case when neither of the signal photons is detected by the measurement systems of Alice and Bob. In this situation one of the dark counts registered by them has to be treated as a 'heralding' click, while the other one is 'heralded' by it (obviously, the exact choice does not matter here, as can be confirmed by checking that t A t Bh ≡ t B t Ah ). On the other hand, the second and third terms on the right-hand side of Eq.(13) correspond to the case when only one of the two photons produced by the source causes a click in one of the measurement systems, but the event is still accepted for the key generation process due to a dark count registered in the other detection system in the appropriately narrowed time window. Finally, the first term on the right-hand side of Eq.(13) accounts for the desired situation in which both SPDC photons from a given pair are detected by Alice's and Bob's measurement systems.

In the case of the simplified QKD scheme considered in our work an error in the raw key can be generated only if at least one of the signal photons from a given SPDC pair is lost, but the event is still accepted for key generation. Since dark counts occur in the detectors of Alice and Bob totally randomly, in all of such situations the error probability is 50%. Therefore, QBER can be calculated using the expression

$$Q = \frac { p _ { \exp } - T _ { A } \eta ( \xi _ { A } ) T _ { B } \eta ( \xi _ { B } ) } { 2 p _ { \exp } }.$$

## References

- 1. Bennett, C. H. &amp; Brassard, G. Quantum cryptography: Public key distribution and coin tossing. In Proceedings of the IEEE International Conference on Computers, Systems, and Signal Processing, Bangalore, India , vol. 11, 175-179 (IEEE, New York, 1984).
- 2. Ekert, A. K. Quantum cryptography based on bell theorem. Phys. Rev. Lett. 67 , 661-663 (1991).
- 3. Hillery, M., Bužek, V. &amp; Berthiaume, A. Quantum secret sharing. Phys. Rev. A 59 , 1829-1834 (1999).
- 4. Bennett, C. H. et al. Teleporting an unknown quantum state via dual classical and einstein-podolsky-rosen channels. Phys. Rev. Lett. 70 , 1895 (1993).
- 5. Bennett, C. H. &amp; Wiesner, S. J. Communication via one- and two-particle operators on einstein-podolsky-rosen states. Phys. Rev. Lett. 69 , 2881 (1992).

- 6. Louisell, W. H., Yariv, A. &amp; Siegman, A. E. Quantum fluctuations and noise in parametric processes. i. Phys. Rev. 124 , 1646-1654 (1961).
- 7. Burnham, D. C. &amp; Weinberg, D. L. Observation of simultaneity in parametric production of optical photon pairs. Phys. Rev. Lett. 25 , 84-87 (1970).
- 8. Fasel, S. et al. High-quality asynchronous heralded single-photon source at telecom wavelength. New J. Phys. 6 , 163 (2004).
- 9. Bock, M., Lenhard, A., Chunnilall, C. &amp; Becher, C. Highly efficient heralded single-photon source for telecom wavelengths based on a ppln waveguide. Opt. Express 24 , 23992-24001 (2016).
- 10. Pomarico, E., Sanguinetti, B., Guerreiro, T., Thew, R. &amp; Zbinden, H. MHz rate and efficient synchronous heralding of single photons at telecom wavelengths. Opt. Express 20 , 23846 (2012).
- 11. Ramelow, S. et al. Highly efficient heralding of entangled single photons. Opt. Express 21 , 6707 (2013).
- 12. Kaneda, F., Garay-Palmett, K., U'Ren, A. B. &amp; Kwiat, P. G. Heralded single-photon source utilizing highly nondegenerate, spectrally factorable spontaneous parametric downconversion. Opt. Express 24 , 10733-10747 (2016).
- 13. Mattle, K., Weinfurter, H., Kwiat, P. G. &amp; Zeilinger, A. Dense coding in experimental quantum communication. Phys. Rev. Lett. 76 , 4656-4659 (1996).
- 14. Bouwmeester, D. et al. Experimental quantum teleportation. Nat. (London) 390 , 575-579 (1997).
- 15. Pan, J.-W., Bouwmeester, D., Weinfurter, H. &amp; Zeilinger, A. Experimental entanglement swapping: Entangling photons that never interacted. Phys. Rev. Lett. 80 , 3891-2894 (1998).
- 16. Jennewein, T., Simon, C., Weihs, G., Weinfurter, H. &amp; Zeilinger, A. Quantum cryptography with entangled photons. Phys. Rev. Lett. 84 , 4729-4732 (2000).
- 17. Tittel, W., Zbinden, H. &amp; Gisin, N. Experimental demonstration of quantum secret sharing. Phys. Rev. A 63 , 042301 (2001).
- 18. Sedziak, K., Lasota, M. &amp; Kolenderski, P. Reducing detection noise of a photon pair in a dispersive medium by controlling its spectral entanglement. Optica 4 , 84 (2017).
- 19. Lasota, M. &amp; Kolenderski, P. Quantum communication improved by spectral entanglement and supplementary chromatic dispersion. Phys. Rev. A 98 , 062310 (2018).
- 20. Sedziak-Kacprowicz, K., Lasota, M. &amp; Kolenderski, P. Remote temporal wavepacket narrowing. Sci. Rep. 9 , 3111 (2019).
- 21. The recalculation can be easily done e.g. using the script placed on the website www.rpphotonics.com/group\_velocity\_dispersion.html.
- 22. See e.g. Corning SMF-28e+ Optical Fiber specifications. www.corning.com/media/worldwide/coc/documents/ Fiber/PI1463\_07-14\_English.pdf.
- 23. Shcheslavskiy, V. et al. Ultrafast time measurements by time-correlated single photon counting coupled with superconducting single photon detector. Rev. Sci. Instrum. 87 , 053117 (2016).
- 24. Yan, Z. et al. An ultra low noise telecom wavelength free running single photon detector using negative feedback avalanche diode. Rev. Sci. Instrum. 83 , 073105 (2012).
- 25. Tosi, A. et al. Ingaas/inp single-photon avalanche diodes show low dark counts and require moderate cooling. In Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series , vol. 7222, 72221G (Proc. of SPIE, 2009).
- 26. Divochiy, A. et al. Single photon detection system for visible and infrared spectrum range. Opt. Lett. 43 , 6085-6088 (2018).
- 27. Gajewski, A. &amp; Kolenderski, P. Spectral correlation control in down-converted photon pairs. Phys. Rev. A 94 , 013838 (2016).
- 28. Kolenderski, P., Wasilewski, W. &amp; Banaszek, K. Modelling and optimization of photon pair sources based on spontaneous parametric down-conversion. Phys. Rev. A 80 , 013811 (2009).
- 29. Sun, P., Guo, J., Wang, T., Liu, L. &amp; Feng, J. Method for measuring the intensity distribution of a small beam spot with phase retrieval. Opt. Laser Eng. 57 , 109 - 113 (2014).
- 30. Hwang, W.-Y. Quantum key distribution with high loss: Toward global secure communication. Phys. Rev. Lett. 91 , 057901 (2003).

- 31. Korzh, B. et al. Provably secure and practical quantum key distribution over 307 km of optical fibre. Nat. Photonics 9 , 163-168 (2015).
- 32. Yin, H.-L. et al. Measurement-device-independent quantum key distribution over a 404 km optical fiber. Opt. Express 117 , 190501 (2016).
- 33. Boaron, A. et al. Secure quantum key distribution over 421 km of optical fiber. Phys. Rev. Lett. 121 , 190502 (2018).
- 34. Wang, Q., Wang, X.-B. &amp; Guo, G.-C. Practical decoy-state method in quantum key distribution with a heralded singlephoton source. Phys. Rev. A 75 , 012312 (2007).
- 35. Wang, Q. &amp; Wang, X.-B. Efficient implementation of the decoy-state measurement-device-independent quantum key distribution with heralded single-photon sources. Phys. Rev. A 88 , 052332 (2013).
- 36. Zhu, J.-R., Li, J., Zhang, C.-M. &amp; Wang, Q. Parameter optimization in biased decoy-state quantum key distribution with both source errors and statistical fluctuations. Quantum Inf. Process. 16 , 238 (2017).
- 37. Eraerds, P., Walenta, N., Legre, M., Gisin, N. &amp; Zbinden, H. Quantum key distribution and 1 gbit/s data encryption over a single fibre. New J. Phys. 12 , 063027 (2010).
- 38. Elkouss, D., Martínez-Mateo, J., Ciurana, A. &amp; Martín, V. Secure optical networks based on quantum key distribution and weakly trusted repeaters. IEEE J. Opt. Commun. Netw. 5 , 316-328 (2013).
- 39. Ciurana, A. et al. Quantum metropolitan optical network based on wavelength division multiplexing. Opt. Express 22 , 1576-1593 (2014).
- 40. Lutz, T., Kolenderski, P. &amp; Jennewein, T. Demonstration of spectral correlation control in a source of polarization entangled photon pairs at telecom wavelength. Opt. Lett. 39 , 1481 (2014).
- 41. Dmitriev, V. G., Gurzadyan, G. G. &amp; Nikogosyan, D. N. Handbook of Nonlinear Optical Crystals (Springer, 1999), 3rd edn.
- 42. Kraus, B., Gisin, N. &amp; Renner, R. Lower and upper bounds on the secret-key rate for quantum key distribution protocols using one-way classical communication. Phys. Rev. Lett. 95 , 080501 (2005).

## Acknowledgements

The authors acknowledge financial support by the Foundation for Polish Science (FNP) (project First Team co-financed by the European Union under the European Regional Development Fund), Ministy of Science and higher Education, Poland (MNiSW) (grant no. 6576/IA/SP/2016) and National Science Centre, Poland (NCN) (Sonata 12 grant no. 2016/23/D/ST2/02064). We wish to thank National Laboratory of Atomic, Molecular and Optical Physics, Torun, Poland for the support.

## Author contributions statement

M.L. performed most of the analytical and numerical calculations. P.K. supervised the research. Both authors contributed to writing the paper.

## Additional information

Competing interests: The authors declare no competing interests.
Miko\{l}aj Lasota Piotr Kolenderski
Optimal photon pairs for quantum communication protocols
We theoretically investigate the problem of finding optimal characteristics
of photon pairs, produced in the spontaneous parametric down-conversion (SPDC)
process, for fiber-based quantum communication (QC) protocols. By using the
accessible setup parameters, the pump pulse duration and the extended
phase-matching function width, we minimize the temporal width of SPDC photons
within the general scenario. This allows one to perform more effectively the
temporal filtering procedure, which aims at reducing the noise acquired by the
measurement devices. Moreover, we compare the obtained results with the
achievable parameter values for SPDC sources based on $\beta$-Barium Borate
(BBO) crystal. We also investigate the influence of non-zero detection timing
jitter. Finally, we apply our optimization strategy to a simple quantum key
distribution scheme to show that the full optimization of an SPDC source can
potentially extend the maximal security distance by several tens of kilometres,
which is around 30\% more as compared to previous approaches.