## Numerical modeling of SNSPD absorption utilizing optical conductivity with quantum corrections

M. Bar´ anek 1 , P. Neilinger 1 2 , , S. Kern , and M. Grajcar 1 1 2 ,

- 1 Department of Experimental Physics, Comenius University, SK-842 48, Bratislava, Slovakia
- 2 Institute of Physics, Slovak Academy of Sciences, D´bravsk´ cesta 9, SK-845 11, Bratislava, Slovakia u a

ABSTRACT. Superconducting nanowire single-photon detectors are widely used in various fields of physics and technology, due to their high efficiency and timing precision. Although, in principle, their detection mechanism offers broadband operation, their wavelength range has to be optimized by the optical cavity parameters for a specific task. We present a study of the optical absorption of a superconducting nanowire single photon detector (SNSPD) with an optical cavity. The optical properties of the niobium nitride films, measured by spectroscopic ellipsometry, were modelled using the Drude-Lorentz model with quantum corrections. The numerical simulations of the optical response of the detectors show that the wavelength range of the detector is not solely determined by its geometry, but the optical conductivity of the disordered thin metallic films contributes considerably. This contribution can be conveniently expressed by the ratio of imaginary and real parts of the optical conductivity. This knowledge can be utilized in detector design.

## INTRODUCTION

Strengthening the security of digital communication channels by building Quantum Key Distribution (QKD) infrastructure is an ongoing effort worldwide[1, 2]. Most QKD systems are based on entangled photon pairs as quantum bits and utilize the existing, low-loss fiber optical infrastructure [3]. To implement an efficient algorithm, a single-photon detector with a high overall detection efficiency, low dark counts, and high speed is required[4]. The State-of-the-art detectors that fulfill these requirements are the superconducting nanowire single-photon detectors (SNSPD)[5]. They are superior in terms of detection efficiency and dark count rates to other types[6], such as the Avalanche Photodiodes (APD)[7], and do not require sub-kelvin temperatures, as the Transition Edge Sensors (TES)[8] do. Moreover, their superior properties make them beneficial in a broad range of applications in different areas of physics, for example, physical chemistry and spectroscopy[9], fluorescent luminescence[10], or fast space-to-ground communication[11] often requiring a lower wavelength range. The detector's wavelength can be optimized by its design depending on the required spectral range. This is also important in the implementation of modern quantum networks[3, 12] as not just the detection efficiency of the detector, but its bandwidth is a relevant parameter.

The working principle of SNSPD is straightforward: a nanowire in the superconducting state is currentbiased close to the edge of the normal metal transition. The absorption of a single photon results in the destruction of superconductivity at some place along the nanowire. Due to the current bias, the normal state region expands, which results in a voltage pulse across the nanowire. These nanowires are fabricated on 615 nm thick, highly disordered superconductors, such as Niobium Nitride[13], Titanium Niobium Nitride[14], Tungsten Silicide[15] and their width is in the range of 30-200 nm.

The overall system detection efficiency (SDE) of

Fig. 1. Scheme of the SNSPD with coupled optical fiber. Nanowire has thickness t , width w , and spacing s . Optical λ /4 resonator consists of a dielectric spacer with thickness d 1 and a gold reflector with d 2 . The gap between optical fiber and nanowire is g . Right: Top-down view of the nanowire of diameter 2 r , with red line showing the position of the cross-sectional view.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-682ed9ae-2ea0-432e-a277-833e377fab06.jpg)

detectors is determined by the efficiency of the photon coupling to the active area of the detector η coupling , the photon absorption efficiency of the superconductor η abs , and by the above-described efficiency of transformation of the absorbed photon to electric signal η intrinsic and is given as[16]:

$$S D E = \eta _ { \text{coupling} } \cdot \eta _ { \text{abs} } \cdot \eta _ { \text{intrinsic} }, \quad \ \ ( 1 )$$

The absorption of a nanowire (thin film) η abs on bare substrate is insufficient (below 30%, Fig. 2). To increase the absorption at the desired wavelength ( λ =1550 nm) the nanowires are commonly fabricated in an optical resonator, at the cost of limiting its bandwidth. A common choice is a λ/ 4 resonant cavity [17].

