www.nature.com/scientificreports OPEN High frequency response of thick REBCO coated conductors in the framework of the FCC study Artur Romanov1*, Patrick Krkotić2,3, Guilherme Telles1, Joan O’Callaghan3, Montse Pont2, Francis Perez2, Xavier Granados1, Sergio Calatroni4, Teresa Puig1* & Joffre Gutierrez1 A thorough microwave response study of high temperature superconductors, considered as an alternative beam screen coating, has become integral in the design decisions for CERN’s future research infrastructure. Here, we present the surface resistance Rs of various REBa2 Cu3 O7−x (RE = rare earth) coated conductors available in large scale as a function of magnetic field in a broad temperature range measured by a Hakki–Coleman type resonator with resonant frequency ν ≈ 8 GHz. Analysis of the high frequency dissipation supported by DC transport characterization reveals the vortex dynamics in thick REBa2 Cu3 O7−x films. Determined microscopic vortex parameters span over a wide range of magnitudes and reflect the relevance of the superconducting layer’s microstructure. We demonstrate that the depinning frequencies ν0 surpass ν , which confirms the operation in high performing, low dissipation pinning regime at measurement conditions. Surface impedance extrapolation to FCC-hh conditions points towards a significant outperformance of copper by coated conductors in terms of surface resistance. The surface resistance margins would open up potential ways for a more efficient frontier circular collider. Shortly after the discovery of high temperature superconductivity (HTS) in oxides by Bednorz and Mueller in ­19861, ceramic cuprates, in particular YBa2 Cu3 O7−x (YBCO), triggered universal euphoria among scientists and industry for the unprecedented high critical temperatures and sound performances under elevated magnetic fields. Decades of research activities later, second generation REBa2 Cu3 O7−x (REBCO, RE = Rare Earth), coated conductors (CCs) mark now the pinnacle in the industrialization of superconducting cuprates. The need for low grain boundary angles in these superconductors is accommodated by biaxially textured growth of ∼ µm thick REBCO on functional oxide buffer layers backed by metallic tapes. It is a quite sophisticated and complex structure with practical use but expensive in ­production2.Most prominent applications take advantage of the well studied direct current (DC) or low frequency alternating current performance profile that identifies CCs as high current conductors at cryogenic temperatures and/or high ­fields3–5. The future circular hadron-hadron collider (FCC-hh), however, considers a novel utilization of superconductors reliant on the high frequency response of HTS at extremely challenging ­conditions6. The FCC-hh is the most ambitious scenario for a post Large Hadron Collider (LHC) machine proposed by CERN. It will operate as an 80−100 km acceleration ring where 16−18 T magnets will steer 1011 protons in 8 cm bunches with revolution period of 0.3 ms. In this machine, collision ener- gies of 100 TeV at center of mass are expected to be reached. The 35.4 W/m/beam synchrotron radiation emitted by the protons will be absorbed by a stainless steel tube, the so called beam screen, held at a temperature window of 40−60 K . The beam screen chamber is located in the beam-steering dipoles to thermally shield them. Given the prospected beam parameters, mirror charges will be induced into the beam screen that peak at magnitudes of 25 A with frequencies up to 1−2 GHz . Governed by the necessity to maintain a stable beam during operation, highly conductive coatings will have to cover the interior of the beam screen chamber. Compared to copper, which has been the subject of a detailed s­ tudy7, HTS promise with their lower surface resistance larger beam stability margins at expected operation t­ emperature8–10. This motivates the need to understand and evaluate the high frequency response of HTS available on an industrial scale at FCC conditions. Our previous studies of six 1 Institut de Ciència de Materials de Barcelona, C.S.I.C., Campus U.A. Barcelona, 08193 Bellaterra, Catalonia, Spain. 2ALBA Synchrotron—CELLS, Carrer de la Llum 2‑26, 08290 Cerdanyola del Vallés, Barcelona, Spain. 3Universitat Politècnica de Catalunya, CommSensLab. c/ Jordi Girona 1, 08034, Barcelona, Catalonia, Spain. 4European Organization for Nuclear Research (CERN), 1211, Geneva 23, Switzerland. *email: aromanov@ icmab.es; teresa.puig@icmab.es Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 1 Vol.:(0123456789) www.nature.com/scientificreports/ commercially available CCs have shown lower Rs than 300 µm thick Cu colaminated on stainless steel (equiva- lent of FCC beam screen) at T = 50 K and ν = 8 GHz up to 9 T . In addition, we demonstrated compatibility with thin a-C layers to mitigate the secondary electron ­yield11. In this work, we extend the temperature range of measured surface resistances to T = 20−70 K . It allows a discussion about the surface resistance and vortex physics of CCs with different microstructures in a wide range of temperatures. The vortex parameters depin- ning frequency ν0, vortex viscosity η and Labusch parameter kp are derived within the Gittleman–Rosenblum and Bardeen–Stephen models. The model proposed by Gittleman and ­Rosenblum12 is a mean-field theory for vortices in a periodic pinning potential driven by high frequency oscillating, subcritical currents without thermal activation. It modulates both the resistive and reactive response of vortices to the driving field. In our case, the lack of complete surface reactance data sets has to be compensated by additional model confinements. With a modified Bardeen–Stephen ­model13, which describes the motion of vortices in a type II superconductor, we can estimate the vortex viscosity η. It reduces the Gittleman–Rosenblum model parameters to only one, the depinning frequency ν0, and thus makes it available through fitting the measured surface resistance Rs . We validate the two model approach for one sample with the ratio of surface reactance and surface resistance at T = 20 K , which lets us determine the depinning frequency sticking exclusively to the Gittleman–Rosenblum model. Finally, having established the microwave vortex parameters, we extrapolate the surface impedance of CCs down to 1 GHz , up to 16 T and compare it to the microwave response of Cu at FCC-hh conditions. Methods Highly sensitive surface impedance measurements are performed with a high frequency, non-destructive, dielec- tric loaded resonator (DR) with resonant frequency ν ≈ 8 GHz . The resonator is based on the Hakki–Coleman ­geometry14, was previously used for non-linear characterization of HTS thin ­films15,16 and has been redesigned as described i­ n11,17 in order to meet the sensitivity and geometrical requirements imposed by characterization of superconductors at cryogenic temperatures (5−120 K) and up to 9 T in the 25 mm bore of a PPMS from Quantum Design. The DR consists of a radiation shielding, cylindrical brass cavity terminated axially by interchangeable 12 mm × 12 mm samples and loaded with a single crystal rutile (TiO2) disk (h = 3 mm and ∅ = 4 mm), which is placed in the center of the cavity. The sample-dielectric-sample structure is sandwiched by copper-beryllium springs required to compensate contractions and expansions in the wide temperature range of measurements. The choice of dielectric material and cavity geometry results in a resonant frequency ν ≈ 8 GHz . Coupling loops through lateral walls are adjusted in such way that the resonator operates in the TE011 mode. This drives micro- wave currents to flow within the superconducting CuO2 planes of the samples under investigation. Estimations of the RF field in the cavity yields Hmax = 28.68 A/m or 0.036 mT18. An RF power dependence study showed negligible (below 2%) increments in Rs with respect to its zero-RF-field values at 0 T external applied magnetic field (see supplementary notes). External magnetic fields are applied parallel to the c-axis of the superconductor. The unloaded quality factor Q0 of the DR relates to the losses in the following ­way17,19,20: 1 RCC,1 RCC,2 Rm = + + + p tan δ. (1) Q0 GCC,1 GCC,2 Gm The first three terms of the sum are losses caused by the exposed conductive surfaces within the cavity, where the subscripts ‘CC,1’, ‘CC,2’ and ‘m’ stand for the coated conductors terminating the cavity and metallic enclosure, respectively. Here, R denotes surface resistances and G are the corresponding geometrical factors. The last term accounts for losses within the dielectric body. The dimensionless filling factor p ≈ 1, since nearly the entire energy stored in the resonator is focused in the dielectric, and tan δ is the loss tangent of rutile. In our geometry, losses of the metallic housing can be n­ eglected17 and GRCC,1 CC,1 ≈ GRCC,2 CC,2 = RGCC , thus   G 1 RCC (T, B) = · − tan δ(T) . (2) 2 Q0 (T, B) The procedures for numerical and empirical determination of G are described i­ n17, whereas tan(δ)(T) is adopted ­from21. In analogy to the surface resistance, the imaginary part of the surface impedance Zs = Rs + iXs can be extracted from microwave measurements. The reactive response of a resonator reflects in the resonant frequency ν according t­ o20,22: �ν �XCC,1 �XCC,2 �Xm �ǫ ′ −2 = + + +p ′ r , (3) νref GCC,1 GCC,2 Gm ǫr,ref where X is the surface reactance, ǫr′ the real part of the relative dielectric permittivity. By x we mean the dif- ference of a quantity from its reference value xref . Using the same approximations made already for the surface resistance, we get �ν(�T, �B) p �ǫr′ (�T)   �XCC (T, �T, B, �B) = −G + . (4) νref (T, B) 2 ǫr,ref (T) Magnetic field measurements at fixed temperatures are executed both in zero-field and field-cooled mode. Demagnetization of the superconducting external field magnet is followed after each field sweep to minimize the effect of trapped fields. DC transport characterization is done either in van der Pauw geometry on 4 mm × 4 mm samples, or with well defined track bars of 30 µm width and 200 µm length. The width of the track bar ensures a cross section Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 2 Vol:.(1234567890) www.nature.com/scientificreports/ large enough to contain a representative number of grain boundaries. The tracks are defined by optical lithog- raphy. Temperature and magnetic field control are achieved, as described above, with a Quantum Design PPMS. Subject of our study are seven commercially available CCs that differ in their microstructure. They are provided by Bruker HTS GmbH, Fujikura Ltd., SuNAM Co. Ltd., SuperOx, SuperPower Inc. and THEVA Dünnschichttech- nik GmbH. An overview of all but one CCs’ microstructure, film thickness and growth method is given i­ n11. In additional to the CCs ­in11, a new CC is provided by Fujikura Ltd. with nominal Ic − width(50 K) = 2 kA/cm − w and film thickness t = 2.5 µm, where Eu is used as the rare-earth with BaHfO3 (BHO) nano rods in the REBCO matrix. All RF characterization is done on CCs whose Ag and Cu protection layers are stripped off. In order to avoid degradation of the superconductor, the metallic layers are either chemically etched right before measure- ments or the CCs are provided stripped by the manufacturer and stored in a low humidity environment. The error on measured surface resistances arising from deviations in rutile centering and inhomogeneities in CC tapes is estimated to be 10 % as discussed ­in11, while the uncertainty in the surface reactance is, in a first approximation, assumed to result from the the error on resonant frequency σν = ν · 1.5 ppm. In the following, CCs are regarded without their protection layers and named according to their manufacturer. Results Surface resistance of thick CCs. In a general case, when exposed to a propagating microwave field, CCs’ entire complex architecture consisting of a thick superconducting layer with tSC ≥ 0.9 µm on top of a buffer layer stack with tdielectric > 0.5 µm backed by a metallic, flexible substrate with tm ∼ 50 µm contributes to the microwave response. tlayer stands here for the thickness of specified layers. Following transmission line analogy and impedance ­transformation23, the effective measured surface impedance of the CCs ZCC can be written as Zstack + iZSC tan(kSC tSC ) ZCC = ZSC , (5) ZSC + iZstack tan(kSC tSC ) with (6)  ZSC = i2πνµ0 ρ̃ being the surface impedance of the superconductor (valid in the local limit) and ρ̃ the complex ­resistivity24,25. Zstack is the effective measured surface impedance of the layers below the superconductor (dielectric buffer lay- ers/metallic substrate) and kSC = µZ0 2πν SC the propagation constant of the superconductor. So, when interested in the characteristic surface impedance of the superconductor ZSC , knowledge about the microwave response of the entire multilayer stack below the superconductor Zstack is necessary. However, Eq. (5) can be simpli- fied if considering the characteristic screening lengths in the material under study. In the framework of the two fluid model, |kSC tSC | ≫ 1 implies ZCC ≈ ZSC , where the influence of the layers below the superconductor can be ­ignored26. This is known as the thick film approximation. |kSC tSC | can be estimated for samples in our experiments. By assuming Tc = 93 K , ρn (70 K) = 60 µ� cm and l (0 K) = 150 nm , typical values for YBCO, and T = 70 K , ν = 8 GHz , as measurement conditions where the smallest propagation constant is expected, we get for our thinnest and thickest CC (tmin ≈ 0.9 µm, tmax ≈ 3.0 µm ) |kmin tSC | ≈ 5 > 1 and |kmax tSC | ≈ 16 > 1, respectively, which does not clarify unambiguously if the thick film approximation is applicable. However, Silva et al.26 provide a numerical example for the deviation of the thick film approximation from Eq. (5). For 1 µm YBCO on a LAO single crystal substrate at ν = 30 GHz and T = 72 K , the error on the characteristic surface resistance by ignoring the substrate RSC+LAO ′ = RSC is smaller than 5%. This error becomes even smaller when decreasing ν . In conclusion, by applying the thick film approximation for the samples in this work, we make a systematic error on the characteristic surface resistance of the superconductor and all values derived from it that is smaller than 5%. With regard to specific microwave applications like the FCC-hh, however, the ZCC is still a very relevant quantity as it determines the effective microwave response of CCs in use. We consider the expression for the experimental surface resistance as a function of RF field Hrf , externally applied flux density B and temperature T8: Rs (Hrf , B, T) = RBCS (Hrf , 0, T) + Rres (Hrf , 0, 0) + Rfl (Hrf , B, T) (7) with BCS formalism deduced component RBCS and the, in general, field and temperature independent resid- ual surface resistance Rres resulting from impurities and defects. RBCS usually considers only the temperature ­dependence27. The externally applied magnetic field dependence enters through Rfl which represents the losses induced by the movement of fluxons when rendered to the mixed ­state8,28. Fluxons or vortices are understood as quantized magnetic field tubes with flux 0 = 2e h and a normal conducting core. For cuprates the vortex cores are of radius ξ ∼ nm, known as the coherence l­ ength29. The last term is subject of many t­ heories12,30 and will be later integral part of our analysis. In Fig. 1 we show Rs of three selected coated conductors as function of applied magnetic field µ0 H at several cryogenic temperatures T. SuNAM (as SuperOx and THEVA) represents a CC whose microstructure relies on an as grown pinning landscape with twin boundaries, stacking faults and point defects (pristine CC)31. Coated conductors with ’nanoengineered’ microstructures, meaning that secondary phases are deliberately incorporated into the REBCO matrix to enhance the pinning strength, are provided by Bruker, Superpower and Fujikura. Typically, nanoengineering can be realized in two ways. Either by adding artificial pinning centers (APC) made from nanosized oxides, e.g. BaZrO3 (BZO) or BHO, in production processes or by establishing specific growth conditions that promote defects. The corresponding CC microstructures control the field behavior of the surface resistance. The SuNAM sample has a steeper increase in Rs with magnetic field B ≈ µ0 H especially at mid and high T as compared to Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 3 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 1.  Absolute surface resistances Rs of three selected coated conductors as a function of magnetic field at T = 20−70 K. The markers (squares, triangles, circles) represent measurement points. Rs (m ) 0 10 Bruker Fujikura Fujikura APC SuNAM SuperOx SuperPower Theva -1 10 0 0.5 1 1.5 2 2.5 3 3.5 4 0 H (T) Figure 2.  Rs (H) of CCs measured at T = 50 K and ν ≈ 8 GHz. The solid lines guide the zero-field cooled measurements, while the dashed lines connect the values recorded with field cooling. Below 2 T, a hysteretic behavior can be observed that points towards influences of trapped fluxons. SuperPower or Bruker. The higher in-field performance of APC based superconductors known from critical current ­measurements32,33 is encountered here, too. This points towards an influence of APC on the surface resistance, hence, vortex dynamics where pinning is meaningful. Comparing Bruker with SuperPower, the field behavior at low and intermediate temperatures up to 50 K is fairly similar. At higher T a slight difference in Rs emerges, visible particularly at 9 T . It is in accordance with Bruker’s ’double disordered’ coined microstruc- ture which is designed to pin vortices strongly at 4.2 K and high fields at expense of performance at higher ­temperatures34. The features of the surface resistance’s field behavior that categorize the CCs into pristine and nanoengineered microstructure are observed for all measured samples. An overview of surface maps for all CCs can be found in the supporting information. Figure 2 shows Rs as a function of magnetic field at T = 50 K for all considered CCs. A hysteretic behavior associated to field cooled/zero field cooled measurement procedures occurs at fields < 2 T, which indicates an Rfc (B=0T)−Rszfc (B=0T) influence of trapped vortices. In Table 1, we present s Rzfc (B=0T) for all measured samples. Interestingly, s the hysteresis appears in different magnitudes for each CCs ranging between 4% and 44%. A study about the effect of trapped field on the surface resistance of CCs and its relation to the superconducting layer’s architecture remains to be conducted. It might allow an optimization of CC technology targeting an improved field quality within the FCC-hh. Evaluation of microscopic vortex parameters. As described by Eq. (7) the microwave response with applied magnetic field is driven by the movement of fluxons in the superconductor. Consider a single vortex Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 4 Vol:.(1234567890) www.nature.com/scientificreports/ Rsfc (B=0T)−Rszfc (B=0T) Provider Rszfc (B=0T) Bruker 44% Fujikura 31% Fujikura APC 15% SuNAM 35% SuperOx 27% SuperPower 4% Theva 10% Table 1.  Surface resistance hysteresis after exposing to 9 T for all CCs at 50 K. from an externally applied magnetic field B = n0, n is the number of fluxons per unit area, in a type-II super- conductor under a harmonic transport current density J = J0 exp (i2πνt) na,b with J0 not close to the critical current density Jc . Here, ν is the current frequency in the range of 1−100 GHz and na,b a unit vector perpendicu- lar to the applied magnetic field with unit vector nc . The vortex will be exposed to a Lorentz force FL = J × 0 nc exerted by the interplay with the electrical current. Under the presence of pinning sites, the vortex motion is constrained by pinning potential wells. This results in a pinning force Fp = −kp x which can be assumed to be linear to its vortex displacement x . kp is the pinning constant or Labusch parameter encoding both the strength of pinning sites and vortex e­ lasticity25. Additionally, the moving normal conducting vortex core gives rise to a viscous force Fv = −ηẋ , with η being the vortex viscosity. It can be understood either as a conversion of Cooper pairs to quasiparticles where energy is l­ost29 or Joule dissipation within the vortex c­ ore13. By neglecting fluxon mass and thermal activation, the equation of motion for a single vortex i­ s24,25 ηẋ + kp x = J × �0 nc . (8) This equilibrium of forces is valid for not too large vortex displacements x in relation to the size of the potential wells within the superconductor. Moreover, nonlinearities are not considered in this ­model24. Gittleman and Rosenblum were the first to derive a response function to equation (8) in form of the complex vortex ­resistivity12 1 B�0 1 ρv = ρff = . (9) 1 + iν0 /ν η 1 + iν0 /ν k ν0 = 2πηp denotes the depinning frequency, that is the material characteristic frequency above which the vortex response changes from the pinning to the flux-flow regime. The former is dominated by reactive, the latter by resistive contributions. At ν ≪ ν0 the vortex resistivity is completely imaginary with ρv = −i2πνB�0 /kp. ν ≫ ν0 reaches the extreme of the flux-flow regime in which ρv = ρff . From Eqs. (6) and (9), we can express the surface impedance originating from fluxon movement as  Zfl = Z(T, B) − Z(T, 0) = �Z = i2πνµ0 ρv . (10) As shown in Figs. 1 and 2, it becomes evident that the vortex response, namely the field dependent part of the surface resistance dominates the total surface resistance under sufficiently high applied magnetic field. Hence, a good grasp of a superconductor’s in-field performance under the influence of microwaves can be gained by estimating the vortex properties η and ν0, that are function parameters of Eq. (9). We start with determining the vortex viscosity η. Bardeen and Stephen derived for the approximation of a superconductor in the ideal local limit (l > ξ0 ) an expression for the flux flow resistivity ρff . They assumed a normal conducting core of finite size with radius ξ0 ∼ √B1 . The flux flow resistivity is governed then by the c2 normal state resistivity ρn and the vortex core size, thus ρff ∼ ρn BBc2 13,29. This approximation was derived for s-wave superconductors. With the time dependent Ginzburg–Landau theory, Ivlev and Kopnin modified the flux flow resistivity for YBCO by the factor a = 1.45 to account for the layered structure of the anisotropic uni- axial superconductor in the range µ0 Hc1 ≪ B ≪ Bc235,36 : B ρff = ρn . (11) a · Bc2 Comparison of Eq. (11) with Eq. (9) yields for the vortex viscosity of REBCO superconductors �0 Bc2 η=a . (12) ρn A possible way to determine the vortex viscosity is by DC transport measurements from which the normal state resistivity ρn and upper critical field Bc2 of the superconductors can be estimated. In Fig. 3a, we show the resis- tivity as a function of temperature in the range T = 85 − 300 K at different applied magnetic fields up to 9 T for one exemplary CC. The magnitude of ρ follows values reported in l­iterature2. ρn below Tc can be characterized in a first approximation from a linear extrapolation of the metallic behavior well above the critical temperature Tc . It was demonstrated ­in37 that the free flux-flow model applies for a considerable part of the resistive transi- tion. For this reason, we can determine the magneto resistance measured close to Tc as ρ(T ≈ Tc , B) = ρff (T, B). Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 5 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 3.  By way of example, data for Fujikura without APC is shown. (a) ρ vs. T at various applied magnetic fields. The blue line is a linear fit to ρ(T) well above Tc which is used to estimate ρn. The inset zooms into the range where ρ = ρff . (b) The ratio Bρn /1.45ρff as a function of temperature, as determined by flux flow scaling. Solid red line in the figure represents the upper critical field Bc2 (T). Bruker Fujikura Fujikura APC SuNAM SuperOx -6 SuperPower 10 (Ns/m2) -7 10 10 20 30 40 50 60 70 80 T (K) Figure 4.  Vortex viscosity for all CCs but Theva determined from Eq. (12) at T = 20−70 K. The solid lines are guides for the eyes. According to Eq. (11), the envelope of Bρaρff (T) represents points where ρ(B) = ρff (B) and gives an estimate for n (T) Bc2 (T) as shown in Fig. 3b for an example CC. As proposed in the theory of 3-D GL ­superconductors29, the upper critical field follows a linear temperature dependence β(Tc − T), with β = 1.6−1.9 T/K depending on characterized CC. The slope coincides with values β = 1.9−2.1 T/K given in previous ­reports37,38. Even below T = 78 K , where the upper critical field of YBCO crosses over to 2-D behaviour with Bc2 ∼ (Tc − T)1/229, our estimates of Bc2 (20 K) = 115−135 T from the linear extrapolation yield reasonable accordance with reported rep. values Bc2 (20 K) ≈ 120 T39. Eq. (12) is fed with estimated ρn and Bc2 in order to calculate the vortex viscosity for measured tempera- tures. The values are presented in Fig. 4 for all CCs but Theva. Vortex parameters of samples from THEVA Dünnschichttechnik GmbH are not presented in this article. The inclined REBCO c-axis in Theva CCs requires the consideration of anisotropy in quasiparticle and/or vortex dynamics which goes beyond the scope of pre- sented elaborations and deserves a separate study. Error bars are derived from the determination of ρn and Bc2 and then propagated to η. The vortex viscosity is for all samples inversely proportional to the temperature. While the magnitudes of η for most of the CCs are within the errors relatively close together, the sample Fujikura with APC significantly undercuts the other CCs by a factor of up to 4 depending on the temperature. Since η relates to the conversion of cooper-pairs to quasiparticles and vice versa of a moving flux ­line40, the entire microstructure of a CC may contribute to η and explain the variety between values of different providers. The second vortex parameter ν0 in applied magnetic field can now be extracted from Eq. (10). When con- straining the Gittleman–Rosenblum model with the determined vortex viscosities (see Fig. 4), we can employ ν0 as a fitting parameter to the measured field response of the surface resistance Rfl (B) = �Rs . By way of example, Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 6 Vol:.(1234567890) www.nature.com/scientificreports/ 10 20 K 30 K 8 40 K 50 K 60 K 70 K 6 Rs (m ) 4 2 0 0 2 4 6 8 10 0 H (T) Figure 5.  The markers show the dissipation under microwaves of ν ≈ 8 GHz with applied magnetic field µ0 H up to 9 T at different temperatures of a Bruker CC. Corresponding fits with real part of equation (10) are represented by solid lines. The lines follow the same color code as the markers. Figure 6.  (a) Depinning frequency determined from fits with real part of Eq. (10) and (b) Labusch-parameter for all CCs determined from the relation kp = 2πν0 · η at T = 20−70 K in 10 K steps. Error bars for ν0 are derived from the quality of the fittings and are propagated together with uncertainty of η to the error of kp. The solid lines are guides for the eye. Rs of a Bruker sample in the range T = 20−70 K together with the real part of Eq. (10) fitted√to each data set of a corresponding temperature is shown in Fig. 5. Up to 50 K , R increases quite steadily with B as described by the Gittleman–Rosenblum model. At higher temperatures, √ especially at high fields, measurement data deviates from the predictions with a steeper increase than B . The Gittleman–Rosenblum model is limited to the mean- field regime, that is where the displacement of vortices from the pinning potential minima is small enough so that the flux tubes stay inside the wells and vortex–vortex interaction is infinite, resulting in a collective movement as a l­ attice24,25. There are two possible explanations for the observed model deviation. On the one hand, thermal activation can excite hops between pinning sites, which a mean-field model cannot account f­ or41. On the other hand, at high fields we are approaching Birr . For example, Birr (T = 70 K) ≈ 16 T for Bruker. In this case, the long-range order gets disturbed and the infinite vortex–vortex repulsion of mean-field models loses its validity. Figure 6a shows the depinning frequency ν0 extracted from fits with Gittleman–Rosenblum model for all CCs at T = 20 − 70K . Except for Fujikura with APC, the deviation of depinning frequencies as function of temperature is quite small and inversely proportional to temperature. A correlation between microstructure and depinning frequency can be observed. Bruker (double disorded), SuperPower (BZO nanorods) and Fujikura APC (BHO nanorods) characterized by the use of artificial pinning centers exhibit increased depinning frequencies as compared to SuNAM or SuperOx CCs that rely on as grown pinning. Remarkably, Fujikura without APC has within the errors the same depinning frequency as Bruker and SuperPower at all temperatures which points Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 7 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 7.  The complex vortex resistivity of the bottom surface of SuNAM as a function of applied magnetic field at T = 20 K. The inset shows the corresponding ratio s = Im(ρ Re(ρv ) = ν as a function of magnetic field. The v) ν0 errors on ρv and s are propagated from errors on the surface resistance Rs and resonant frequency ν. towards a nanoengineered microstructure that resembles the pinning properties of both CCs with APC. The determined values for all but one CC remain in accordance with literature which report mostly ν0 = 10−50 GHz in the given temperature range (­ see25,42 and references within). The depinning frequency of Fujikura with APC exceeds conventionally measured ν0 both in magnitude and steepness of the temperature dependence. As already mentioned above, Fujikura with APC is the only CC that incorporates Eu as the rare earth and BHO as the nano additive. Whether these two changes are the governing factors for the increased depinning frequencies cannot be concluded unambiguously without additional studies. In order to determine the depinning frequency ν0 within the Gittleman–Rosenblum model only, the resistive and reactive part of the surface impedance are needed. It gives rise to the ratio s = Im(ρ Re(ρv ) = ν and thus the deter- v) ν0 mination of depinning frequency ν0 independent of the vortex viscosity η. We measured the increase of surface reactance with magnetic field �Xs = Xs (T, H) − Xs (T, 0) for one CC at T = 20 K. As already discussed and seen from Eq. (3), contributions to the reactive response of a resonator can also arise from the dielectric. During the measurement of Xs, we made carefully sure to eliminate those contributions by reaching a temperature stability T with which �ǫr′ (�T) ≈ 0 (see Eq. (4)). Figure 7 presents ρv for a SuNAM CC at T = 20 K and high fields. The real part of the vortex resistivity follows, as expected, a linear field behavior. The upward curvature in Im(ρv ) suggests a field behavior in the vortex ­parameters43. From the inset of Fig. 7, the depinning frequency can be esti- mated as a function of magnetic field and is in the range ν0SuNAM (20 K, 4−9 T) ≈ 16.3−19.3 GHz . The excellent agreement with ν0SuNAM (20 K) ≈ (17.31 ± 1.25) GHz determined from a Gittleman–Rosenblum model fit with constrained η, as presented in Fig. 6a, reinforces the legitimacy of our approach for vortex viscosity estimations. Both, the obtained values for the depinning frequency as presented in Fig. 6a and the ratio s = Im(ρ v) Re(ρv ) > 1 point towards pinning based vortex dynamics of the CCs at measurement conditions. For all characterized CCs at least up to T = 60 K and 9 T, the measurement frequency ν is well below the depinning frequency ν0. Thus, we would expect the reactive part of the surface impedance to be dominant over the resistive part, which is, in fact, observed for SuNAM at T = 20 K in Fig. 7. These conditions indicate operation in the so called pinning regime where the microwave response is low in dissipation and the surface resistance is governed by the Labusch parameter kp , consequently by the pinning forces. Phenomenologically speaking, at low frequencies vortices are displaced from their potential minimum far enough to ’see’ the pinning well. It contrasts with the flux-flow regime. In the flux-flow regime, where ν > ν0, oscillations of flux tubes are so small that the pinning potential becomes invisible to the nanometric fluxon movements. The fluxons dissipate in their potential minimum as if no pinning potential was present. In this case, the vortex viscosity η mainly contributes to R25. The Labusch parameter kp = 2πν0 · η is presented as a function of temperature for all CCs in Fig. 6b. Here, the CCs’ different microstructures are reflected again quite well. Bruker, SuperPower and both Fujikura samples, defined as CCs with high pinning microstructure, exhibit over the entire T range increased kp as compared to as grown pinning CCs. As described above, the magnitude of measured surface resistances Rs are mostly gov- erned by the Labusch parameter kp at measurement conditions, meaning CCs with small kp display high surface resistances. Uncertainties in the determination of pinning parameters given in the figures originate from estimated errors on Rs , ν and goodness of fits. While they are significant in the comparison of the CCs, they do not account for the uncertainties of the intrinsic vortex parameters introduced by the pure choice of a model. In all the analy- ses, we have neglected flux creep effects which may be relevant especially at elevated temperatures. It has been reported that by neglecting the flux creep, the real part of the intrinsic microwave response is overestimated. Thus, the depinning frequency and Labusch parameter derived here within the Gittleman–Rosenblum model may represent lower ­limits43. Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 8 Vol:.(1234567890) www.nature.com/scientificreports/ 100 Zs (m ) 10-1 Bruker Fujikura 10-2 Fujikura APC SuNAM SuperOx SuperPower -3 FCC Cu 10 0 2 4 6 8 10 12 14 16 0 H (T) Figure 8.  To T = 50 K and ν = 1 GHz extrapolated surface impedance for 300 µm thick Cu sputtered on stainless steel, denoted here as FCC Cu, and several CCs as a function of magnetic field. The stars represent both the surface resistance and surface reactance values of FCC Cu. The solid lines are surface resistances, the dashed lines surface reactances of CCs. The provider of the CCs are color coded. Extrapolation of surface impedance to FCC conditions. Transverse beam coupling between acceler- ated hadrons and the introduced image currents within the beam screen coating pose one of the performance limitations in the current FCC-hh design. The transverse beam coupling impedance is directly proportional to the surface impedance of the material facing the particle b­ eam9,44. In a previous work, we compared the surface resistance of the potential FCC-hh beam screen coating Cu with the CCs presented here. At T = 50 K , µ0 H = 9 T and ν ≈ 8 GHz , the best performing CC has a 2.5 smaller surface resistance than ­Cu11. Aforemen- tioned work and previous discussions enable extrapolation of the surface impedance of Cu and CCs to FCC-hh conditions. The microwave response of 300 µm thick Cu colaminated on stainless steel by CERN is measured with the presented dielectric resonator and amounts RsFCC Cu (50 K, 0 T, 8 GHz) ≈ 7.5 m�11. Since the surface ­ epth44,45 impedance of a normal, non-magnetic metal is related to its skin d  2ρ δ= (13) 2πνµ0 by ρ Zsmetal = (1 + i) , (14) δ w e c a n d e d u c e t h e s u r f a c e r e s i s t a n c e o f C u a t F C C c o n d i t i o n s t o b e RsFCC Cu (50 K, 0 T, 1 GHz) = 18 GHz FCC Cu (50 K, 0 T, 8 GHz) ≈ 2.65 m�. The change with applied magnetic GHz · Rs field is insignificant. As indicated by Eq. (14), Rs = Xs for normal conductors. Returning to Eq. (7), the field independent part for the surface resistance RBCS of a superconductor in the dirty limit (short mean free path) has a quadratic frequency dependence for T < Tc and not too high frequen- cies ( 2πν ≪ 2�/ℏ ,  is the energy gap)46. Rres can exhibit different frequency ­dependencies47, however, in high Tc superconductors quadratic frequency behaviors were ­observed46,48. In this setting, it can be considered CC (50 K, 0 T, 1 GHz)+RCC (50 K, 0 T, 1 GHz) = 1 GHz 2 ·RCC (50 K, 0 T, 8 GHz) ≈ 10 µ�. The field independ-   RBCS res 8 GHz s ent part of the surface reactance can be written within the BCS theory as XBCS = 2πνµ0 l , where l stands for the London penetration depth. Assuming the conventional temperature d ­ ependence29 2l (0K) 2l = (15) 1 − (T/Tc )4 with l (0K) = 150 nm2 and Tc the critical temperature measured by transport measurement of corresponding CCs, we can estimate the expected surface reactance XBCS (T = 50 K, ν = 1GHz) ≈ 1.2 m�. Finally, the field dependent part of the superconductor’s surface impedance at FCC conditions is calculated from Eq. (10) and estimated pinning parameters at T = 50 K . A discussion on the limits of this approach is given in the supple- mentary notes. In Fig. 8, we present surface impedance extrapolated to FCC conditions as a function of magnetic field for FCC Cu in comparison with all CCs covered in this work. The surface resistance of the CCs at 16 T is expected to be smaller than that of Cu by a factor of 15–70 depending on the provider. The surface reactance of CCs exceeds that of Cu by a factor of 1.5–2.3. So, the use of current state of the art CCs in the beam screen chamber results in highly decreased surface resistance with the trade off of slightly increased reactive microwave Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 9 Vol.:(0123456789) www.nature.com/scientificreports/ response as compared to Cu. Realized in practice, it would mean significant benefits in terms of beam stability margins, resistive power losses and possibly a minor maximum beam current limitation within the ­FCC49. Fur- thermore, the CC technology offers room for optimization. As pointed out in the previous section, CCs under study perform at FCC conditions in the pinning regime. Hence, from an engineering point of view, it is desir- able to use superconductors with large pinning constants kp = 2πν0 · η, that is, large depinning frequency and vortex viscosity for this particular application. As shown in Fig. 6a, incorporation of artificial pinning centers have positive impacts on the depinning frequency. The vortex viscosity η can be increased by either decreasing the normal state resistivity ρn or increasing the upper critical field Bc2 as indicated by Eq. (12). Conclusions We have studied the surface resistance, Rs , and vortex properties (ν0, η and kp) at microwave frequencies of seven commercially available REBCO coated conductors as a function of magnetic field up to 9 T at ν = 8 GHz and 20−70K. We have revealed that these magnitudes correlate strongly with the microstructure of the superconduct- ing layer, elucidated the vortex dynamics at high frequencies and confirmed the high microwave performance of CCs for large scale applications. We have calculated the surface impedance of CCs at FCC conditions (T = 50 K, ν = 1 GHz and 16 T ) by means of the Gittleman–Rosenblum model. Results suggest that as compared to Cu colaminated on stainless steel, the surface resistance is expected to be a factor 15–70 smaller (depending on the provider), while the surface reactance stays in the same order of magnitude. The estimated outstanding microwave performance of CCs at FCC conditions encourage to take next steps in the feasibility study. Namely, to determine experimentally the surface impedances at 1 GHz and 16 T, study the different ways to attach tens of meter-long CCs to the beam screen chamber taking into account the mechanical properties of the stainless-steel/CC ensemble, ensure magnetic field quality inside the beam-screen chamber, and study possible hazards such as a dipole magnet quench. These factors will be examined in detail as they represent essential puzzle pieces for a final, reliable cost-benefit assessment for this solution. In summary, our manuscript highlights explicit relations between superconducting properties and micro- scopic vortex parameters to tailor CCs for specific microwave applications, and promotes further considerations for employing CCs as a beam screen coating alternative to copper for the FCC-hh. Data availability Raw data was generated at ICMAB-CSIC. Raw and derived data supporting the findings of this study are avail- able from the corresponding author T.P. on request. Received: 3 April 2020; Accepted: 16 June 2020 References 1. Bednorz, J. G. & Müller, K. A. Possible high T c superconductivity in the Ba–La–Cu–O system. Z. Phys. B Condens. Matter 64, 189–193 (1986). 2. Larbalestier, D., Gurevich, A., Feldmann, D. M. & Polyanskii, A. High-T c superconducting materials for electric power applications. In Materials For Sustainable Energy: A Collection of Peer-Reviewed Research and Review Articles from Nature Publishing Group, 311–320 (World Scientific, 2011). 3. Hahn, S. et al. 45.5-tesla direct-current magnetic field generated with a high-temperature superconducting magnet. Nature 570, 496–499 (2019). 4. Stemmle, M., Merschel, F., Noe, M. & Hobl, A. Ampacity—installation of advanced superconducting 10 kV system in city center replaces conventional 110 kV cables. In 2013 IEEE International Conference on Applied Superconductivity and Electromagnetic Devices, 323–326 (IEEE, 2013). 5. Winkler, T., Consortium, E. et al. The ecoswing project. In IOP Conference Series: Materials Science and Engineering, vol. 502, 012004 (IOP Publishing, 2019). 6. Benedikt, M., Zimmermann, F., Syphers, M., Goddard, B. & Schulte, D. FCC-hh hadron collider-parameter scenarios and staging options (2015). 7. Bellafont, I. et al. Design of the future circular hadron collider beam vacuum chamber. Phys. Rev. Accel. Beams 23, 033201 (2020). 8. Calatroni, S. HTS coatings for impedance reduction in particle accelerators: case study for the FCC at CERN. IEEE Trans. Appl. Supercond. 26, 1–4 (2016). 9. Calatroni, S. et al. Thallium-based high-temperature superconductors for beam impedance mitigation in the Future Circular Collider. Supercond. Sci. Technol.. https​://doi.org/10.1088/1361-6668/aa6bd​0 (2017). 10. Benedikt, M. et al. FCC-hh: The Hadron Collider: future circular collider conceptual design report volume 3. Eur. Phys. J. Spec. Top. 228, 755–1107. https​://doi.org/10.1140/epjst​/e2019​-90008​7-0 (2019). 11. Puig, T. et al. Coated conductor technology for the beamscreen chamber of future high energy circular colliders. Supercond. Sci. Technol. 32, 094006. https​://doi.org/10.1088/1361-6668/ab2e6​6 (2019). 12. Gittleman, J. I. & Rosenblum, B. Radio-frequency resistance in the mixed state for subcritical currents. Phys. Rev. Lett. 16, 734–736. https​://doi.org/10.1103/PhysR​evLet​t.16.734 (1966). 13. Bardeen, J. & Stephen, M. Theory of the motion of vortices in superconductors. Phys. Rev. 140, A1197 (1965). 14. Hakki, B. & Coleman, P. A dielectric resonator method of measuring inductive capacities in the millimeter range. IRE Trans. Microw. Theory Tech. 8, 402–410 (1960). 15. Mateu, J. & Collado, C. A general approach for the calculation of intermodulation distortion in cavities with superconducting endplates. Applied Phys. Lett. 82, 97–99 (2003). 16. Mateu, J., Collado, C., Menéndez, O. & O’Callaghan, J. Experiments and model of intermodulation distortion in a Rutile resonator with Y-Ba 2-Cu 3-O 7−δ endplates. J. Supercond. 16, 873–880 (2003). 17. Krkotić, P., Aguasca, A. & O’Callaghan, J. M. Small footprint evaluation of metal coatings for additive manufacturing. In 2018 48th European Microwave Conference (EuMC), 882–885 (IEEE, 2018). 18. Krkotić, P. et al. RF Characterisation of HTS-CC tapes as alternative coating for the FCC-hh beam screen (2019). FCC week 2019. 19. Shen, Z. Y. et al. High T c superconductor-sapphire microwave resonator with extremely High Q-values up to 90 K. IEEE Trans. Microw. Theory Tech. 40, 2424–2432. https​://doi.org/10.1109/22.17991​2 (1992). Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 10 Vol:.(1234567890) www.nature.com/scientificreports/ 20. Alimenti, A., Torokhtii, K., Silva, E. & Pompeo, N. Challenging microwave resonant measurement techniques for conducting material characterization. Meas. Sci. Technol.. https​://doi.org/10.1088/1361-6501/ab0e6​5 (2019). 21. Klein, N. et al. Dielectric properties of rutile and its use in high temperature superconducting resonators. J. Appl. Phys. 78, 6683–6686. https​://doi.org/10.1063/1.36049​0 (1995). 22. Pompeo, N., Torokhtii, K. & Silva, E. Dielectric resonators for the measurements of the surface impedance of superconducting films. Meas. Sci. Rev. 14, 164–170. https​://doi.org/10.2478/msr-2014-0022 (2014). 23. Pompeo, N. et al. Vortex pinning and flux flow microwave studies of coated conductors. IEEE Transactions on Applied Supercon- ductivity (2019). 24. Silva, E., Pompeo, N. & Dobrovolskiy, O. V. Vortices at microwave frequencies. Superconductors at the Nanoscale: From Basic Research to Applications 253–278 (2017). https​://doi.org/10.1515/97831​10456​806-009. 25. Golosovsky, M., Tsindlekht, M. & Davidov, D. High-frequency vortex dynamics in. Supercond. Sci. Technol. 9, 1–15. https​://doi. org/10.1088/0953-2048/9/1/001 (1996). 26. Silva, E., Lanucara, M. & Marcon, R. The effective surface resistance of superconductor/dielectric/metal structures. Supercond. Sci. Technol. 9, 934–941. https​://doi.org/10.1088/0953-2048/9/11/003 (1996). 27. Hein, M. High-temperature-superconductor thin films at microwave frequencies Vol. 155 (Springer, Berlin, 1999). 28. Benvenuti, C. et al. Study of the surface resistance of superconducting niobium films at 1.5 GHz. Phys. C: Supercond. 316, 153–188 (1999). 29. Tinkham, M. Introduction to Superconductivity (Courier Corporation, Mineola, 2004). 30. Coffey, M. W. & Clem, J. R. Unified theory of effects of vortex pinning and flux creep upon the rf surface impedance of type-II superconductors. Phys. Rev. Lett. 67, 386–389. https​://doi.org/10.1103/PhysR​evLet​t.67.386 (1991). 31. Palau, A. et al. Disentangling vortex pinning landscape in chemical solution deposited superconducting YBa 2 Cu 3 O 7−x films and nanocomposites. Supercond. Sci. Technol. 31, 034004 (2018). 32. Gutiérrez, J. et al. Strong isotropic flux pinning in solution-derived YBa 2 Cu 3 O 7−x nanocomposite superconductor films. Nat. Mater. 6, 367–373. https​://doi.org/10.1038/nmat1​893 (2007). 33. MacManus-Driscoll, J. et al. Strongly enhanced current densities in superconducting coated conductors of YBa 2 Cu 3 O 7−x + BaZrO 3. Nat. Mater. 3, 439 (2004). 34. Usoskin, A. et al. Double-disordered HTS-coated conductors and their assemblies aimed for ultra-high fields: large area tapes. IEEE Trans. Appl. Supercond. 28, 1–6 (2018). 35. Ivlev, B. & Kopnin, N. Flux flow in layered high-T c superconductors. Phys. Rev. B 42, 10052 (1990). 36. Ivlev, B. & Kopnin, N. Flux-flow conductivity in anisotropic and layered high-T c superconductors. EPL (Europhys. Lett.) 15, 349 (1991). 37. Langan, R. et al. Scaling behaviour of magneto-resistance in clean YBa 2 Cu 3 O 7−δ single crystals in the vortex liquid state. Phys. C: Supercond. 313, 294–304 (1999). 38. Welp, U., Kwok, W., Crabtree, G., Vandervoort, K. & Liu, J. Magnetic measurements of the upper critical field of YBa 2 Cu 3 O 7−δ single crystals. Phys. Rev. Lett. 62, 1908 (1989). 39. Sekitani, T., Miura, N., Ikeda, S., Matsuda, Y. & Shiohara, Y. Upper critical field for optimally-doped YBa 2 Cu 3 O 7−δ. Phys. B: Condens. Matter 346, 319–324 (2004). 40. Abrikosov, A. A. Fundamentals of the theory of metals (Courier Dover Publications, Mineola, 2017). 41. Wu, D. H., Booth, J. & Anlage, S. M. Frequency and field variation of vortex dynamics in YBa 2 Cu 3 O 7−δ. Phys. Rev. Lett. 75, 525 (1995). 42. Tsuchiya, Y. et al. Estimation of vortex viscosity from the complex surface impedance measurement in the mixed state of YBa2C- u3Oy. Phys. C: Supercond. Appl. 362, 273–276. https​://doi.org/10.1016/S0921​-4534(01)00686​-4 (2001) (0102021). 43. Pompeo, N. & Silva, E. Reliable determination of vortex parameters from measurements of the microwave complex resistivity. Phys. Rev. B Condens. Matter Mater. Phys. 78, 1–10. https​://doi.org/10.1103/PhysR​evB.78.09450​3 (2008). 44. Krkotić, P., Niedermayer, U. & Boine-Frankenheim, O. High-temperature superconductor coating for coupling impedance reduc- tion in the FCC-hh beam screen. Nuclear Instrum. Methods Phys. Res. Sect. A: Accel Spectrom. Detect. Assoc. Equip. 895, 56–61 (2018). 45. Collin, R. Foundations for Microwave Engineering. IEEE Press Series on Electromagnetic Wave Theory (Wiley, 2001). 46. Attanasio, C., Maritato, L. & Vaglio, R. Residual surface resistance of polycrystalline superconductors. Phys. Rev. B 43, 6128 (1991). 47. Halbritter, J. On RF residual losses in superconducting cavities (1984). 48. Delayen, J. & Bohn, C. Temperature, frequency, and rf field dependence of the surface resistance of polycrystalline YBa 2 Cu 3 O 7−x . Phys. Rev. B 40, 5151 (1989). 49. Calatroni, S. & Vaglio, R. Surface resistance of superconductors in the presence of a DC magnetic field: frequency and field intensity limits. IEEE Trans. Appl. Supercond. 27, 1–6. https​://doi.org/10.1109/TASC.2017.26916​04 (2017). Acknowledgements Authors acknowledge the support and samples provided by Bruker HTS GmbH, Fujikura Ltd., Sunam Co., SuperOx, SuperPower Inc. and THEVA Dünnschichttechnik GmbH. Authors acknowledge CERN funding FCC- GOV-CC-0073/ 1724666/KE3359, MAT2014-51778-C2 COACHSUPENERGY, RTI2018-095853-B-C21 SuMaTe from MICINN and co-financing by the European Regional Development Fund; 2017-SGR 1519 from Generalitat de Catalunya and COST Action NANOCO-HYBRI (CA16218) from EU. Authors acknowledge Prof. Ruggero Vaglio for sharing his knowledge on vortex physics inferable from surface resistance measurements in several project meetings. ICMAB authors acknowledge the Center of Excellence award Severo Ochoa SEV-2015-0496 and its Future Interdisciplinary Projects action. UPC funding was provided through the Unit of Excellence María de Maeztu MDM2016-0600. JG acknowledges MSCA-COFUND-2014-665919 for financial support. AR acknowledges MSCA-COFUND-2016-754397 for PhD grant and would like to thank Prof. Enrico Silva, Prof. Nicola Pompeo, Dr. Kostiantyn Torokhtii and Andrea Alimenti from University Roma Tre for hosting him in their laboratory part of which were most helpful discussions on RF measurement techniques. Author Contributions J.G., T.P. and A.R. conceived the idea of the study. J.O., M.P. and P.K. designed and implemented the system to perform microwave surface impedance measurements. A.R., P.K. and G.T. conducted the experiments. A.R. carried out the data treatment. A.R., T.P and J.G. interpreted the results and wrote the manuscript with contri- butions from all authors. Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 11 Vol.:(0123456789) www.nature.com/scientificreports/ Competing interests The authors declare no competing interests. Additional information Supplementary information is available for this paper at https​://doi.org/10.1038/s4159​8-020-69004​-z. Correspondence and requests for materials should be addressed to A.R. or T.P. Reprints and permissions information is available at www.nature.com/reprints. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creat​iveco​mmons​.org/licen​ses/by/4.0/. © The Author(s) 2020 Scientific Reports | (2020) 10:12325 | https://doi.org/10.1038/s41598-020-69004-z 12 Vol:.(1234567890)