Measuring Complex Refractive Indices of a Nanometer-Thick Superconducting Film Using Terahertz Time-Domain Spectroscopy with a 10 Femtoseconds Pulse Laser

Superconducting thin films are widely applied in various fields, including switching devices, because of their phase transition behaviors in relation to temperature changes. Therefore, it is important to quantitatively determine the optical constant of a superconducting material in the thin-film state. We performed a terahertz time-domain spectroscopy, based on a 10 femtoseconds pulse laser, to measure the optical constant of a superconducting GdBa2Cu3O7−x (GdBCO) thin film in the terahertz region. We then estimated the terahertz refractive indices of the 70 nm-thick GdBCO film using a numerical extraction process, even though the film thickness was approximately 1/10,000 times smaller than the terahertz wavelength range of 200 μm to 1 mm. The resulting refractive indices of the GdBCO thin film were consistent with the theoretical results using the two-fluid model. Our work will help to further understand the terahertz optical properties of superconducting thin films with thicknesses under 100 nm, as well as provide a standard platform for characterizing the optical properties of thin films without the need of Kramers–Kronig transformation at the terahertz frequencies.

Recently, many studies have been carried out in order to develop new switching devices combined with phase transition materials modulated via external stimuli, such as temperature and electrical voltage [9][10][11][12]. In addition, superconducting materials are widely applied in the field of switching devices because of their phase transition below the critical temperature T c [13][14][15][16][17]. In the terahertz frequency region, such materials are much more promising for developing metamaterials and switching devices because the superconductor properties remain robust and well-maintained [18]. Furthermore, it is important to ensure that the superconducting property is well maintained in the desired thin film form, because superconductors are used in the form of thin films in most terahertz applications [19,20].
One advantage of employing time-domain spectroscopy is that the optical constants of a target material can be determined via a transmission or reflection experiment without the need of Kramers-Kronig transformation, because the amplitude and phase information can be obtained simultaneously [21]. Of these, analysis of the results is relatively simpler in the transmission experiment than in the reflection experiment. In the case of a thick substrate there is a sufficient phase difference; hence, the refractive index can be easily obtained.
However, phase retardation of a terahertz-transparent thin film is not sufficient, and the refractive index is difficult to obtain owing to the pulse width resolution of the probe beam [22]. In contrast, because the film in a metal has a sufficiently large refractive index, the phase difference can be distinguished even in the case of a thin film. Therefore, the complex refractive index of the metal sample can be determined [23][24][25]. More specifically, if nd/c > PW, where n is the real part of the refractive index of the material, d is the thickness, c is the speed of light, and PW is the pulse width of the probe beam, we can estimate the complex refractive indices of the material. For example, a probe beam with a pulse width of 10 fs can be used to analyze a material with the optical path (=nd) of 3 µm. In many terahertz time-domain spectroscopy studies, lasers with a pulse width of about 100 fs were used [19,26,27]. The corresponding resolution of the optical path by the 100 fs pulse is about 30 µm, which is too long to directly measure any change in the phase of the terahertz pulse through the thin film with a short optical path (nd/c PW). One study estimated the complex conductivities of metal thin films with a thicknesses of 10 to 30 nm, using a laser with a pulse width of 12 fs [28]. In the case of a conducting material, it would be relatively easier to estimate its refractive index, even within the form of thin films, but we need one more procedure to accurately estimate complex refractive indices of superconductor thin films, which retain relatively low refractive indices below T c .
In this work, we introduced a method for measuring the complex refractive indices of a target material using terahertz time-domain spectroscopy combined with a numerical extraction process. Consistent with previous studies, we double-checked the complex refractive indices of an optically thick insulating substrate (undoped silicon and LaAlO 3 ) and metallic thin films [25,27,29,30]. Based on our ability to analyze the refractive indices of both insulating and conducting materials, we estimated the terahertz refractive indices of a 70 nm-thick superconducting GdBa 2 Cu 3 O 7−x (GdBCO) thin film with a critical temperature of around 90 K in a broad temperature range of 20-180 K.