This resonator consists of a thin dielectric layer with thickness d 1 , matching its optical length Λ with λ/ 4 of the desired light's wavelength. The nanowire is fabricated on one side of the dielectric and the other side is covered with a low-loss mirror metallic layer (as shown in Fig. 1). The resonator enhances the absorption of the photons in the nanowire close to 100% (Fig. 2) by creating an antinode of the electric field of TE mode in the superconductor. The detectable light is

guided by an optical fiber to this structure.

To properly model the optical properties of this structure, and thus to maximize the absorption η abs in the nanowire, not only has the dielectric thickness d 1 to be optimized, but the precise optical properties of the disordered metal have to be considered. However, their optical properties in the metallic state are not trivial due to the presence of quantum corrections, which the standard Drude model fails to describe [18]. This means that, to model their optical absorption, the optical properties of the specific thin film have to be either measured in the infrared spectra, or the right optical model - or an extrapolation method - has to be chosen, if measurements in the required range are not accessible.

In this paper, we investigate the dependence of the optical absorption of NbN nanowires in λ/ 4 resonant cavity on film thickness and nanowire design. The optical conductivity of a set of NbN films was determined from spectroscopic ellipsometry in the visible range, and they were fitted by the modified Drude-Lorentz model, which takes the quantum corrections into account [19] and provides an excellent fit of these conductivities. The presence of quantum corrections in the IR range results in a strong wavelength and thickness dependence of the optical properties of NbN films. The obtained optical fits and thickness dependencies of its parameters are used to model the absorption spectra of nanowires. The presented approach can be applied to other disordered films, such as MoC and NbTiN. The optical properties of the latter[20] are almost identical to the properties of NbN.





Fig. 2. Nanowire optical absorption η abs on bare sapphire substrate and on λ/ 4 resonator.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-fe1068f7-bc30-40dd-a63e-46f8094a6a66.jpg)

## Optical properties of NbN films

Niobium nitride samples were deposited by means of Pulsed Laser Deposition[21] from a 99 % pure niobium target on top of single-side polished sapphire substrate in a nitrogen atmosphere with 1 % of hydrogen. The film thickness varied from 8 up to 22 nm. The optical conductivities of our films were determined from spectroscopic ellipsometry measured in wavelength range from 400 to 1000 nm at room temperature.

In Ref.[19] the quantum corrections to the optical conductivity of these NbN films were thoroughly studied. These corrections are arising from localization and interaction effects both having a square root energy dependence[22]. The optical conductivity ˜ = σ σ 1 + iσ 2 is described by the modified Drude-Lorentz model in the form:

$$\intertext { w - 1 } _ { \begin{array} { c } \ s o n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \begin{array} { c } \ s u n d \end{array} } \ s u n d \end{array}$$

$$\sigma _ { 2 } ( \omega ) = \mathcal { H } [ \sigma _ { 1 } ( \omega ) ] - ( \varepsilon _ { \infty } - 1 ) \varepsilon _ { 0 } \omega \quad \ \ ( 3 )$$

respectively. Here, σ D is the Drude conductivity, Γ is the electron relaxation rate and Q is the strength of the quantum corrections, also referred to as quantumness. The second term in eq. 2 describes the transition peak at Ω L , with strength σ L , and width Γ L . This formula describes the anomalous suppression in the IR range and below, which is characteristic of these films. The optical conductivity for films with different thicknesses is shown in Fig.3 The same result can be represented by σ 1 vs. σ 2 plots in Fig.4, where the wavelength dependence is color-coded. The solid spirals are the fits with eqs. 2,3. The fitted parameters are listed in Tab. 1.

As can be seen in Fig. 3a, the Drude conductivity σ 0 varies weakly with thickness, but the relaxation rate Γ and the quantumness Q are strongly thickness dependent. This reflects the commonly observed suppression of DC conductivity and, most importantly, the significant changes in the optical properties of films, especially in the IR range. By lowering the thickness, the real part of the IR range conductivity is suppressed and the imaginary part changes its sign. The dielectric function

$$\varepsilon _ { r } ( \omega ) = \varepsilon _ { 1 } + i \varepsilon _ { 2 } = 1 + i \tilde { \sigma } ( \omega ) / ( \epsilon _ { 0 } \omega ) \quad ( 4 )$$

