### Keywords

### Abstract

Terahertz time-domain spectroscopy is used to perform the first detailed studies of the electrodynamic properties of MoRe (60\( \% \)/40\( \% \)) films with thicknesses ranging from 10 to 100 nm. Films are prepared by magnetron sputtering technique on silicon substrates. The critical temperatures vary from \( 6.5\, \textrm{K} \) (for 10 nm film) to \( 9.5\, \textrm{K} \) (for 100 nm film). Spectra of complex permittivity, conductivity, refraction index, surface impedance and reflection coefficient for the films are acquired at frequencies \( 0.15 - 2.4\, \textrm{THz} \) (wavenumbers \( 5 – 80\, \textrm{cm}^{-1} \)) and in the temperature interval \( T=5\,–\, 300\,\textrm{K} \). For all films, temperature dependencies of the superconducting energy gap, penetration depth, superconducting condensate plasma frequency, and normalised superfluid density are obtained on a quantitative level. It is shown that the reduction of film thickness leads to a strong decrease of the critical temperature and magnitude of the energy gap. The observed suppression of superconductivity is assigned to reduction of the superconducting order parameter due to the contribution to the free energy of the electronic energy states at the surface of superconductor. The MoRe films with the obtained characteristics can be used in designing advanced superconducting electronic devices.

### Introduction

Molybdenum-rhenium superconducting alloys continue to attract the attention of researchers from a fundamental viewpoint because of their intriguing properties in the superconducting (SC) state, as well as their significant potential for use in modern electronics. Mo_{1−x}Re_{x} films exhibit superconductivity with a superconducting critical temperature reaching up to \( T_c=15\, \textrm{K} \) [1-4], which is of the order of magnitude higher than \( T_c =0.9\, \textrm{K} \) in Mo and \( T_c =1\, \textrm{K} \) in Re [5]. The compounds are s-wave type-\( II \) superconductors, with the upper critical field of \( 8\, \textrm{T} \) at \( 4\, \textrm{K} \) [6-8]. There are indications of a two-band/two-gap superconductivity in the MoRe alloys [9, 10] allowing to consider them as materials for the studies of multiband superconductivity. As an intriguing property, one can consider an electronic topological transition that is observed around the critical concentration \( x_c = 0.11 \) [11-13] making MoRe an ideal playground for study of time-reversal symmetry breaking in unconventional superconductors [14]. Thin MoRe films exhibit excellent superconducting properties [15], stable chemical, electrical and mechanical characteristics [15-17], high radiation stability [18], ability to make good electrical contact with carbon-based materials (graphene, carbon nanotubes) [19], stability for electron beam lithography [20] can be easily fabricated. All of these properties have revived a lively interest in MoRe alloys in the context of their use in developing of elements for modern electronic systems. As a transparent for carriers interface to graphene and carbon nanotubes [8, 19], MoRe films are promising for the development of elements for quantum computation [21, 22] and quantum memory [23, 24]. Few nanometers thick MoRe films can be used for the fabrication of SQUID-on-tip devices with enhanced characteristics [20, 25] and spin-qubit elements [26, 27].

The growing interest in thin MoRe films has emerged among communities working on developing sensitive single-photon detectors [15, 28]. These detectors commonly use NbN thin films with the critical temperature \( T_c \) close to \( 9\, \textrm{K} \) (for a 4nm thick film) [28]. With the SC doubled energy gap \( 2\Delta ≈ 6\, \textrm{meV} \) (for 15 nm thick film) [29], high quantum efficiency NbN-based detectors operating in visible and middle infrared spectra regions can be realized [28]. To achieve even lower, terahertz (THz) working frequencies, superconductors with smaller gap values are necessary. Numerous compounds fit this requirement, such as MoSi (\( 2\Delta(0) = 2.28\, \textrm{meV} \)), MoGe (\( 2\Delta(0) = 2.2\, \textrm{meV} \)), WSi (\( 2\Delta(0) = 1.52\, \textrm{meV} \)), and NbSi (\( 2\Delta(0) = 0.94\, \textrm{meV} \)), as reported in Ref. [30] (here (\( \Delta(0) \)) stands for the SC energy gap at zero temperature). However, their critical temperatures are rather low even for bulk samples. Thin MoRe films stands out because of their relatively high \( T_c \) values, which is close to that for NbN flims. Thus, detailed studies of the electrodynamic properties of MoRe films in sub-terahertz and terahertz frequency domains are needed. In this frequency range the temperature evolution of the superconducting gap can be spectroscopically traced in detail and analyzed. In addition, operating frequencies of electronic communication systems of the next generation indeed belong to this range. Since such studies have not been reported so far, our goal was to perform the first systematic measurements and analyses of the electrodynamic characteristics of Mo_{0.6}Re_{0.4} films with various thicknesses at frequencies \( f = 0.15 – 2.4\, \textrm{THz} \) (wavenumbers \( \nu=1/\lambda_{cm}=5–80\, \textrm{cm}^{-1} \), where \( \lambda_{cm} \) is radiation wavelength in cm) and in the temperature interval \( T=5 – 300\, \textrm{K} \).