Terahertz Time-Domain Spectroscopy
To experimentally determine the complex refractive indices of the samples, terahertz time-domain spectroscopy was performed (Figure 1). A terahertz pulse with a singlecycle picosecond pulse was generated from a commercial GaAs emitter (Tera-SED3, Laser Quantum, Stockport, UK), illuminated with a Ti:sapphire laser pulse train with the center wavelength of 780 nm, repetition rate of 75.1 MHz, and pulse width of 10 fs (Synergy, Spectra-Physics, Andover, MA, USA). For the detection, we used the electro-optic sampling method with a 1 mm-thick ZnTe crystal [31], and the resulting signal-to-noise ratio (SNR) is 5000:1. As mentioned in the introduction, the 10 fs probe pulse in the electro-optic sampling method has limited to an optical path resolution of 3 um.

Sample Fabrication
Silver thin films with thicknesses of 15, 30 and 50 nm, which are sub-skin depth in the terahertz region, were deposited via e-beam evaporation onto a 500 µm-thick undoped silicon substrate. The deposition conditions for the silver thin film included a degree of vacuum less than 5.0 × 10 −6 Torr and deposition rate of 1 Å/s. Additionally, a 70 nm-thick GdBCO thin film was deposited onto a 500 µm-thick single-crystal LaAlO 3 (LAO) substrate using pulsed laser deposition (KrF excimer laser, LPX pro240, 248 nm, Coherent, Santa Clara, USA). The deposition conditions for the GdBCO thin film were optimized for the best epitaxial quality as follows: GdBCO was synthesized at 780 • C under the partial O 2 pressure of 350 mTorr with the laser energy of 130 mJ/pulse; subsequently, oxygen annealing was performed at 500 mTorr and 500 • C for 1 h. As shown in Figure 2, the GdBCO film thickness is 70 ± 6 nm confirmed by scanning electron microscope (SEM).

Sample Fabrication
Silver thin films with thicknesses of 15, 30 and 50 nm, which are sub-skin depth in the terahertz region, were deposited via e-beam evaporation onto a 500 μm-thick undoped silicon substrate. The deposition conditions for the silver thin film included a degree of vacuum less than 5.0 × 10 −6 Torr and deposition rate of 1 Å/s. Additionally, a 70 nm-thick GdBCO thin film was deposited onto a 500 μm-thick single-crystal LaAlO3 (LAO) substrate using pulsed laser deposition (KrF excimer laser, LPX pro240, 248 nm, Coherent, Santa Clara, USA). The deposition conditions for the GdBCO thin film were optimized for the best epitaxial quality as follows: GdBCO was synthesized at 780 °C under the partial O2 pressure of 350 mTorr with the laser energy of 130 mJ/pulse; subsequently, oxygen annealing was performed at 500 mTorr and 500 °C for 1 h. As shown in Figure 2, the GdBCO film thickness is 70 ± 6 nm confirmed by scanning electron microscope (SEM).

Sample Fabrication
Silver thin films with thicknesses of 15, 30 and 50 nm, which are sub-skin depth in the terahertz region, were deposited via e-beam evaporation onto a 500 μm-thick undoped silicon substrate. The deposition conditions for the silver thin film included a degree of vacuum less than 5.0 × 10 −6 Torr and deposition rate of 1 Å/s. Additionally, a 70 nm-thick GdBCO thin film was deposited onto a 500 μm-thick single-crystal LaAlO3 (LAO) substrate using pulsed laser deposition (KrF excimer laser, LPX pro240, 248 nm, Coherent, Santa Clara, USA). The deposition conditions for the GdBCO thin film were optimized for the best epitaxial quality as follows: GdBCO was synthesized at 780 °C under the partial O2 pressure of 350 mTorr with the laser energy of 130 mJ/pulse; subsequently, oxygen annealing was performed at 500 mTorr and 500 °C for 1 h. As shown in Figure 2, the GdBCO film thickness is 70 ± 6 nm confirmed by scanning electron microscope (SEM).

Transmission Coefficients of Bulk and Thin Film Samples
Using the Fresnel equation, we can obtain the transmission and reflection of light of the media using the refractive index, and vice versa. As shown in Figure 3, the Fresnel equation is defined as follows [32].