for various film thicknesses is plotted in Fig.3b. The known thickness dependence of these parameters allows us to conveniently model the optical properties of SNSPDs. The complex refraction index can be expressed in terms of permittivity, as:

$$\tilde { n } = \sqrt { \varepsilon _ { r } ( \omega ) }, & & ( 5 ) \\ \text{ant to note that not in} \, \text{the magnitude}$$

It is important to note, that not just the magnitude of the conductivity changes with the thickness, but the ratio of real to imaginary parts changes as well. This is easily visible by the change of the angle of the arrow pointing to the conductivity at λ 0 =1550 nm in Fig. 4. This angle is defined by the tangent of real and imaginary parts of the optical conductivity as tan θ 0 = σ 2 ( λ 0 ) /σ 1 ( λ 0 ) and is a characteristic of the films for a given thickness - it does not depend on the fill factor of the nanowire, since scaling the conductivity scales both real and imaginary part.

## Optical absorption simulations

|   d [nm] |   σ D [ MS/ m] |   /planckover2pi1 Γ[eV] |    Q |   σ 1 [ MS/ m] |   /planckover2pi1 Γ 1 [eV] |   /planckover2pi1 Ω 1 [eV] |
|----------|----------------|-------------------------|------|----------------|----------------------------|----------------------------|
|        8 |           0.77 |                    2.1  | 0.7  |           0.9  |                       1    |                       6.13 |
|        9 |           0.84 |                    2    | 0.6  |           0.95 |                       1.1  |                       6.14 |
|       14 |           1    |                    1.76 | 0.66 |           1.12 |                       0.56 |                       5.19 |
|       22 |           0.94 |                    1.77 | 0.5  |           1.07 |                       0.68 |                       5.54 |

TAB. 1. Parameters of optical model eq.(2) obtained from the ellipsometric data fit.

Fig. 3. Specific conductivity (top) and permittivity (bottom) of thin NbN films, consisting of real (solid line) and imaginary (dashed line) parts. Thick lines are the result of spectroscopic ellipsometry measurement, and thin lines result from the model described above. Dots at zero frequency are the room temperature DC conductivities.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-7e624c4d-b464-40cf-b3ce-9e3c7c0ba775.jpg)

Fig. 4. Wavelength-dependent real and imaginary part of conductivity for films with various thicknesses. The arrows show the mismatch in the simple rescaling of the conductivities, easily visible by the change in the angle θ 0 of arrows, which point at the conductivity at 1550 nm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-1fac224c-8c5e-42fc-8c6f-0158756e2eb4.jpg)

We studied the effect of the film's conductivity on the absorption spectra of SNSPD with a similar detector design to [14]. The numerical simulations were carried out in a commercial frequency-domain EM solver[23]. The air gap between the fiber and the surface of the detector leads to additional Fabry-Per´t reso onances superposing the absorption spectra of the λ/ 4 resonator. This could be expressed in the wavelength dependence of the coupling parameter η coupling ( λ ). In our 2D simulations, the optical fiber output is modelled as a plane wave coupled to air and we assume infinite size of the detector in the plane perpendicular to the light propagation, which corresponds to the ideal coupling η coupling = 1 and allows us to focus on the absorption in the nanowire η abs .

The optical absorption of the SNSPD can be described by the effective impedance model[24], which approximates the thin superconducting meander as a lumped element effective homogeneous medium, with effective sheet conductance ˜ = G G 1 + iG 2 :

$$\tilde { G } = \tilde { \sigma } \frac { w } { w + s } t \ [ S / \Box ],$$

where ˜( ) σ t is the optical conductivity (in units of S/m) of the NbN film with thickness t , and the ratio f = w/ w ( + ) is the fill s factor. This approximation is valid for the width and the thickness of the nanowires w, t /lessmuch λ . As these parameters are usually w ∼ 30-200 nm, t ∼ 5-22 nm, and s ∼ 50-300 nm, this approximation holds well for λ &gt; 1 µ m. The validity of this approximation is demonstrated by numerical modeling of the absorption spectra of a detector