### Results and Discussion

**Experimental procedure**

The Mo_{0.6}Re_{0.4} films were prepared by magnetron sputtering on highly insulating and thus transparent for terahertz radiation silicon substrates of thickness ≈0.5 mm, as described in Ref. [20]. Films thicknesses were 10, 15, 20, 30, 50, 70, and 100 nm. Terahertz (THz) characteristics of the films were determined using the commercial TeraView TPS Spectra 3000 time-domain spectrometer by measuring complex-valued (amplitude and phase) transmission coefficient spectra of the films on the substrates. Temperature-dependent THz spectra of the real parts of permittivity \( \varepsilon^{\prime}(\nu, T) \) and real conductivity \( \sigma^{\prime}(\nu, T) \) of the MoRe films were calculated with the TeraCalc software. This employs numerical solution of the system of two essentially nonlinear equations for the amplitude \( Tr(\varepsilon^{\prime}, \sigma^{\prime}) \) and the phase \( \varphi_{Tr} (\varepsilon^{\prime}, \sigma^{\prime}) \) of the complex transmission coefficient \( Tr(\varepsilon^{\prime}, \sigma^{\prime}) exp[i\varphi_{Tr} (\varepsilon^{\prime}, \sigma^{\prime})] \) of a two-layered system (film on a substrate), see Eqs. (1)-(5) in Supporting Information. Other quantities, such as the real and imaginary parts of the complex refractive index \( n^{∗} = n+ik \), and of the surface impedance \( Z_s = R_s +iX_s \), were also calculated using standard expressions (see, e. g., textbook [31]). THz parameters of the Si substrates were determined beforehand. Measurements were performed in a quasi-optical configuration of the measurement scheme of the spectrometer, with a homemade helium-flow cryostat equipped with Mylar windows.

**Normal state**

Figure 1(a) shows room temperature spectra of transmission coefficient \( Tr \) of the bare silicon substrate and of substrates with Mo_{0.6}Re_{0.4} films of 10, 20, 50 and 100 nm thick. The periodic oscillations in the spectra arise from the interference of the electromagnetic waves inside the plane-parallel transparent subresonators. The interval \( \Delta \nu \) between oscillations is determined mainly by the substrate refraction index \( n \) and thickness \( d, \Delta \nu = (2nd)^{−1} \), while the amplitudes of these peaks are determined by the transparency of the films. For the bare substrate, the transmission of the interferometric peaks reaches nearly unity, while the transmission of the peaks for the samples with MoRe films is significantly lower (at about two orders of magnitude for the 100 nm thick film). There is a phase shift at about \( \pi \) between the maxima observed for the bare substrate and the substrates covered with MoRe films. This is a consequence of impedance matching of the silicon-air interface by relatively thin conducting Mo_{0.6}Re_{0.4} film with a thickness smaller than the skin depth. The phenomenon was considered in detail in Refs. [32, 33].

Spectra of transmission coefficient of all samples as well as frequency dependences of conductivity and permittivity of all films did not reveal any significant changes during cooling from 300 K down to the critical temperature \( T_c \), in accordance with the literature data [35]. At all temperatures in the normal state, \( T > T_c \), the conductivity spectra of all films were dispersionless, indicating the metal-type response in the low frequency limit, \( \nu \ll \gamma \), where (\( \gamma \) is the carriers scattering rate) of the Drude conductivity model [36, 37]. This is demonstrated in Figure 1(b), where the conductivity spectra of two films, 10 nm and 50 nm thick, are shown as examples. Figure 1(c) shows the dependence of the normal state AC conductivity \( \sigma^{\prime} \), AC resistivity \( \rho^{\prime} = 1/\sigma^{\prime} \), and surface impedance \( (\sigma^{\prime} d)^{-1} \) on the thickness \( d \) of the MoRe films. It is seen that the dependences are well described by the following expression obtained for granular films accounting additional scattering effects at film surface and grain boundaries [38]