Transmission Coefficients of Bulk and Thin Film Samples
Using the Fresnel equation, we can obtain the transmission and reflection of light of the media using the refractive index, and vice versa. As shown in Figure 3, the Fresnel equation is defined as follows [32]. When the electric field is normally incident to material j from i, the Fresnel coefficients are: The electric field of light passing through the interface is given by: When the electric field passes through the inside of the medium, the electric field is exponentially decayed as follows where is a wavevector in the medium, expressed as: According to Equations (4) and (5), it can be seen that the real part of the refractive index mainly contributes to the phase changes, and the imaginary part of the refractive index mainly contributes to the amplitude changes. Using these relationships, the transmitted electric fields for the bulk and thin film samples can be calculated as follows.
In time-domain spectroscopy experiments, the first and multiple signals for bulk samples, such as single crystals, can be distinguished from the time data. Therefore, there is no need to consider multiple reflections. When the electric field is normally incident to material j from i, the Fresnel coefficients are: The electric field of light passing through the interface is given by: When the electric field passes through the inside of the medium, the electric field is exponentially decayed as follows where k k is a wavevector in the medium, expressed as: According to Equations (4) and (5), it can be seen that the real part of the refractive index mainly contributes to the phase changes, and the imaginary part of the refractive index mainly contributes to the amplitude changes. Using these relationships, the transmitted electric fields for the bulk and thin film samples can be calculated as follows.
In time-domain spectroscopy experiments, the first and multiple signals for bulk samples, such as single crystals, can be distinguished from the time data. Therefore, there is no need to consider multiple reflections.
In the experiment, the case was first measured with only the aperture, and then the sample was measured. Subsequently, the data measured on the sample was divided by the data measured with only the aperture. For the transmitted electric field of the aperture case, as shown in Figure 4a, the sample thickness was considered.

of 18
In the experiment, the case was first measured with only the aperture, and then the sample was measured. Subsequently, the data measured on the sample was divided by the data measured with only the aperture. For the transmitted electric field of the aperture case, as shown in Figure 4a, the sample thickness was considered.
Thus, the transmission coefficient of the bulk sample is: For thin films, the terms reflected from the interface (so-called multiple reflections or Fabry-Perot interference) must be considered as the multiple reflections, and the first signal cannot be distinguished from the time data. Thin films are generally deposited onto substrates, and the interfaces of the materials are shown in Figure 5.
Thus, the transmission coefficient of the bulk sample is: For thin films, the terms reflected from the interface (so-called multiple reflections or Fabry-Perot interference) must be considered as the multiple reflections, and the first signal cannot be distinguished from the time data. Thin films are generally deposited onto substrates, and the interfaces of the materials are shown in Figure 5. the data measured with only the aperture. For the transmitted electric field of the apertur case, as shown in Figure 4a, the sample thickness was considered.
Thus, the transmission coefficient of the bulk sample is: For thin films, the terms reflected from the interface (so-called multiple reflections o Fabry-Perot interference) must be considered as the multiple reflections, and the first sig nal cannot be distinguished from the time data. Thin films are generally deposited ont substrates, and the interfaces of the materials are shown in Figure 5.  Fabry-Perot interference is defined as the multiple reflections that occur at both interfaces of a thin sample, as shown in Equation (9).
Crystals 2021, 11, 651 6 of 18 Thus, the transmitted electric field through the thin sample can be established as: In the experiment, the situation will be as shown in Figure 6.
Fabry-Perot interference is defined as the multiple reflections that occur at both interfaces of a thin sample, as shown in Equation (9).
Thus, the transmitted electric field through the thin sample can be established as: In the experiment, the situation will be as shown in Figure 6. Typical data processing for calculating the transmission coefficient in an experiment involves dividing the transmitted signal through the thin film on the substrate ( Figure 6) by that of the substrate only ( Figure 4b). Therefore, the transmission coefficient of the thin film is given by:

Numerical Method
Most of the materials used as a substrate have very little absorption in the terahertz region. In particular, when the refractive index of the substrate is = + , the K3 value is very small compared with the N3 value. Then, the following approximation can be obtained [33]: The transmitted electric field can be expressed as Equation (13), including the amplitude and phase components.
Combining Equations (12) and (13) provides the following: Typical data processing for calculating the transmission coefficient in an experiment involves dividing the transmitted signal through the thin film on the substrate ( Figure 6) by that of the substrate only ( Figure 4b). Therefore, the transmission coefficient of the thin film is given by:

Numerical Method
Most of the materials used as a substrate have very little absorption in the terahertz region. In particular, when the refractive index of the substrate is n 3 = N 3 + iK 3 , the K 3 value is very small compared with the N 3 value. Then, the following approximation can be obtained [33]: The transmitted electric field can be expressed as Equation (13), including the amplitude and phase components.
Combining Equations (12) and (13) provides the following: Therefore, the equations for the refractive indices of thick and insulating materials can be obtained, respectively, as: These values are almost correct; however, to improve the accuracy, these values are used as the initial value, and the numerical method is performed to determine the result faster and more accurately, as shown in Figure 7. We used the fsolve function in MATLAB for our analysis. The fsolve function is an optimization tool that finds the optimal solution using the Levenberg-Marquardt algorithm [34,35]. Using the optimization tool, we finally obtained the most accurate refractive indices, which were perfectly matched with the amplitude and phase of the transmitted signal through the sample.
Therefore, the equations for the refractive indices of thick and insulating materials can be obtained, respectively, as: These values are almost correct; however, to improve the accuracy, these values are used as the initial value, and the numerical method is performed to determine the result faster and more accurately, as shown in Figure 7. We used the fsolve function in MATLAB for our analysis. The fsolve function is an optimization tool that finds the optimal solution using the Levenberg-Marquardt algorithm [34,35]. Using the optimization tool, we finally obtained the most accurate refractive indices, which were perfectly matched with the amplitude and phase of the transmitted signal through the sample. In the case of conducting samples, the n and values of the complex refractive index are generally comparable in the terahertz frequencies. For this reason, to set an appropriate initial value when calculating the numerical solution, the refractive indices are considered to be n = when the refractive index of the film is = + . If n = , we can assume that there is only one variable. Therefore, Equation (18) was established, which can be solved using only the amplitude data.
Following this, a reasonable solution of the complex refractive index can be obtained using the same optimization tool with this initial value. In the case of conducting samples, the n and κ values of the complex refractive index are generally comparable in the terahertz frequencies. For this reason, to set an appropriate initial value when calculating the numerical solution, the refractive indices are considered to be n = κ when the refractive index of the film is n f ilm = n + iκ. If n = κ, we can assume that there is only one variable. Therefore, Equation (18) was established, which can be solved using only the amplitude data.
Following this, a reasonable solution of the complex refractive index can be obtained using the same optimization tool with this initial value. Silicon wafers are one of the most commonly used substrates. In Section 4.2, an undoped silicon wafer (2000 Ω) was used as a substrate for the silver thin films. The time data of the terahertz transmitted through the aperture and the Si wafer are shown in Figure 8, and the normalized Fourier transform data are shown in Figure 9. By substituting the amplitude and phase of the Fourier transform data into Equations (16) and (17), the refractive index values can be obtained. Using these as an initial value, we obtained the numerical solution of Equation (6) to determine the refractive indices, as shown in Figure 10 (n ≈ 3.4, κ ≈ 0.005). This result is consistent with the generally known refractive indices of high resistivity: float-zone silicon [29].

Undoped Silicon and Single Crystal LaAlO3
Silicon wafers are one of the most commonly used substrates. In Section 4.2, an undoped silicon wafer (2000 Ω) was used as a substrate for the silver thin films. The time data of the terahertz transmitted through the aperture and the Si wafer are shown in Figure 8, and the normalized Fourier transform data are shown in Figure 9. By substituting the amplitude and phase of the Fourier transform data into Equations (16) and (17), the refractive index values can be obtained. Using these as an initial value, we obtained the numerical solution of Equation (6) to determine the refractive indices, as shown in Figure  10 ( ≈ 3.4, ≈ 0.005). This result is consistent with the generally known refractive indices of high resistivity: float-zone silicon [29].