with different nanowire widths and constant fill factor f = 0 5, shown in Fig. 5. The absorption spectra . are almost identical for w /lessorsimilar 400nm. The negligible differences in the absorption maximum ( &lt; 0 5%) and . its wavelength ( &lt; 50 nm) may originate in simulation mesh discretization. The ≈ 1% losses originate in the Au layer with d 2 = 100 nm.

Fig. 5. Simulated absorption spectra of the detector for different nanowire widths and constant fill factor f = 0 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-c1f1f6f2-c4c5-4c0b-a84b-80687884602f.jpg)

In Fig. 6, we present the absorption spectra of a detector with 100 nm width NbN nanowire of different thicknesses t = 8 -22 nm on SiO 2 layer with thickness d 1 = 250 nm. The absorption maximum corresponding to the optical length of an ideal λ/ 4 resonator is 4Λ = 4 nd 1 ≈ 1600 nm. For each thickness, the fill factor f was tuned to maximize the absorption at λ 0 = 1550 nm. The maximal absorptions at λ 0 , and the wavelength of the absolute maximal absorption λ are listed in Table 2. Note, that for this analysis, the optical properties of NbN films were extracted from ellipsometric measurements on films deposited on c-cut sapphire substrate. The optical properties of films deposited on different substrates, for example, SiO 2 or Si 3 N , can differ. 4 However, the presented results can be easily generalized for any combination of film and dielectric material, as we will show below.





Fig. 6. Simulated absorption spectra of nanowires with different thicknesses for fill factor maximizing absorption at 1550 nm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-472f000c-0517-44f1-b7d0-e9a5db8fac72.jpg)

The absorption spectra for t = 8, 9, 14, 22 nm are shown in Fig. 6. As it is visible, due to the thickness dependence of the optical conductivity (refraction index), the thickness of the film significantly shifts the position of the absorption maxima and changes the shape of the absorption spectra. The maximum of the absorption value max( η abs ), at its peak wavelength λ max , is also

TAB. 2. Simulated maximum absorption at 1550 nm and the optimized fill factor for 1550 nm, wavelength of maximum absorption, and maximum absorption value for different thickness

| 250 nm SiO 2 , 4Λ = 1600 nm   | 250 nm SiO 2 , 4Λ = 1600 nm   | 250 nm SiO 2 , 4Λ = 1600 nm   | 250 nm SiO 2 , 4Λ = 1600 nm   | 250 nm SiO 2 , 4Λ = 1600 nm   | 250 nm SiO 2 , 4Λ = 1600 nm   |
|-------------------------------|-------------------------------|-------------------------------|-------------------------------|-------------------------------|-------------------------------|
| d [nm]                        | σ 2 /σ 1 @1550 nm             | max( η abs ) @1550 nm [%]     | fill fac- tor                 | λ max [nm]                    | max ( η abs ) [%]             |
| 8                             | 0.025                         | 99.1                          | 0.6                           | 1650                          | 99.3                          |
| 9                             | 0.092                         | 99.2                          | 0.45                          | 1550                          | 99.2                          |
| 14                            | 0.242                         | 96.8                          | 0.3                           | 1350                          | 99.1                          |
| 22                            | 0.314                         | 91.5                          | 0.225                         | 1250                          | 97.9                          |

Fig. 7. Optical absorption in nanowire with fixed G 2 = 0 mS (left) and varying G 1 , and with fixed real effective conductance G 1 = 2 5 . mS (right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-56fe811a-bc4a-411f-9b34-d0302c8fa11d.jpg)

lowered for thicker films due to the increased reflection, but it is negligible compared to the drop of absorption at the desired wavelength λ 0 = 1550 nm. To study how G 1 and G 2 separately affect the absorption spectra, we simulated the absorption of the nanowire with w = 100 nm, s = 100 nm ( f = 0 5), . t = 10 nm on d 1 = 250 nm SiO 2 dielectric and d 2 = 100 nm Au reflector.