Here \( \rho_0 \) and \( l_0 \) are the bulk resistivity and the bulk mean free path of the charge carriers, respectively, and A is a constant of the order of unity. Fitting the \( \rho(d) \) dependence by Eq. ((1)) allows to determine \( l_0=10\, \textrm{nm} \) and \( \rho_0=75\, \mu \Omega \cdot \textrm{cm} \). The value \( l_0=10\, \textrm{nm} \) is more than two times larger than \( l=4.2\, \textrm{nm} \) obtained for MoRe alloy [14]. Taking the Fermi velocity \( \nu_F =2.34 \cdot 10^7\, \textrm{cm/s} \) [14], we estimate the scattering rate of charge carriers \( \gamma = \nu_F/(2\pi cl_0) \approx 120\, \textrm{cm}^{−1} \) (c is the speed of light). This value lies above our working frequency interval, thus confirming that the normal state response of the films corresponds to the low-frequency limit of the Drude conductivity model.

**Superconducting state**

Spectra of transmission coefficient, conductivity and permittivity of all studied films undergo strong changes when these films enter the SC state. This is demonstrated in Figure 2 by the spectra obtained for 10 nm-thick MoRe film on silicon substrate (see also the spectra for 50 nm-thick MoRe film in figure S1 in Supporting Information). These changes are typical for SC films as first it was demonstrated by Tinkham and coworkers [41-44] and reproduced later by other researchers in numerous experiments on conventional and high-\( T_c \) superconductors. As described below, by analyzing the spectra of the real parts of conductivity (\( \sigma^{\prime} \)) and permittivity (\( \varepsilon^{\prime} \)) through model fitting, we enable accurate determination of the critical temperature \( T_c \) (with an accuracy \( \pm 0.2\, \textrm{K}) \) and the energy gap \( \Delta \) (with an accuracy \( \pm 0.14\, \textrm{meV} \)), as well as other characteristics of the SC state of MoRe films. At the lowest achievable temperature of \( T=4.4\, \textrm{K} \), the envelope for interferometric peaks in the transmission spectrum (dashed red line in Figure 2(a)) displays broad maxima around the value \( 2\Delta /hc \) (where \( h \) is Planck’s constant), corresponding to the SC gap. The THz conductivity that characterizes absorption of radiation becomes suppressed at below \( 50 − 60\, \textrm{cm}^{−1} \), showing a kink around \( 2\Delta /hc \) (Figure 2(b)). Also, the dielectric permittivity strongly decreases below \( 20\, \textrm{cm}^{−1} \) (see Ref. [45]) according to the relationship \( \varepsilon^{\prime}=-\left(\nu_{p l}^{S C} / \nu\right)^2 \) representing dielectric response of SC condensate, where \( \nu_{p l}^{S C}=\frac{1}{2 \pi c} \sqrt{\frac{4 \pi n_e e^2}{m^*}} \) is the plasma frequency of SC condenstae, \( n_e \) is the density electrons, \( e \) is the electron charge, and \( m^∗ \) is the electron effective mass. The dielectric permittivity reaches large negative values that characterise the inductive response of Cooper-pairs condensate under the zero-frequency delta function in the conductivity spectrum. The superconducting transition also leads to strong changes in other electrodynamic quantities.

Figure 3 shows the temperature evolution of the spectra of real and imaginary parts of the complex refraction coefficient \( n^∗ = n + ik \) of the 10 nm-thick film (similar data for 50 nm-thick film is shown in figure S2). In the normal state \( n \approx k \sim (\sigma^{\prime}/\nu)^{0.5} \), as expected for a good conductor [36, 37]. Below \( T_c \), the extinction coefficient \( k \) strongly increases and refraction index \( n \) decreases, leading to a strong enhancement in the bulk reflection coefficient \( R \) below \( 2\Delta /hc \) (inset in Figure 3(b)). (*Note that an ideal s-wave superconductor at \( T = 0 \) and \( \nu < 2\Delta /hc \) would have zero AC conductivity \( \sigma^{\prime} \) and large negative permittivity \( \varepsilon^{\prime} \) values leading to \( n=0 \) and \( k \gg 1 \). This case would correspond to an ideal reflector with a bulk reflection coefficient [31] \( R_{bulk} = [(n − 1)^2 + k^2]/[(n+1)^2+k^2] \equiv (1+k^2)^2/(1+k^2)^2 = 1] \)*). Similar spectra for the 50 nm-thick film in figure S2. One can see in Figure 2(b) and Figure 3(a) that in the SC state at low wavenumbers, \( \nu < 2\Delta /hc \), there is significant scattering in the experimental data for the the conductivity \( \sigma^{\prime} \) and refraction index \( n \); it is especially pronounced for the 50 nm-thick film (figures S1, S2). The reason is that the transmission spectra measured in the SC state are mainly determined by the large negative permittivity \( \varepsilon^{\prime} \) and the large positive extinction coefficient \( k \). As a result, both quantities \( \varepsilon^{\prime} \) and \( k \) are determined much more precisely compared to \( \sigma^{\prime} \) and \( n \) [33, 35, 46].