Undoped Silicon and Single Crystal LaAlO3
Silicon wafers are one of the most commonly used substrates. In Section 4.2, an undoped silicon wafer (2000 Ω) was used as a substrate for the silver thin films. The time data of the terahertz transmitted through the aperture and the Si wafer are shown in Figure 8, and the normalized Fourier transform data are shown in Figure 9. By substituting the amplitude and phase of the Fourier transform data into Equations (16) and (17), the refractive index values can be obtained. Using these as an initial value, we obtained the numerical solution of Equation (6) to determine the refractive indices, as shown in Figure  10 ( ≈ 3.4, ≈ 0.005). This result is consistent with the generally known refractive indices of high resistivity: float-zone silicon [29].    A 500 μm-thick single-crystal LAO was used to deposit the GdBCO film. Usually an SrTiO3 (STO) substrate is used, but in the terahertz region transmission is very low, owing to the influence of its phonon, thus an STO is not suitable [36,37]. The time data of the transmission through the aperture and LAO are shown in Figure 11, and the normalized Fourier transform data are shown in Figure 12. Meanwhile, Figure 13 presents the refractive indices of LAO, obtained using the same method as that used to calculate the refractive indices of the Si wafer. Therefore, we obtained the refractive indices of LAO ( = 4.8, ≈ 0.01). This result is consistent with the generally known refractive indices of single crystal LAO [30].  A 500 µm-thick single-crystal LAO was used to deposit the GdBCO film. Usually an SrTiO 3 (STO) substrate is used, but in the terahertz region transmission is very low, owing to the influence of its phonon, thus an STO is not suitable [36,37]. The time data of the transmission through the aperture and LAO are shown in Figure 11, and the normalized Fourier transform data are shown in Figure 12. Meanwhile, Figure 13 presents the refractive indices of LAO, obtained using the same method as that used to calculate the refractive indices of the Si wafer. Therefore, we obtained the refractive indices of LAO (n = 4.8, κ ≈ 0.01). This result is consistent with the generally known refractive indices of single crystal LAO [30]. A 500 μm-thick single-crystal LAO was used to deposit the GdBCO film. Usually an SrTiO3 (STO) substrate is used, but in the terahertz region transmission is very low, owing to the influence of its phonon, thus an STO is not suitable [36,37]. The time data of the transmission through the aperture and LAO are shown in Figure 11, and the normalized Fourier transform data are shown in Figure 12. Meanwhile, Figure 13 presents the refractive indices of LAO, obtained using the same method as that used to calculate the refractive indices of the Si wafer. Therefore, we obtained the refractive indices of LAO ( = 4.8, ≈ 0.01). This result is consistent with the generally known refractive indices of single crystal LAO [30].

Silver Thin Films
Silver is a noble metal with very high conductivity. Thus, it is widely applied in the field of metamaterials and plasmonics as a perfect metal in the terahertz and lower frequency region. Noble metals have very large refractive indices in the terahertz region, so even if it is thin, the optical path is large enough. As a result, there is no problem in estimating the refractive index, even with the form of the thin film. Figure 14 shows the measured time trace data of the silver films and the substrate. Using the data in Figure 15, obtained from the Fourier transform of the data in Figure 14, the refractive indices of the silver thin film can be estimated, as shown in Figure 16, by applying the numerical extraction method presented in Section 3.

Silver Thin Films
Silver is a noble metal with very high conductivity. Thus, it is widely applied in the field of metamaterials and plasmonics as a perfect metal in the terahertz and lower frequency region. Noble metals have very large refractive indices in the terahertz region, so even if it is thin, the optical path is large enough. As a result, there is no problem in estimating the refractive index, even with the form of the thin film. Figure 14 shows the measured time trace data of the silver films and the substrate. Using the data in Figure 15, obtained from the Fourier transform of the data in Figure 14, the refractive indices of the silver thin film can be estimated, as shown in Figure 16, by applying the numerical extraction method presented in Section 3.