The absorption spectra for varying G 1 and constant G 2 = 0 reveal that the real part of the effective sheet conductance G 1 ( σ 1 , respectively) mainly affects the value of the absorption maxima, but not their position λ max , shown in Fig.7a. On the other hand, variation of G 2 ( σ 2 , respectively) for constant G 1 results in the shift of absorption peak λ max , shown in Fig.7b. The simulated range of tan θ is approx. 0 - 0.4, which is comparable to the values of our films (0 - 0.32). Thus, the resonant absorption peak of the detector is shifted by the imaginary part of the conductance towards lower values and thus, its wavelength λ max is not solely determined by the thickness of the dielectric layer, as one could simply assume. This effect, governed by tan θ , is universal and should be present in any thin film with similar values of tan θ .

The dependence of 1 /λ max on the conductivity of the NbN films, represented by tan θ , is linear for a broad range of tan θ , as is visible in Fig. 8, and can be approximated as:

$$\frac { 1 } { \lambda _ { m a x } ( \sigma _ { 2 } ) } = \frac { 1 } { \lambda _ { m a x } ( \sigma _ { 2 } = 0 ) } + \frac { 1 } { \lambda _ { c } } \tan \theta, \quad ( 7 )$$

where λ c is a fitting parameter. The value λ c is governed by the optical properties of the dielectric layer and its thickness. For the studied geometry of the detector on SiO , 2 λ c ∼ 4 2 . µ m. Here, we assumed tan θ to be constant and equal to tan θ 0 , which is a feasible approximation in the vicinity of λ 0 , as is shown in Fig. 10.

Since tan θ increases with decreasing wavelength, the positive feedback further lowers the resonance wavelength.

Fig. 8. Dependance of wavelengths of the simulated absorption maxima on the conductivity ratio tan θ = σ /σ 2 1 = G /G 2 1 , given by samples with different thicknesses. Left: On top of the 250nm SiO 2 substrate, Right: On top of the 200 nm Si N 3 4 membrane

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-0e825130-1848-4642-bc2d-2612c119c2e3.jpg)

This simple relation can be directly implemented in the detector design. The absorption maximum shift due to the imaginary part of the conductivity (contribution of the thin film), for a given dielectric layer, limits the possible range of film thicknesses.

However, using e.g. a commercially available 200 nm thick Si N 3 4 membrane, the resonance wavelength determined by the membrane is 4Λ ≈ 1800 nm and the corresponding absorption at λ 0 is thus lower, as can be seen in Fig. 8. Using thicker NbN films, the resonance wavelength can be lowered and the absorption can be significantly increased, see Tab.3. For each film thickness, the fill factor was optimized to obtain the maximum absorption at λ 0 . Thus, fine-tuning the absorption peak position up to 200 nm is possible, without significant loss of total optical absorption η abs .

TAB. 3. Table of sample parameters with simulated maximum absorption at 1550 nm, optimized fill factor for 1550 nm, the wavelength of maximum absorption, and maximum absorption value, respectively.

| 200 nm Si 3 N 4 , 4Λ = 1800 nm   | 200 nm Si 3 N 4 , 4Λ = 1800 nm   | 200 nm Si 3 N 4 , 4Λ = 1800 nm   | 200 nm Si 3 N 4 , 4Λ = 1800 nm   | 200 nm Si 3 N 4 , 4Λ = 1800 nm   | 200 nm Si 3 N 4 , 4Λ = 1800 nm   |
|----------------------------------|----------------------------------|----------------------------------|----------------------------------|----------------------------------|----------------------------------|
| d [nm]                           | σ 2 /σ 1 @1550 nm                | max( η abs ) @1550 nm [%]        | fill fac- tor                    | λ max [nm]                       | max ( η abs ) [%]                |
| 8                                | 0.025                            | 93.3                             | 0.65                             | 1800                             | 98.8                             |
| 9                                | 0.092                            | 95.3                             | 0.5                              | 1750                             | 98.7                             |
| 14                               | 0.242                            | 98.6                             | 0.3                              | 1600                             | 98.6                             |
| 22                               | 0.314                            | 97.5                             | 0.225                            | 1450                             | 98.1                             |

The fit parameter for the 200 nm Si N 3 4 membrane is λ c = 6 3 . µ m, estimated from the linear fit in Fig. 8.

## Discussion