To add more to the THz electrodynamic properties of the studied Mo_{0.6}Re_{0.4} films, we present in Figure 4 the spectra of the real \( R_s \) and imaginary \( X_s \) parts of surface impedance \( Z_s = R_s +iX_s \). For a thin film, it reads \( Z_s = (\sigma^∗d)^{−1} = \) \( \sigma^{\prime}[d(\sigma^{\prime 2}+\) \(\sigma^{\prime \prime 2})]^{−1}+i\sigma^{\prime \prime 2}[d(\sigma^{\prime 2}+\sigma^{\prime \prime 2})]^{−1} \), where \( \sigma^* = \sigma^{\prime} + i\sigma^{\prime \prime} \) is the complex-valued conductivity of the film material. At \( T > T_c \) we observe dispersionless dependences \( R_s \) and \( X_s \approx 0 \) on the wavenumber, as expected for a good conductor [37, 48]. In the SC state, the real impedance \( R_s \) (related to absorption) is strongly suppressed, while the frequency dependence of \( X_s \) has a broad maximum around the doubled gap frequency. These findings are typical fingerprints of a superconducting response [47].

Although MoRe compounds are considered as two-band superconductors [9, 10], we do not observe any corresponding feature in the measured THz spectra (Figure 2-Figure 4, as well as figure S1 and S2). The reason might be that the contribution to the electrodynamic response of one of the bands considerably exceeds the contribution of the other band. Thus, in order to model the experimental spectra, we used the BCS expressions of Zimmermann et al. [40] derived for a single-band superconductor. The fitting results shown by the solid lines in Figure 2-Figure 4, S1-S2 describe the experimental spectra quite well. From Figure 2(b) we see that the measured conductivity values are slightly above the BCS results, indicating extra below-gap absorption that may be due to bulk and/or surface defects. The performed fitting of the spectra allows us to determine temperature dependencies of characteristics of the SC condensate – namely, the plasma frequency \( \nu_{p l}^{S C} \), the London penetration depth \( \lambda_L = (2\pi \nu_{p l}^{S C})^{−1} \), the normalised superfluid density \( 1/\lambda_{L}^{2} \) and the energy gap Figure 5.

The temperature dependences of \( \nu_{p l}^{S C} \), \( \lambda_L \) and \( 1/\lambda_{L}^{2} \) are in agreement with the expressions of two-fluid superconductivity model:

The temperature dependencies of the gap energy values for the MoRe films of different thicknesses are well described by the BCS-like empirical expression [49]

In the above expressions ((2))-((4)) \( \nu_{p l}^{S C}(0) \), \( \lambda_L(0) \) and \( 2\Delta(0) \) represent zero-temperature values of the corresponding quantities. The doubled SC gap values obtained for all films are below \( 3\, \textrm{meV} \), so that the corresponding gap frequencies \( 2\Delta(0)/hc < 24\, \textrm{cm}^{−1} \) are smaller than the estimated scattering rate that was estimated above to be \( \gamma \approx 120\, \textrm{cm}^{−1} \). This means that we deal with the dirty limit superconductivity for all tested MoRe films.

Earlier, it was shown that the superfluid density in Mo_{0.6}Re_{0.4} polycrystalline sample cannot be explained within the framework of a single-band approach, since they show clear deviation from single-band behavior at low temperatures \( 2-6\, \textrm{K} \) [9]. Those data were described considering a two-band approach with the gaps \( 2\Delta_1=3.9\, \textrm{meV} \) and \( 2\Delta_2=1\, \textrm{meV} \). The gap value obtained in our experiment for the thickest (100 nm) film, \( 2\Delta(0)=2.9\, \textrm{meV} \), falls between \( 2\Delta_1 \) and \( 2\Delta_2 \).

**Thickness dependence of THz characteristics of Mo _{0.6}Re_{0.4} films in superconducting state. **