Silver Thin Films
Silver is a noble metal with very high conductivity. Thus, it is widely applied in the field of metamaterials and plasmonics as a perfect metal in the terahertz and lower frequency region. Noble metals have very large refractive indices in the terahertz region, so even if it is thin, the optical path is large enough. As a result, there is no problem in estimating the refractive index, even with the form of the thin film. Figure 14 shows the measured time trace data of the silver films and the substrate. Using the data in Figure 15, obtained from the Fourier transform of the data in Figure 14, the refractive indices of the silver thin film can be estimated, as shown in Figure 16, by applying the numerical extraction method presented in Section 3.     The simple Drude model for conductors is useful for describing the properties of metals in the low-frequency region. The Drude model dielectric function is given by The simple Drude model for conductors is useful for describing the properties of metals in the low-frequency region. The Drude model dielectric function is given by where ω p = 1.37 × 10 16 Hz and ω τ = 27.4 × 10 12 Hz are the Drude parameters of a bulk silver [38]. As shown in Figure 16, it is possible to estimate the refractive indices of the sub-50 nmthick silver films, which are in a sub-skin depth region. It can be observed that the feature for this thin film tends to differ from that of bulk metal. Compared with the n and κ values extracted from the Drude model of the bulk characteristic, it is apparent that both these values decrease as the film becomes even thinner. In particular, κ decreases relatively more. In addition, the refractive indices of the 15 nm-thick silver film are approximately half that of the bulk. It should be noted that the terahertz transmittance through the metal film of the sub-skin depth thickness increases even more as its thickness decreases, because the refractive indices of the metal film also decline with the decreasing thickness. As the optical path of the silver thin film with a thickness of 15 nm is approximately 400 × 15 nm = 6 µm, the refractive indices of the silver film as thin as 10 nm can be experimentally extracted accurately. The estimated indices of n and κ in the sub-skin depth thick silver films are qualitatively consistent with previous results [25,27].