In similar studies, the properties of a single deposited film are often measured (for example by spectroscopic ellipsometry), and the determined specific conductivity is considered to be thickness-independent in the simulations. Consequently, the nanowire's thickness and width are then tuned to maximize the optical absorption in the nanowire at the aimed wavelength λ 0 . For example, in Ref. [25], or Ref. [26] the same refraction index was used for samples with thicknesses ranging from 4 to 10 nm, in case of Ref. [13] even from 1 to 16 nm. Rescaling the thickness of a specific film to a requested value results in the inaccuracy of the simulated absorption. This fact is acknowledged in [24], however, no further analysis is present. Rarely, the refractive index is also considered to be wavelength independent[25].

To illustrate the error arising from this approximation, we simulated the optical absorption in the above-described nanowire for the respective conductivity ˜( ), and the corresponding optical indices, given by σ t eq. 5. For example, using the conductivity of the film with t meas =14 nm and simply rescaling it to t = 9 nm for the simulation would result in a decrease of the simulated absorption maxima by approximately 200 nm, as can be seen in Fig.9.

Fig. 9. Simulated optical absorption of 9 nm thick sample on top of the resonant cavity. Other curves show the change in optical absorption of the optical conductivity (refraction index) taken from a sample with different thicknesses.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-b1bd6b21-9cf4-4065-90ee-4ca962c5f27e.jpg)

If one assumes the optical conductivity (refraction index) of these films to be wavelength-independent, the error in the absorption spectra will accumulate with the distance from the reference wavelength, as is shown in Fig. 10.

Fig. 10. Simulated absorption of a sample with wavelength dependent n λ ( ), constant n λ ( = 1000 1300 1550nm) and , , n corresponding to constant σ λ ( = 1550nm), respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00623v1-figure-11cf2060-51d0-4b9f-a609-0c27a1a16596.jpg)

## CONCLUSION

We modelled the absorption spectra of NbN nanowires in a λ/ 4 resonator for different film thicknesses and showed that the optical conductivity of thin films significantly affects the absorption spectra. The optical conductivity of disordered metallic films is strongly affected by the presence of quantum corrections, resulting in nontrivial wavelength and thickness dependence of the optical conductivity. The real part of the conductivity determines the amplitude of absorption of the nanowire in the λ/ 4 optical resonator, whereas the imaginary part contributes to the wavelength shift of the maxima. The wavelength shift of the absorption maxima can be expressed in terms of the ratio of imaginary and real parts of the optical conductivity, which is a characteristic of a film with a particular thickness. This knowledge can be utilized in detector design and optimization of its optical properties. Thus, to properly simulate and optimize the absorption of the detectors at a required wavelength, it is necessary to work with the precise optical conductivity of the particular layers. These can be obtained by direct measurements of the samples in the required range, or by utilizing a model that describes the optical conductivity of the disordered samples well, as presented in Ref. [19].

## ACKNOWLEDGMENT

This work was supported by the Slovak Research and Development Agency under the contract APVV-200425, and by the project skQCI (101091548), founded by the European Union (DIGITAL) and the Recovery and Resilience Plan of the Slovak Republic.

## REFERENCES