Figure 6 presents the dependencies of the penetration depth \( \lambda_L(0) \), the critical temperature Tc and the SC gap energy \( 2\Delta(0) \) on the thickness \( d \) of the MoRe. Reduction of the film thickness leads to an increase of the penetration depth and the decrease in both the critical temperature and zero-temperature energy gap. The ratio \( 2\Delta(0)/k_BT_c = 3.8 \pm 0.2 \) is slightly higher than the standard BCS value 3.52 for weak-coupling regime, and this indicates moderately strong electron-phonon coupling in the MoRe films (\( k_B \) is the Boltzmann constant). By fitting the experimental dependence of \( \lambda_L(0) \) on \( d \) with using the expression [50]

which is valid in the dirty limit, we are able to estimate the bulk BCS coherence length \( \xi_0 = 26\, \textrm{nm} \); \( \lambda_0 \) corresponds to the bulk penetration depth. The result of the fitting is shown in Figure 6(a) by solid line. Taking \( 2\Delta(0) = 3\, \textrm{meV} \) as typical energy gap of a bulk material (according to Figure 6(c)) and the Fermi velocity [14] \( \nu_F = 2.34 \cdot 10^7\, \textrm{cm/s} \) and using the well-known expression [50] \( \xi_0 = ħ\nu_F /(\pi \Delta) \) we get the value of the BCS coherence length \( \xi_0 \approx 33\, \textrm{nm} \) that is close to the above estimate \( \xi_0=26\, \textrm{nm} \). It makes possible to estimate the Ginsburg-Landau parameter, we obtain \( κ = \lambda_L/\xi \sim 10 \), indicating type-\( II \) superconductivity in MoRe films. Using the dirty-limit expression [50]

we can evaluate the dependence of the coherence length \( \xi \) on film thickness. The result is shown in the inset of Figure 6(a) together with the thickness dependence of the mean free path \( l \). For all films, we have \( \xi(d) > l(d) \), what again confirm the dirty limit SC. Smallness of the coherence length \( \xi \) and the mean free path \( l \) with respect to both the penetration depth in normal state \( \delta = (2\pi k \nu)^{-1} \) and in SC state \( \lambda_L \) (see insets in Figure 2(c), S1c, Figure 5(b)) means that the local approach holds in the normal and SC states of the studied Mo_{0.6}Re_{0.4} films.

Let us consider the suppression of superconductivity in thin films that is accompanied by a reduction of both the critical temperature and the energy gap. Since Mo_{0.6}Re_{0.4} films are binary alloy, such suppression could be due to disorder. Maekawa and Fukuyama [51] were the first to analyze such suppression caused by Cooper pairs localization in disordered thin films that leading to an enhanced Coulomb repulsion and suppressed electronic density of states. They obtained the following expression

where \( g^{\prime} < 0 \) stands for the attractive interaction responsible for the Cooper pair formation, \( g_1 > 0 \) characterizes Coulomb repulsion, \( N(0) \) is electronic density of states at the Fermi level per spin, \( R_{\square} \) is sheet resistance, \( T^{bulk}_c \) is the critical temperature of a bulk sample. The expression (7) with the parameters \( g^{\prime} = -0.2 \) and \( g_1 = 7 \) describes the experimental \( T_c(d) \) data quite well (dashed line in Figure 6(b)), but the value \( g_1 = 7 \) seems to be too large as compared with \( g_1 \approx 1 \) considered in Maekawa and Fukuyama’s paper.

Alternatively, we can use the theory of Finkel’shtein [52] based on the renormalization group approach and leading to the following expression

where \( t=\frac{e^2}{2 \pi^2 \hbar} R_{\square} \), \( \gamma = 1/log(T^{bulk}_c\tau) \), and \tau is the transit time. Taking the Fermi velocity \( \nu_F=2.34 \cdot 10^{7}\, \textrm{cm/s} \) [14] and thickness dependent mean free path (inset of Figure 6(a)), we conclude that the expression ((8)) cannot reproduce our experimental data (see red line in the inset of Figure 6(b)), since we should consider unreasonably short bulk mean free path, \( l_0 = 0.56\, \textrm{Å} \) (see blue line in the inset of Figure 6(b)).

Following the paper Simonin [53] we consider a surface contribution (*The author adds the term proportional to an integral \( \int_S\psi^2ds \) over the surface of a superconductor into the free-energy functional, where \( \psi \) is the normalized order parameter. This leads to the modification of the boundary conditions and results in the suppression of \( T_c \)*) to the Ginzburg-Landau free energy functional that leads to decrease in the density of states near the surface. This effect becomes significant in thin films where the surface-to-volume ratio is dominant. The theory of Simonin predicts the following dependence