GdBCO Thin Film
GdBCO is a high-T c cuprate superconducting material with a T c of about 90 K. Such materials are more intriguing in various fields because they exhibit phase transition characteristics at moderately high temperatures. We measured transmitted signals through the 70 nm-thick GdBCO film on the LAO substrate in the temperature range between 20 and 180 K (Figure 17). By Fourier-transforming the time-domain data, we obtained the transmitted amplitude and phase spectra for the sample and a bare substrate. For normalization of the measured spectra, we used a reference signal for the bare LAO substrate. We obtained the normalized transmitted amplitude and phase spectra as the ratio between the measured data for the sample and the bare substrate, as shown in Figure 18. Using the same method as was used for the analysis of the silver thin films described in Section 4.2., we estimated the refractive indices of the GdBCO thin film at the selected temperatures, as shown in Figure 19.     As shown in Figure 19, the GdBCO thin films exhibit a Drude feature above Tc, like a general conducting material, and superconducting features at temperatures below Tc. We will further explain the details of the two-fluid model in the following paragraph.
To analyze the optical properties of the superconducting thin film, a two-fluid model was used (Figure 20a,b). The two-fluid model is similar to the Drude model, except that it explains the conductivity characteristics of a superconductor such that the ratio of carriers (normal electrons and superconducting electrons) changes with temperature [39]. The two-fluid model conductivity is expressed as follows As shown in Figure 19, the GdBCO thin films exhibit a Drude feature above T c , like a general conducting material, and superconducting features at temperatures below T c . We will further explain the details of the two-fluid model in the following paragraph.
To analyze the optical properties of the superconducting thin film, a two-fluid model was used (Figure 20a,b). The two-fluid model is similar to the Drude model, except that it explains the conductivity characteristics of a superconductor such that the ratio of carriers (normal electrons and superconducting electrons) changes with temperature [39]. The two-fluid model conductivity is expressed as follows where f n and f s are ratios of the normal and superconducting carriers, respectively, given by where τ and t are the relaxation time and the reduced temperature T/T c , respectively, and γ is an empirical parameter.  Here, α = 10, γ = 1.5, T c = 90 K, ω p = 2 × 10 15 Hz, and τ(1) = 3.57 × 10 −14 s. The superconductivity emerges on the CuO 2 plane in cuprate superconductors. The A site ion (Y, Gd), located in between two CuO 2 layers, plays the role of the insulating layer. The two-fluid model is based on normal and superconducting carriers, and since the two-fluid model has been mainly applied for YBa 2 Cu 3 O 7 (YBCO), we used the parameters of YBCO, which is similar to the GdBCO [40], in the analysis using the two-fluid model [41].
The complex dielectric function is expressed as: The complex refractive index is given by: In the time data of Figure 17, the tendency of the terahertz transient appears at an earlier delay time, as the temperature decrease is contrary to the belief that the conductivity increases and the refractive index increases as the temperature decreases. This is because, as the temperature decreases, the real part of the refractive index decreases, which shortens the optical path. In addition, the imaginary part of the refractive index will greatly affect the amplitude of the exponential term in the Fresnel coefficient. However, as it also exists in the term outside the exponential, it may affect the optical path, although it is less than the real part of the refractive index. As shown in Figure 20, the 70 nm-thick superconducting GdBCO film exhibited qualitatively similar properties. However, the reason why the quantitative values are different not only owes to the use of the YBCO parameters, but also the GdBCO film, which is too thin to sustain the quality of bulk [37,42]. Moreover, contrary to the theory (Figure 20a,b), there was no significant difference in the real part of the refractive index from 20 to 60 K throughout the experiment (Figure 20c,d). The pulse width of the probe beam was approximately 10 fs; therefore, it was difficult to distinguish the time delay when the optical path was less than 3 µm. As the thickness of the thin film is 70 nm and the real part of the refractive index is approximately 50, the optical path is approximately 3.5 µm, which is close to the resolution limit. Therefore, it is difficult to accurately estimate a refractive index value that is smaller than this. However, as previously mentioned, the imaginary part of the refractive index corresponds more with amplitude, so it would be qualitatively good agreement with the results obtained by the two-fluid model.
In the experiment, the SNR is 5000:1. Therefore, it hardly affects the error in the estimated values, and the most critical error is the uncertainty in the thickness of the GdBCO film, as shown in the SEM image in Figure 2. Accordingly, the refractive index and the related error bars that can occur due to the thickness at the highest and lowest temperatures are shown in Figure 21. The error of the refractive indices due to the thickness uncertainty of about ±10% is about ±5%.
In summary, we experimentally estimated the refractive indices of a 70 nm-thick superconducting GdBCO film via transmission-type terahertz time-domain spectroscopy combined with the numerical extraction method. Although many studies were performed with superconducting thin films, especially in the terahertz and lower frequency regions, they were limited to analysis of the actual index of refraction of the superconducting thin film, where its thickness is in the sub-skin depth region. Our terahertz time-domain spectroscopy, based on the 10-fs pulse laser, showed a great possibility of estimating the terahertz dielectric constants of sub-100 nm thick superconducting films, even in the superconducting state. In addition, this technique, combined with the numerical extraction method, can be applied to experimentally analyze the terahertz optical properties of various thin film materials with the one-shot measurement. In summary, we experimentally estimated the refractive indices of a 70 nm-thick superconducting GdBCO film via transmission-type terahertz time-domain spectroscopy combined with the numerical extraction method. Although many studies were performed with superconducting thin films, especially in the terahertz and lower frequency regions, they were limited to analysis of the actual index of refraction of the superconducting thin film, where its thickness is in the sub-skin depth region. Our terahertz time-domain spectroscopy, based on the 10-fs pulse laser, showed a great possibility of estimating the terahertz dielectric constants of sub-100 nm thick superconducting films, even in the superconducting state. In addition, this technique, combined with the numerical extraction method, can be applied to experimentally analyze the terahertz optical properties of various thin film materials with the one-shot measurement.  Acknowledgments: Part of this study has been performed using facilities at IBS center for Correlated Electron Systems, Seoul National University.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Data Availability Statement:
The data presented in this study are available from the corresponding author.