- 1. Yu-Ao Chen, Qiang Zhang, Teng-Yun Chen, WenQi Cai, Sheng-Kai Liao, Jun Zhang, Kai Chen, Juan Yin, Ji-Gang Ren, Zhu Chen, Sheng-Long Han, Qing Yu, Ken Liang, Fei Zhou, Xiao Yuan, Mei-Sheng Zhao, Tian-Yin Wang, Xiao Jiang, Liang Zhang, Wei-Yue Liu, Yang Li, Qi Shen, Yuan Cao, Chao-Yang Lu, Rong Shu, Jian-Yu Wang, Li Li, Nai-Le Liu, Feihu Xu, Xiang-Bin Wang, Cheng-Zhi Peng, and Jian-Wei Pan. An integrated space-to-ground quantum communication network over 4,600 kilometres. Nature , 589(7841):214-219, Jan 2021.
- 2. Domenico Ribezzo, Mujtaba Zahidy, Ilaria Vagniluca, Nicola Biagi, Saverio Francesconi, Tommaso Occhipinti, Leif K. Oxenløwe, Martin Lonˇ cari´, c Ivan Cviti´, c Mario Stipˇevi´, c c ˇ iga Z Puˇ savec, Rainer Kaltenbaek, Anton Ramˇ sak, Francesco Cesa, Giorgio Giorgetti, Francesco Scazza, Angelo Bassi, Paolo De Natale, Francesco Saverio Cataliotti, Massimo Inguscio, Davide Bacco, and Alessandro Zavatta. Deploying an inter-european quantum network. Advanced Quantum Technologies , 6(2):2200061, 2023.
- 3. Djeylan Aktas, Bruno Fedrici, Florian Kaiser, Tommaso Lunghi, Laurent Labont´ e, and S´bastien e Tanzilli. Entanglement distribution over 150 km in wavelength division multiplexed channels for quantum cryptography. Laser &amp; Photonics Reviews , 10(3):451-457, 2016.
- 4. Hiroyuki Shibata, Toshimori Honjo, and Kaoru Shimizu. Quantum key distribution over a 72db channel loss using ultralow dark count superconducting single-photon detectors. Opt. Lett. , 39(17):5078-5081, Sep 2014.
- 5. G. N. Gol'tsman, O. Okunev, G. Chulkova, A. Lipatov, A. Semenov, K. Smirnov, B. Voronov, A. Dzardanov, C. Williams, and Roman Sobolewski. Picosecond superconducting singlephoton optical detector. Applied Physics Letters , 79(6):705-707, 08 2001.
- 6. Iman Esmaeil Zadeh, J. Chang, Johannes W. N. Los, Samuel Gyger, Ali W. Elshaari, Stephan Steinhauer, Sander N. Dorenbos, and Val Zwiller. Superconducting nanowire single-photon detectors: A perspective on evolution, state-of-the-art, future developments, and applications. Applied Physics Letters , 118(19):190502, 05 2021.
- 7. D. Renker. Geiger-mode avalanche photodiodes, history, properties and problems. Nuclear Instruments and Methods in Physics Research Section A:

- Accelerators, Spectrometers, Detectors and Associated Equipment , 567(1):48-56, 2006. Proceedings of the 4th International Conference on New Developments in Photodetection.
- 8. B. Cabrera, R. M. Clarke, P. Colling, A. J. Miller, S. Nam, and R. W. Romani. Detection of single infrared, optical, and ultraviolet photons using superconducting transition edge sensors. Applied Physics Letters , 73(6):735-737, 08 1998.
- 9. Jascha A. Lau, Varun B. Verma, Dirk Schwarzer, and Alec M. Wodtke. Superconducting singlephoton detectors in the mid-infrared for physical chemistry and spectroscopy. Chem. Soc. Rev. , 52:921-941, 2023.
- 10. Heqing Wang, Hao Li, Lixing You, Yong Wang, Lu Zhang, Xiaoyan Yang, Weijun Zhang, Zhen Wang, and Xiaoming Xie. Fast and high efficiency superconducting nanowire single-photon detector at 630 nm wavelength. Appl. Opt. , 58(8):18681872, Mar 2019.
- 11. H. Hao, QY. Zhao, and YH. et al. Huang. A compact multi-pixel superconducting nanowire singlephoton detector array supporting gigabit space-toground communications. Light: Science &amp; Applications , 13(1):25, Jan 2024.
- 12. Han Chuen Lim, Akio Yoshizawa, Hidemi Tsuchida, and Kazuro Kikuchi. Wavelengthmultiplexed entanglement distribution. Optical Fiber Technology , 16(4):225-235, 2010.
- 13. WeiJun Zhang, LiXing You, Hao Li, Jia Huang, ChaoLin Lv, Lu Zhang, XiaoYu Liu, JunJie Wu, Zhen Wang, and XiaoMing Xie. Nbn superconducting nanowire single photon detector with efficiency over 90% at 1550 nm wavelength operational at compact cryocooler temperature. Science China Physics, Mechanics &amp; Astronomy , 60(12):120314, Oct 2017.
- 14. J. Chang, J. W. N. Los, J. O. Tenorio-Pearl, N. Noordzij, R. Gourgues, A. Guardiani, J. R. Zichi, S. F. Pereira, H. P. Urbach, V. Zwiller, S. N. Dorenbos, and I. Esmaeil Zadeh. Detecting telecom single photons with 99.5-2.07+0.5% system detection efficiency and high time resolution. APL Photonics , 6(3):036114, 03 2021.
- 15. F. Marsili, V. B. Verma, J. A. Stern, S. Harrington, A. E. Lita, T. Gerrits, I. Vayshenker, B. Baek, M. D. Shaw, R. P. Mirin, and S. W. Nam. Detecting single infrared photons with 93% system efficiency. Nature Photonics , 7(3):210-214, March 2013.
- 16. Lixing You. Superconducting nanowire singlephoton detectors for quantum information. Nanophotonics , 9(9):2673-2692, 2020.
- 17. L Redaelli, G Bulgarini, S Dobrovolskiy, S N Dorenbos, V Zwiller, E Monroy, and J M G´ erard. Design of broadband high-efficiency superconducting-nanowire single photon detectors. Superconductor Science and Technology , 29(6):065016, may 2016.
- 18. P. Neilinger, J. Greguˇ s, D. Manca, B. Granˇ ciˇ, c
- M. Kopˇ c´ ık, P. Szab´, P. Samuely, R. Hlubina, and o M. Grajcar. Observation of quantum corrections to conductivity up to optical frequencies. Phys. Rev. B , 100:241106, Dec 2019.
- 19. Samuel Kern, Pavol Neilinger, Magdal´ ena Pol´ˇ ckov´, a a Martin Bar´ anek, Tom´ aˇ s Plecenik, Tom´ aˇ Roch, and Miroslav Grajcar. s Optical and transport properties of nbn thin films revisited. arXiv preprint arXiv:2405.03704 , 2024.
- 20. Archan Banerjee, Robert M. Heath, Dmitry Morozov, Dilini Hemakumara, Umberto Nasti, Iain Thayne, and Robert H. Hadfield. Optical properties of refractory metal based thin films. Opt. Mater. Express , 8(8):2072-2088, Aug 2018.
- 21. Serhii Volkov, Maros Gregor, Tomas Roch, Leonid Satrapinskyy, Branislav Granˇ ciˇ, Tomas Fiantok, c and Andrej Plecenik. Superconducting properties of very high quality nbn thin films grown by pulsed laser deposition. Journal of Electrical Engineering , 70(7):89-94, 2019.
- 22. Alex L Efros and Michael Pollak. Electron-electron interactions in disordered systems . Elsevier, 2012.
- 23. COMSOL Multiphysics ® v. 6.2. www.comsol.com. COMSOL AB, Stockholm, Sweden.
- 24. Driessen, E. F. C., Braakman, F. R., Reiger, E. M., Dorenbos, S. N., Zwiller, V., and de Dood, M. J. A. Impedance model for the polarizationdependent optical absorption of superconducting single-photon detectors. Eur. Phys. J. Appl. Phys. , 47(1):10701, 2009.
- 25. Vikas Anant, Andrew J. Kerman, Eric A. Dauler, Joel K. W. Yang, Kristine M. Rosfjord, and Karl K. Berggren. Optical properties of superconducting nanowire single-photon detectors. Opt. Express , 16(14):10750-10761, Jul 2008.
- 26. Taro Yamashita, Shigehito Miki, Hirotaka Terai, and Zhen Wang. Low-filling-factor superconducting single photon detector with high system detection efficiency. Opt. Express , 21(22):27177-27184, November 2013.
Martin Baránek Pavol Neilinger Samuel Kern Miroslav Grajcar
Numerical modeling of SNSPD absorption utilizing optical conductivity with quantum corrections
Superconducting nanowire single-photon detectors are widely used in various
fields of physics and technology, due to their high efficiency and timing
precision. Although, in principle, their detection mechanism offers broadband
operation, their wavelength range has to be optimized by the optical cavity
parameters for a specific task. We present a study of the optical absorption of
a superconducting nanowire single photon detector (SNSPD) with an optical
cavity. The optical properties of the niobium nitride films, measured by
spectroscopic ellipsometry, were modelled using the Drude-Lorentz model with
quantum corrections. The numerical simulations of the optical response of the
detectors show that the wavelength range of the detector is not solely
determined by its geometry, but the optical conductivity of the disordered thin
metallic films contributes considerably. This contribution can be conveniently
expressed by the ratio of imaginary and real parts of the optical conductivity.
This knowledge can be utilized in detector design.