where \( d_c = 2a/N(0)V \) is the critical thickness, corresponding to \( T_c = 0 \), \( a \) is the Thomas-Fermi screening length that is of the order of lattice parameter, \( V \) is the interaction potential and \( N(0)V \) is the bulk interaction potential. The expression (9) describes well the experimental data with the ratio \( 2a/N(0)V = 1.2\, \textrm{nm} \) (solid line in Figure 6(b)). Taking for Mo_{0.6}Re_{0.4} \( a \approx 0.3\, \textrm{nm} \), we obtain a quite reasonable value for the interaction potential \( N(0)V \approx 0.5 \). This allows us to conclude that the reduction of the order parameter near the surfaces of SC films can be considered as a more appropriate mechanism of suppression of superconductivity in the studied films. We note that similar results have been obtained for pure molybdenum films [54].

According to Figure 6(c), the ratio \( 2\Delta/k_BT_c = 3.8 \pm 0.2 \) seems to be thickness-independent, suggesting that the critical temperature depends on the film thickness in the same way as the energy gap. This is demonstrated by a fit of the \( 2\Delta(d) \) dependence with an expression \( 2\Delta(d) \sim (1 − const/d) \) that is similar to Eq. (9). We thus suggest that the suppression of both, the critical temperature and the SC energy gap are governed by the same mechanism – the alteration of the boundary conditions caused by the surface term introduced in the GL functional, resulting in change of the interaction potential \( N(0)V \).

### Conclusion

The first systematic studies of terahertz electrodynamic properties of Mo_{0.6}Re_{0.4} films of thicknesses ranging from 10 nm to 100 nm and corresponding critical temperatures from \( 6.5\, \textrm{K} \) to \( 9.5\, \textrm{K} \) are performed at frequencies \( 0.15 – 2.4\, \textrm{THz} \) and in the temperature interval \( T=5 – 300\, \textrm{K} \). The spectra of conductivity and permittivity are described within the BCS single-band approach. The obtained ratio \( 2\Delta/k_BT_c = 3.8 \pm 0.2 \) slightly exceeds the BCS value of 3.52, demonstrating moderately strong electron-phonon coupling. The temperature dependencies of the superconducting energy gap, the penetration depth, the superconducting condensate plasma frequency and the normalized superfluid density are reproduced within the two-fluid superconductivity model. A strong decrease of the critical temperature and the energy gap in thin MoRe films is associated with the suppression of the superconducting order parameter due to a decrease in the local electronic density of states in MoRe films near their surfaces. Although our thinnest films have slightly smaller values of the critical temperature than those obtained earlier [28], the values \( T_c(d) \) are still higher than the critical temperatures of the alloys [30] MoSi, MoGe, WSi, NbSi, mentioned in the Introduction section. This makes MoRe alloy a perspective material for the design of elements and devices of advanced electronics and optoelectronics, for example, single-photon detectors, high-Q resonators for frequency locking systems (see figure S3), etc.

### Acknowledgement

This work was supported by RSF 23-72-30004 (terahertz experiments and data analysis). Magnetron sputtering fabrication was supported by the Ministry of Science and Higher Education of the Russian Federation (No. FSMG-2022-0014).

### Supporting information and data availability

The Supporting Information is available free of charge at https://jmsn.colab.ws/publications/01-01002/supporting-file.

The data that support the findings of this study are available from the corresponding authors upon request.

### Contact information

Corresponding author: Boris M. Nekrasov, https://orcid.org/0000-0003-3324-2893, e-mail nekrasov.bm@phystech.edu.

Corresponding author: Vasily S. Stolyarov, https://orcid.org/0000-0002-5317-0818, e-mail stolyarov.vs@phystech.edu.

### Author contribution

Elena S. Zhukova, Alexander V. Melentev and Lenar S. Kadyrov carried out the THz experiments, spectra processing and data analysis. Boris M. Nekrasov, Anton S. Shaimardanov, Andrey G. Shishkin, Aleksander A. Golubov, Mikhail Yu. Kupriyanov performed data analysis, Boris P. Gorshunov and Vasily S. Stolyarov conceived and supervised the work. All authors contributed to the manuscript.

### Competing Interests

The Authors declare no Competing Financial or Non-Financial Interests.

### References

Optical Properties of Metals

Introduction to superconductivity