www.nature.com/scientificreports OPEN Ultra‑high critical current densities of superconducting ­YBa2Cu3O7‑δ thin films in the overdoped state A. Stangl1*, A. Palau1, G. Deutscher2, X. Obradors1 & T. Puig1* The functional properties of cuprates are strongly determined by the doping state and carrier density. We present an oxygen doping study of ­YBa2Cu3O7-δ (YBCO) thin films from underdoped to overdoped state, correlating the measured charge carrier density, nH, the hole doping, p, and the critical current density, Jc . Our results show experimental demonstration of strong increase of Jc with nH, up to Quantum Critical Point (QCP), due to an increase of the superconducting condensation energy. The ultra-high Jc achieved, 90 MA ­cm−2 at 5 K corresponds to about a fifth of the depairing current, i.e. a value among the highest ever reported in YBCO films. The overdoped regime is confirmed by a sudden increase of nH, associated to the reconstruction of the Fermi-surface at the QCP. Overdoping YBCO opens a promising route to extend the current carrying capabilities of rare-earth barium copper oxide (REBCO) coated conductors for applications. High temperature superconductors (HTS) are key materials for many applications, with special interest in those related to the energy ­sector1,2. However, to reach these opportunities, outstanding research had to be made within the last 30 years from the discovery of high temperature superconductivity, to overcome the most criti- cal materials i­ ssues3,4. The successful development of superconducting Coated Conductors (CC) is based on the uniform textured deposition of thick, homogeneous nanoengineered structures in kilometres length (a desired correlation over twelve orders of magnitude!). Cuprate materials in the form of CCs have opened up new oppor- tunities in the field of superconducting ­applications5 and nowadays offer the performance to boost technological innovations such as dissipation-free energy transmission in superconducting grids, highly efficient engines for electrical aviation or compact fusion reactors beyond ITER. However, current carrying capacities are still far from theoretical l­ imits6,7. While on one hand dissipation free current transport is intrinsically limited by the depairing current density, Jd , beyond which Cooper pair condensation is energetically not favourable, in real systems the dissipation free, critical current density, Jc is smaller due to the motion of magnetic vortices. One major goal in the development of superconducting tapes is to merge Jc with the theoretical limit, Jd . Nowadays, the main approach to increase Jc is the introduction of nanoscaled defects (artificial pinning centres) in the YBCO matrix to immobilize magnetic ­vortices8–13. The total pinning force, Fp = Jc × B, is enhanced by an increased number of elementary pinning sites of optimised size and ­distribution14. However, the normal and superconducting (SC) state properties of cuprate materials are strongly governed by hole doping of the superconducting C ­ uO2-planes. Likewise, the elementary pinning strength itself varies with the condensation energy, Ec , per coherence volume which was predicted to peak at the critical doping p* = 0.19 where the pseudogap c­ loses15. ­In16 it was reported that the mean field value of the heat capacity jump at the critical temperature, Tc , and therefore the condensation energy, increases up to the maximum oxygen stoichiometry ­O7, beyond optimal doping, in bulk YBCO. More recently, a strong increase of Ec , between the optimal doping (p = 0.16) and p* was found by measurements of the critical fi ­ elds17,18. These results further emphasize the enhancement of Jc by overdoping YBCO (as confirmed by results of the doping level), which motivated the interest of early studies to achieve overdoped YBCO ­films19–23. Nevertheless, dem- onstration of overdoping of YBCO films by oxygen doping has been quite scare up to now. In this work, we have reached the overdoped state by oxygen post-processing of YBCO thin films with a thickness of 200 and 250 nm, grown from pulsed laser deposition (PLD) and chemical solution deposition (CSD), respectively. We show the influence of oxygen doping on the charge carrier concentration determined by Hall effect measurements all the way to the overdoped state. We find that Jc strongly increases with charge carrier density, far into the overdoped state where a large Fermi surface of well-defined quasiparticles e­ xists24. 1 Institut de Ciència de Materials de Barcelona (ICMAB-CSIC) Campus de Bellaterra, 08193 Bellaterra, Barcelona, Spain. 2Department of Physics and Astronomy, Tel Aviv University, 69978 Tel Aviv, Israel. *email: alexander.stangl@grenoble-inp.fr; teresa.puig@icmab.es Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 1 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 1.  Electrical analysis of YBCO thin films. (a) Hall constant RH obtained at 3 T as a function of temperature. RH is perfectly linear within the analysed field range (− 9 to 9 T) and vanishes in zero field condition at all T as shown in the inset. Charge carrier density is obtained via Hall effect measurements using nH (T) = RH (T)q1 . (b) In-plane resistivity, ρ(T, H), as a function of temperature at different magnetic fields (0, 0.5, 1, 3, 5, 7, 9 T, H  c) to evaluate the field dependent superconducting transition temperature, T0 (H). Inset shows normalized resistivity, ρnorm = (ρ(T) − ρ0 )/bT , where b is the linear slope at high temperatures (T > 150 K). Doping dependent deviation from unity is observed, as shown for an underdoped (UD, downwards deviation), optimally doped (OpD) and overdoped (OD, upwards bending) 200 nm thick YBCO film, respectively. We demonstrate that overdoped YBCO films achieve record Jc values of 90 MA/cm2 at 5 K self-field, reaching a fifth of the depairing current density. Our results are in line with recent proposals of the controversial extrapola- tion of the pseudogap line, its relationship with the Quantum Critical Point and a Fermi Surface reconstruction ­ oint25–27, which are crucial issues to understand the nature of high-Tc superconducting cuprates, at the critical p particularly the overdoped state. We foresee the hybridization of overdoping and nanoengineering of HTS-CC as an emerging opportunity to significantly improve CC performances for a broad range of superconducting energy applications. Results Normal state electrical properties and doping state of YBCO thin films. Normal state electrical properties of cuprate superconductors strongly vary with doping. Underdoped YBCO thin films were obtained by using oxygen partial pressure below 1 bar during the post growth oxygenation process. Overdoped YBCO films were achieved by post-growth oxygen heat treatments at low temperatures, enabled by catalytically acti- vated surface oxygen exchange, using a thin Ag surface decoration layer (for details see Supplementary Infor- mation SI-Fig. 1). The catalytic effect occurs via a job sharing mechanism where Ag dissociates the molecular oxygen and YBCO incorporates the oxygen ion into the ­bulk28. Electrical resistance and Hall measurements as a function of temperature have been employed in addition to room temperature c-parameter analysis to determine the doping state of the studied films. Representative examples of the temperature dependence of the Hall resist- ance RH (T) and the electrical resistivity ρ(T, H) at different magnetic fields are given in Fig. 1 for our films. We have inferred the charge carrier density nH (T) = eRH1(T) by Hall measurements at constant magnetic field (3 T) and sweeping temperature (SI-Fig. 2). We observe a strong temperature dependence of RH, respectively nH, as previously reported for HTS c­ uprates29,30, with a maximum of RH around 100 K as shown in Fig. 1a. A similar drop of RH (T) below 100 K in underdoped YBCO was recently explained by a reconstruction of the Fermi surface in the normal state due to charge-density-wave ­order31,32. However, in the present case the downturn Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 2 Vol:.(1234567890) www.nature.com/scientificreports/ Figure 2.  Phase diagram of YBCO thin films. Zero field critical temperature as a function of charge carrier density, nH , obtained by Hall effect measurements at 100 K for thin YBCO films grown by PLD (200 nm thick, cyan diamonds) and CSD (250 nm thick, red circles). Vertical line marks optimally doping, while arrows indicate underdoped (UD) and overdoped (OD) regime. Inset magnifies Tc in the overdoped regime, showing a weak but distinct decrease with increasing charge carrier density, from above 91 K to around 89 K. in RH, which is also observed in optimally and overdoped YBCO, is caused by the onset of superconductivity, as the applied magnetic field (µHa = 3 T) is too weak to suppress the occurring phase transition to the super- conducting state. In this work, we use the charge carrier density extracted at 100 K to study the correlation with superconducting properties, such as the critical current density. The resistivity, ρ(T, H) in Fig. 1b reveals a broadening of the superconducting transition width with increasing magnetic field. From ρ(T, H) measurements we have inferred the coherence length, ξ(0), and a characteristic magnetic field, H0, as discussed below. The deviation of ρ(T) from a linear temperature dependence at low T is shown in the inset of Fig. 1b which correlates with under- and overdoping, respectively, consistent with previ- ous ­reports33. Figure 2 displays the critical temperature as a function of the Hall number, nH (100 K), measured at 100 K. Tc rapidly increases at low values of nH, as expected from the parabolic dependence of Tc on doping generally found in cuprate s­ uperconductors34 (assuming direct proportionality between the charge carrier density and doping in the underdoped region). Above the optimal doping of about nH (100 K) = 3 · 1021/cm3, Tc decreases, but following a much weaker dependence, as shown in the inset. This deviation from a parabolic dependence of Tc was reported ­previously33 and can be explained by an analysis of the doping dependence of the charge carrier density, as performed in the following. In the underdoped regime, 2 we calculated p from measurements of Tc using the parabolic dependence Tc = 82.6 p − 0.16 , with Tc,max = 92 K, generally found for cuprates. However, this method was  1 − Tc,max found insufficient for optimally and overdoped thin films. Hence, for these higher doped films we determined the doping number from measurements of the c-parameter by HR-XRD, as described in the methods (see also SI-Fig. 3a,b). In Fig. 3 we plot the evolution of the charge carrier density, n, at T = 100 K with doping, p. Here we use the charge carrier density per Cu in the ­CuO2-planes, n = nH2V , with the volume of the unit cell V (the factor ½ is owing to the fact that YBCO has two C ­ uO2-planes per unit cell). In this calculation, we assumed a constant V = 0.173 ­nm3, owing to the fact that in the range of oxygen doping explored, the variation in V is rather small (< 1%)35 while large variations in nH occur. We notice that although the same general behaviour is observed for both types of films, both reaching the overdoped state, PLD films could usually reach a higher oxygenation state than CSD. This could probably be associated to different oxygenation kinetics affecting the microstructure of the film (strain, defects, etc.). For p < 0.16, we find n = p in a broad doping range (underdoped regime). Above optimal doping (p > 0.16), n sharply increases. A similar result was reported previously, where nH was measured at low temperatures using very high magnetic fields suppressing the superconducting state (e.g. at 50 K and up to 88 T)31. This can be understood by a Fermi surface reconstruction (FSR) in proximity to the pseudogap (PG) critical point at p*, resulting in a non-unique relation between the charge carrier density, n, and doping, p, over the full range of the cuprate phase ­diagram27. In the far overdoped regime, cuprate superconductors have been recently shown to exhibit a large cylindrical Fermi Surface (FS) with n(T → 0) = 1 + p in the zero temperature ­limit31,32,36,37. Below the critical doping, the volume of the FS reduces by one hole per Cu in the ­CuO2 plane to n(T → 0) = p, e.g. due to the introduction of an antinodal gap opening which seems to arise from short range antiferromagnetic correlations in the ­system25. This transition is expected to occur within a narrow doping range between the optimal doping p = 0.16 and the closing of the PG at the critical doping p∗ = 0.19, as observed in our study. The PG critical point at p* has been shown to have all the features of a quantum critical point at T = 038. It is remarkable that the charge carrier density measured at 100 K, above the onset of superconductivity, preserves the expected behaviour for n(T → 0) of a single band metal with a Fermi-surface containing small hole-like pockets on the underdoped site and a sharp transition towards a large Fermi-volume, due to a recon- struction of the Fermi surface, in overdoped YBCO. The latter is also the reason for the observed deviation from a parabolic doping dependence of Tc on nH, shown in Fig. 2. Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 3 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 3.  Evolution of the charge carrier density with doping. YBCO normal state charge carrier density per ­CuO2-plane, n, is drawn as a function of doping p. The charge carrier density is given by n = nH2V at 100 K. The doping p is obtained for optimally and overdoped samples (full symbols) via HR-XRD measurements of the c -parameter and for underdoped films (open symbols) from the parabolic doping dependence of Tc . The vertical lines indicate optimally and critical doping. Below popt we find n = p, corresponding to a Fermi-surface with small hole and/or electron pockets in the underdoped regime. For p > p∗ a large, cylindrical FS is expected in the metallic overdoped regime, with n = 1 + p. A transition between a small and a large Fermi-surface occurs above p = 0.16. This is in good agreement with previous reports, but remarkable, as within this work n is obtained by Hall effect measurements using small fields above the onset of superconductivity at 100 K. Superconducting properties of YBCO films: Correlation with doping and condensation energy. The main result of this work is shown in Fig. 4, where we plot the inductive critical current density, Jc, as a function of nH (100 K) in (a) at self-field and (b) at an applied magnetic field of 7 T, determined by mag- netisation measurements at 5 K (SI-Fig. 4). For both cases, we find a strong increase of Jc with nH, extending beyond optimal doping, far into the overdoped regime up to the equivalent value of p* = 0.19 (corresponding nH values for optimal and critical doping are indicated with vertical dashed lines). Notice the ultrahigh values of inductive Jc at self-field and 5 K, Jc (5 K), achieved beyond optimal doping up to p*. An evaluation of the depairing current within the Ginzburg–Landau t­ heory39 is done in the following in order to compare the obtained experimental values with this theoretical limit, though strictly valid only near Tc: φ0 JdGL (T) = (1) 33/2 πµ0 2 (T)ξ (T) √ with the temperature and material dependent magnetic penetration depth (T) = (0)/ 1 − t and coherence 4 length ξ (T) = ξ (0) (1 + t 2 )/(1 − t 2 ), the flux quantum φ0 and the reduced temperature t = T/Tc. Using values for (0) as reported i­ n40 for YBCO bulk and ξ (0) as obtained within this work, we can estimate JdGL (5 K) for opti- mally and overdoped YBCO, as shown in Table 1. The ratio Jc /Jd for YBCO thin films analysed within this work, approximately doubles from ~ 9% to ~ 18% going from optimally to overdoped (p = p*), and further evidences the relevance of going to the overdoped state. The maximum Jc (5 K) obtained within this study is 89.4 MA/cm2, reaching about a fifth of the fundamental limit. To the best of the authors knowledge, this is among the highest inductive Jc (5 K) values reported for REBCO thin film superconductors at zero fi ­ eld41–45. We suggest, that this increase of Jc in YBCO films with the charge carrier density up to p* is a consequence of the increase of the condensation energy with charge carrier density, as shown in the following. Vortex matter in HTS cuprates has demonstrated to be extremely rich with new vortex phases that were not expected from the low temperature superconductivity ­knowledge46–48, specially related to the high thermal energy and flux creep ­phenomena49 that these cuprates experience. In this context, the pinning energy (which is pro- portional to the condensation energy, Ec41) can be related to an effective activation  energy  determined from m  β 50 magnetoresistance measurements, which can be written generally as U(H, T) = 1 − TTc H0 H , where β is a constant close to unity and m a material dependent parameter. The characteristic magnetic field H0 is pro- portional to the pinning energy, and thus closely related to the condensation ­energy23. We  have obtained H0 by  1/β 51–53 analysing the vortex glass transition line, given by the empirical formula HG = H0 1−t(H) t(H) , with t(H) = TG (H)/Tc , where Tc is the critical temperature at zero field and TG (H) the field dependent vortex-glass transition temperature. We find that H0 depends on the Hall number, as shown in Fig. 5a. As nH increases from 2 to 9 · 1021/cm3, H0 doubles from 35 to > 70 T. Thus, suggesting an increase of the condensation energy,Ec , of the order of a factor two within this doping range. This is in agreement with the trend observed from the analysis of the condensation energy, Ec , determined from measurements of the specific heat jump in ­Bi221254 and ­YBCO55, and the upper and lower critical fields, Hc1 and Hc2 in ­YBCO17, which can be correlated, according to BCS theory, to the zero temperature condensation energy per unit volume, Ec ∝ C and Ec ∝ Hc2. In these two cases, Ec also approximately doubles from optimal to the critical ­doping17,56. In the inset of Fig. 5a we report a Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 4 Vol:.(1234567890) www.nature.com/scientificreports/ Figure 4.  Dependence of Jc on charge carrier density: Self-field inductive critical current density, Jc , at 5 K versus charge carrier density nH (100 K) in (a) self-field and (b) an applied magnetic field of 7 T of YBCO thin films obtained by CSD (red circles, 250 nm) and PLD (cyan diamonds, 200 nm). The critical current density is determined by SQUID magnetisation measurements. Jc is strongly enhanced by increasing the charge carrier density far into the overdoped regime. Optimal and critical doping, popt and p∗, are marked with vertical lines, while the shadowed area around p∗ indicates the uncertainty of defining the critical doping in terms of nH via Fig. 3. Error bars correspond to uncertainties in film thickness. ξ(0) (nm) λ(0) (nm) Jc (5 K)(MA/cm2) Jd (5 K)(MA/cm2) Jc /Jd (5 K) Optimal doped 1.7 134 30 330 9% Overdoped 1.6 112 90 500 18% (p = p*) Table 1.  Estimation of the ratio of critical current density of this work, Jc, over depairing critical current density, Jd . ξ(0) and Jc obtained within this work, (0) ­from40 and Jd calculated using Eq. (1). linear dependence of Jc (5 K) on H0 , which reinforce the condensation energy,Ec , as the underlying quantity governing both parameters, thus further suggesting H0 ∝ E c. To verify the plausibility of our results, we analyse the expected relationship between the depairing current density, Jd , the condensation energy, Ec and the charge carrier density, nH. Therefore, we consider the condensa- tion energy per pair u(T), which is equal to the critical kinetic energy at zero temperature: 1 2 u(0) = mv , (2) 2 c with the depairing critical velocity of a pair vc and m the mass of the electron pair. The condensation energy per pair is equal to the condensation energy per unit volume divided by the superconducting pair density ns , hence u(0) = Ec (0)/ns . Using the critical velocity, we can generally write the depairing critical current density as Jd = e(2ns )vc . (3) Combining Eq. (2) and (3) we obtain Jd2 ∝ ns Ec . (4) Therefore, we expect that the square of the depairing critical current density varies as the superconducting pair density multiplied by the condensation energy. We further assume that ns varies with the measured Hall carrier density nH ∝ 2ns . Therefore, we expect that going from optimal doping to the highest achieved over-doping, nH Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 5 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 5.  Charge carrier density dependence of superconducting parameters: (a) Characteristic magnetic field, H0, versus charge carrier density at 100 K. H0 is obtained by electrical measurements of the vortex glass transition line. It provides a measure of the pinning energy and thus is closely linked to the superconducting condensation energy, Ec. The red diamond is reproduced ­from23 (280 nm, PLD), falling on the same line as our results. Inset shows Jc (5 K) versus H0, revealing a linear relation between these independently measured quantities. This emphasizes the strong correlation of the critical current density and the condensation energy. Dashed line is a guide to the eye. (b) Superconducting coherence length ξ(0) as a function of Hall number at 100 K. ξ(0) is obtained by electrical measurements of Hc2, determined from flux flow resistivity analysis up to 9 T (SI-Fig. 5b). A weak decrease with increasing charge carrier doping is observed, which alone cannot account for the strongly enhanced Jc in the overdoped region. y-error bars correspond to uncertainties from the underlying fit procedures. increases by a factor of 3, while H0 increases by a factor of 2. From the derived expression, the depairing critical current density is expected to increase over that range by a factor of 2.5, which is consistent with our results reported in Fig. 4a, suggesting that the measured self-field inductive critical current density varies similarly to the depairing critical current density. Equation 4 is experimentally verified in Fig. 6 with Jc , H0 and nH experi- mentally determined values. Thus, we believe that the increase of Jc with ultrahigh critical currents up to p* is fully consistent with an increase of the condensation energy, Ec , in YBCO films, in agreement with previous results determined from specific h­ eat56 and critical fi ­ elds17 measurements. Some of these works have suggested the existence of a peak in the condensation energy Ec , and thus in the density of states, NF, at p* when the PG closes and cuprates enter the strange metal state before reaching the Fermi liquid ­behaviour15,17,27,56. This is one of the features of a QCP irrespective of the details of the microscopic model used to describe the PG f­ ormation23,27,57,58. Unfortunately, this has been measured up to now only in Ca-doped YBCO and ­La1.8-xEu0.2SrxCuO427 and iron- pnictide ­superconductors59–62, but it could suggest a universal behaviour. It would be very interesting to be able to overdope these YBCO films beyond the present p* value to confirm also the peak in Ec. In this work, we have evaluated the self-field inductive Jc from saturated SQUID magnetization curves, to ensure being in the vortex pinning ­regime63,64, which is further confirmed by the strong increase of Jc with ns also for 7 T (Fig. 4b). Hence, variations of the coherence length might contribute to the increase of the critical current, as the pinning energy, Up (0) = Ec (0)ξ(0)3, scales with ξ (0)356 and the coherence length itself exhibits a dependence on ­doping65–67. To be able to exclude that a major contribution arises from changes in ξ (0), we have analysed ξ (0) for several samples spanning a broad doping range, as shown in Fig. 5b. ξ (0) is extracted from electrical measurements of the flux flow resistivity up to 9 T as described i­ n68–70. We observe only a small variation of ξ(0), while the Hall number changes by more than a factor 4. Thus, the Jc enhancement due to this modified coherence length would account for an upper limit of 25%, much less than the observed increase of > 300% within the doping range from 2 to 8 · 1021/cm3. Thus, we conclude that the enhancement of Jc observed is mainly governed by the modified condensation energy due to the increase of oxygen doping. Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 6 Vol:.(1234567890) www.nature.com/scientificreports/ Figure 6.  Experimental verification of Eq. (4), Jc2 ∝ nH Ec, using Ec ∝ H0. The three parameters Jc, nH and EC are derived experimentally, in this work, from independent measurements. Consequences in pinning behaviour and comparison with nanoengineered films. YBCO films are very sensitive to atomic and nanoscale defects due to the small coherence length in high temperature superconductors, which in addition leads to strong thermal fluctuations, especially at high temperatures, which strongly hinder vortex pinning. The consequence is that very different behaviour is observed for different pin- ning centres at different temperatures. While columnar defects are more advantageous at high temperatures, weak uncorrelated pinning sites strongly contribute to the overall pinning force at reduced t­emperatures41,71. The effect of random point defects is observed to be more significant at low ­temperatures10. Jc (77 K) vs Jc (5 K) is shown in Fig. 7 for films grown by PLD and CSD. In both types of samples, we observe a strong correla- tion between Jc at the two different temperatures up to about Jc (5 K) = 50 MA/cm2. However, CSD and PLD films follow different trends, which may be attributed to different pinning efficiencies associated to two dif- ferent microstructures. CSD films typically preserve a much stronger distorted matrix, giving raise to strong pinning in strained regions, being more efficient at higher t­emperatures12. Above Jc (5 K) = 50 MA/cm2, PLD films show a deviation from the initial linear relation with Jc (77 K), resulting in a much weaker dependence ( Jc (77K) increases from 3 to 4 MA/cm2, while Jc (5 K) almost doubles from 50 to 90 MA/cm2). We propose that this saturation is caused by the reduced efficiency of weak pinning sites at high temperatures. Notice, in Fig. 7 upper panel, the consistency of Tc with self-field Jc (5 K), further demonstrating the intrinsic relationship of these two magnitudes with the charge carrier density (Figs. 7, 8). The strong increase of Jc with doping motivates for a comparison of the in-field Jc (H) properties with reported values from other optimized approaches, including the more recent nanoengineering of the microstructure by the embedding of nanoparticles and nanorods in HTS films and coated conductors. Jc (H  c) at low tempera- tures is shown in Fig. 8 for various different samples (see Figure caption for details on sample composition). The most striking feature is the strong enhancement of the inductive Jc at self-field (µ0 H = 0) in the overdoped film, compared to conventional approaches by nanodefecting the matrix, as discussed above. Up to µ0 H = 2 T, overdoping enables an outstanding increase of Jc (H). Weak pinning, the expected main contribution in pristine YBCO films, is known to rapidly decrease with small magnetic fi ­ elds10, explaining the rapid decrease of Jc (H) below 1 T for the overdoped PLD sample. However, the extraordinary high self-field inductive Jc asserts high Jc (H) at even intermediate fields, comparable to the best performing nanocomposites, whose pinning is governed by 1D nanorods and strained regions around 3D nanoparticles. We also show a representative example of a pris- tine optimally doped YBCO film, as reported i­n72, highlighting the potential of overdoping YBCO. We notice that in this study, Jc was obtained by magnetisation measurements, so these results should better be compared with other inductive measurements, but we want to stress that typically, inductive measurements do suffer from extensible flux-creep effects specially at high magnetic field. Therefore, we encourage future work achieving the overdoped state in nanoengineered films to reach ultrahigh critical current densities also at high magnetic fields. Discussions We have fabricated overdoped YBCO thin films by means of different post-processing oxygen heat treatments being able to tune the doping state and reach the critical doping value p* = 0.19, close to the Quantum Critical Point. The overdoped state is confirmed by a small decrease of Tc and a transition of the normal state charge carrier density from n = p to n = 1 + p, where p is the doping state of the C ­ u2O planes, in agreement with the proposed reconstruction of the Fermi surface above optimal doping (p = 0.16). However, this result is remarkable, as n was obtained above the onset of superconductivity at 100 K, preserving the expected behaviour from the limit T → 0. We suggest that low temperature measurements at high fields, necessary to suppress superconductivity, would be highly interesting to reveal the further temperature evolution of nH (T) in these thin films. The evaluated overdoped regime is characterised by an increase of the condensation energy, leading to extraordinary self-field Jc values at low temperatures, up to one fifth of the depairing current, reaching 90 MA/ cm2 at 5 K. We experimentally demonstrate a strong increase of Jc with the charge carrier density, nH, at self-field and high magnetic fields of 7 T up to the critical point p*, and we verified the theoretical relation between Jc , Ec Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 7 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 7.  Superconducting physical properties correlations for PLD and CSD films. Tc dependence with self- field inductive Jc (5 K) (upper panel) showing the intrinsic relation of both quantities with the charge carrier density. Self-field Jc (77 K) as a function of self-field Jc (5 K) for YBCO films obtained by PLD and CSD (lower panel). The different linear trends for Jc (5 K)<50 MA/cm2 reveal growth dependent pinning defect landscapes, with possibly higher strong-pinning contribution in CSD films. With increasing Jc (5 K) in PLD films, Jc (77 K) saturates, probably due to a reduced contribution of weak pinning at high temperatures and the proximity to the superconducting transition in the overdoped regime, as indicated in the upper panel by a decreased Tc for high Jc (5 K) films. All lines are guide to the eye. Figure 8.  Jc (H  c) at 5 K of some best performing YBCO films. Field dependence of reported high critical current densities for several nanocomposite (NC) YBCO thin films in comparison with overdoped (OD) YBCO obtained in this work (PLD). Overdoped YBCO exhibits a remarkable self-field inductive Jc , almost 60% higher than previous record films, compensating the fast decrease of Jc (H) at low fields as typical for pristine YBCO films. From literature reproduced results cover the currently best-practice strategies of nanoengineering YBCO of coated conductors: 15% Zr doped (Gd,Y)BCO (Xu NC, MOCVD, at 4.2 K, reproduced under CC-BY)42, nanoscale defected REBCO with 4% BZO (Goyal NC, PLD, at 5 K, reproduced under CC-BY)78, REBCO with 15% Zr addition (Majkic NC, MOCVD, at 4.2 K, raw data was kindly provided by the author)79. Additionally we show a pristine YBCO film (Xu, MOCDV, at 4.2 K, reproduced with permission from AIP publishing) as reported ­in72. and nH (Eq. 4) from three independent experimental parameters. The distinct behaviour of Jc observed at 5 K and 77 K, suggests that we have been able to modify the weak pinning individual strength through the modifi- cation of the condensation energy by doping. p-doping strategies with oxygen post-processing treatments are expected to be scalable and uniform in long length, therefore, we envisage a viable hybridisation of overdoping and nanoengineering of YBCO films, which offers powerful prospects to further push prevailing limitations of dissipation-free current transport in cuprate superconductors and Coated Conductors at high magnetic fields of interest for applications. We reinforce the interest to find ways to overdope YBCO films beyond the present Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 8 Vol:.(1234567890) www.nature.com/scientificreports/ critical value p* to confirm the existence of a peak in condensation energy, Ec , at higher doping levels and shine light on the consequences of crossing the Quantum Critical Point. Methods Film fabrication. The ­YBa2Cu3O7-δ thin films are grown using chemical solution deposition (CSD) and pulsed laser deposition (PLD) on 5 × 5 ­mm2 ­LaAlO3 (100) and SrTi O ­ 3 (100) single crystal substrates with thick- nesses of 200 nm (PLD) and 250 nm (CSD), respectively. In case of CSD, the stoichiometric amount of precursor metal trifluoroacetate salts is dissolved in an alcoholic solution and deposited by spin-coating, followed by a pyrolysis (∼300 °C) and growth (∼800 °C) temperature treatment at PO2 = 0.2 ­mbar73,74. PLD layers are deposited at 800 °C at PO2 = 0.3 mbar with a pulse frequency of 5 ­Hz75. The PLD-targets were fabricated by Oxolutia SL (Spain) and consist of pressed and sintered, stoichiometric YBCO powder at 87% density. After growth, a 100 nm thick surface decoration layer of patterned Ag is deposited on the surface by sputtering at RT (SI-Fig. 1a). This Ag layer catalytically enhances oxygen exchange activity of YBCO during the following oxygenation process, as a job sharing mechanism facilitates the dissociation and incorporation of O ­ 2 (SI-Fig. 1b)28,76. Additionally, the silver coating provides good electrical contact for electrical measurements. Dewetting of the Ag layer into small islands with diameters of about 0.1 µm is observed at around 300 °C. However, this effect is not detrimental to its catalytic activity nor electrode functionality. Hole doping of YBCO is achieved by oxygen incorporation during the post growth oxygenation at 1 bar with an oxygen flux density of 0.16 l ­min-1 ­cm-2 at different intermediate temperatures (280–550 °C) with dwell times between 30 and 240 min. Structural characterisation. Layers obtained by either growth technique are epitaxial textured, twinned and highly c-axis oriented with no trace of secondary phases, as determined by X-ray diffraction measurements (Bruker D8 Discover), as shown in SI-Fig. 3a. In the presented films with thicknesses above 200 nm, no macro- scopic strain due to lattice mismatch with the used substrates (­ LaAlO3 and S­ rTiO3) was observed. The c-lattice parameter is obtained by HR-XRD measurements using the Nelson–Riley method, which allows the determination of the lattice parameter with very high precision, as aberration errors are minimised at very high angles (see SI-Fig. 3b). For optimally and overdoped films, we determined the doping number via the c-parameter, using the empirical equation p = c1 y + c2 y 6 + p0 with y = 1 − c/c020,77. The prefactors ci depend on the sample type and growth process. In this work we have used values reported i­ n77 (c0 = 11.8447, c1 = 11.491, c2 = 5.17 · 109 ). A small systematic constant offset was corrected by introducing the additional parameter p0 = −0.02. Magnetic and electrical analysis. SQUID magnetometry (Quantum Design) was used to determine critical current densities via the width of the magnetisation loop, as explained in SI-Fig. 4a, using the Bean criti- cal state model for thin discs. In-depth electrical analysis was performed using a Physical Property Measurement System (Quantum Design) over a broad temperature range. Contacts for electrical measurements were glued with silver paint on top of 400 µm squared Ag electrodes sputtered at the corners of the films. Electrical measurements were performed in Van der Pauw and Hall configuration in fields up to 9 T (H  c), averaged over two permutations of the electrical contacts and positive and negative excitation current in DC mode. The studied YBCO films are highly twinned, therefore the influence of metallic CuO chains short circuiting the Hall voltage along the b-direction is neglected in the calculation of the Hall coefficient RH. The critical temperature at zero field, Tc , and the field dependent vortex-glass transition temperature, TG (H), are determined from ρ(T, H) measurements to the point where the electricalresistance  in response to a small excitation current vanishes. The vortex-glass transition line 1/β HG = H0 1−t(H) t(H) is used to obtain the characteristic magnetic field, H0, by linear fitting as explained in SI- Fig. 5a. From the same measurements, we were able to determine the upper critical field, Hc2 (T), and in the zero temperature limit, Hc2 (0), using the classical Werthamer-Helfand-Hohenberg relation, which in turn defines the coherence length, ξ (0). The determination of Hc2 (T) is shown in SI-Fig. 5b. Data availability The data that support the findings of this study are available from the corresponding authors on reasonable request. Received: 22 December 2020; Accepted: 24 March 2021 References 1. Rogalla, H. & Kes, P. H. 100 Years of Superconductivity (Taylor & Francis, 2011). 2. Shiohara, Y., Taneda, T. & Yoshizumi, M. Overview of materials and power applications of coated conductors project. Jpn. J. Appl. Phys. 51, 010007 (2012). 3. Haugan, T. J., Puig, T., Matsumoto, K. & Wu, J. Artificial pinning centers in (Y, RE)-Ba-Cu-O superconductors: Recent progress and future perspective. Supercond. Sci. Technol. 33, 040301 (2020). 4. Obradors, X. et al. Nanostructured superconductors with efficient vortex pinning. In Comprehensive Nanoscience and Technology 303–349 (Elsevier, 2011). https://​doi.​org/​10.​1016/​B978-0-​12-​374396-​1.​00112-4. 5. Obradors, X. & Puig, T. Coated conductors for power applications: Materials challenges. Supercond. Sci. Technol. 27, 044003 (2014). 6. Larbalestier, D. C., Gurevich, A., Feldmann, D. M. & Polyanskii, A. High-Tc superconducting materials for electric power applica- tions. Nature 414, 368–377 (2001). 7. Civale, L. Pushing the limits for the highest critical currents in superconductors. Proc. Natl. Acad. Sci. 116, 10201–10203 (2019). Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 9 Vol.:(0123456789) www.nature.com/scientificreports/ 8. Haugan, T., Barnes, P. N., Wheeler, R., Melsenkothen, F. & Sumption, M. Addition of nanoparticle dispersions to enhance flux pinning of the YBa2Cu3O7-x superconductor. Nature 430, 867–870 (2004). 9. MacManus-Driscoll, J. L. et al. Strongly enhanced current densities in superconducting coated conductors of YBa2Cu3O7-x + BaZrO3. Nat. Mater. 3, 439–443 (2004). 10. Gutiérrez, J. et al. Strong isotropic flux pinning in solution-derived YBa2Cu3O7-x nanocomposite superconductor films. Nat. Mater. 6, 367–373 (2007). 11. Kang, S. et al. High-performance high-Tc superconducting wires. Science 311, 1911–1914 (2006). 12. Llordés, A. et al. Nanoscale strain-induced pair suppression as a vortex-pinning mechanism in high-temperature superconductors. Nat. Mater. 11, 329–336 (2012). 13. Soler, L. et al. Ultrafast transient liquid assisted growth of high current density superconducting films. Nat. Commun. 11, 2 (2020). 14. Sadovskyy, I. A., Koshelev, A. E., Kwok, W., Welp, U. & Glatz, A. Targeted evolution of pinning landscapes for large superconduct- ing critical currents. Proc. Natl. Acad. Sci. 116, 10291–10296 (2019). 15. Tallon, J. L. & Loram, J. W. Doping dependence of T* - what is the real high-Tc phase diagram?. Phys. C Supercond. its Appl. 349, 53–68 (2001). 16. Deutscher, G. New Superconductors: From Granular to High T C (World Scientific, 2006). https://​doi.​org/​10.​1142/​3432 17. Grissonnanche, G. et al. Direct measurement of the upper critical field in cuprate superconductors. Nat. Commun. 5, 1–8 (2014). 18. Talantsev, E. F. et al. Oxygen deficiency, stacking faults and calcium substitution in MOD YBCO coated conductors. IEEE Trans. Appl. Supercond. 23, 1–4 (2013). 19. Talantsev, E. F. et al. Hole doping dependence of critical current density in YBa2Cu3O7-δ conductors. Appl. Phys. Lett. 104, 242601 (2014). 20. Arpaia, R., Andersson, E., Trabaldo, E., Bauch, T. & Lombardi, F. Probing the phase diagram of cuprates with YBa2Cu3O7-δ thin films and nanowires. Phys. Rev. Mater. 2, 024804 (2018). 21. Baghdadi, R. et al. Fabricating nanogaps in YBa2Cu3O7−δ for hybrid proximity-based josephson junctions. Phys. Rev. Appl. 4, 014022 (2015). 22. Strickland, N. M., Semwal, A., Williams, G. V. M., Verebelyi, D. T. & Zhang, W. Optimizing the doping state of YBCO coated conductors. Supercond. Sci. Technol. 17, S473–S476 (2004). 23. Naqib, S. H. & Islam, R. S. Possible quantum critical behavior revealed by the critical current density of hole doped high-T c cuprates in comparison to heavy fermion superconductors. Sci. Rep. 9, 1–8 (2019). 24. Vignolle, B. et al. Quantum oscillations in an overdoped high-Tc superconductor. Nature 455, 952–955 (2008). 25. Doiron-Leyraud, N. et al. Pseudogap phase of cuprate superconductors confined by Fermi surface topology. Nat. Commun. 8, 1–7 (2017). 26. Storey, J. G. Hall effect and Fermi surface reconstruction via electron pockets in the high-Tc cuprates. EPL 113, 2 (2016). 27. Proust, C. & Taillefer, L. The remarkable underlying ground states of cuprate superconductors. Annu. Rev. Condens. Matter Phys. 10, 409–429 (2019). 28. Stangl, A. PhD Thesis on Oxygen kinetics and charge doping for high critical current YBCO films. (Universitat Autònoma de Barcelona, 2019). ISBN: 9788449087103. URL: http://​hdl.​handle.​net/​10803/​667212. 29. Ando, Y., Lavrov, A. N., Komiya, S., Segawa, K. & Sun, X. F. Mobility of the doped holes and the antiferromagnetic correlations in underdoped high-Tc cuprates. Phys. Rev. Lett. 87, 1–4 (2001). 30. Wuyts, B., Moshchalkov, V. V. & Bruynseraede, Y. Resistivity and Hall effect of metallic oxygen-deficient YBa2Cu3Ox films in the normal state. Phys. Rev. B 53, 9418–9432 (1996). 31. Badoux, S. et al. Change of carrier density at the pseudogap critical point of a cuprate superconductor. Nature 531, 210–214 (2016). 32. Collignon, C. et al. Fermi-surface transformation across the pseudogap critical point of the cuprate superconductor La1.6- xNd0.4SrxCuO4. Phys. Rev. B 95, 224517 (2017). 33. Castro, H. & Deutscher, G. Anomalous Fermi liquid behavior of overdoped high-Tc superconductors. Phys. Rev. B Condens. Matter Mater. Phys. 70, 1–8 (2004). 34. Tallon, J. L., Bernhard, C., Shaked, H., Hitterman, R. L. & Jorgensen, J. D. Generic superconducting phase behavior in high-T, cuprates: T, variation with hole concentration in YBa2C3O7-delta. Phys. Rev. B 51, 911–914 (1995). 35. Jorgensen, J. D. et al. Structural properties of oxygen-deficient YBa2Cu3O7-d. Phys. Rev. B 41, 1863 (1990). 36. He, Y. et al. Fermi surface and pseudogap evolution in a cuprate superconductor. Science 344, 608–611 (2014). 37. Fujita, K. et al. Simultaneous transitions in cuprate momentum-space topology and electronic symmetry breaking. Science 344, 612–616 (2014). 38. Tallon, J. L. et al. Critical doping in overdoped high-Tc superconductors: A quantum critical point?. Phys. Status Solidi 215, 531–540 (1999). 39. Tinkham, M. Introduction to Superconductivity (McGraw-Hill, 1996). 40. Sonier, J. E. et al. Hole-doping dependence of the magnetic penetration depth and vortex core size in YBa2Cu3Oy: Evidence for stripe correlations near 1/8 hole doping. Phys. Rev. B 76, 134518 (2007). 41. Kwok, W. et al. Vortices in high-performance high-temperature superconductors. Rep. Prog. Phys. 79, 116501 (2016). 42. Xu, A. et al. Strongly enhanced vortex pinning from 4 to 77 K in magnetic fields up to 31 T in 15 mol.% Zr-added (Gd, Y)-Ba-Cu-O superconducting tapes. APL Mater. 2, 046111 (2014). 43. Huang, J. et al. Enhanced flux pinning properties in YBa2Cu3O7−δ/(CoFe2O4)0.3(CeO2)0.7 multilayer thin films. IEEE Trans. Appl. Supercond. 25, 1–4 (2015). 44. Wang, F. & Tian, H. BaZrO3 (BZO) nanoparticles as effective pinning centers for YBa2Cu3O7 − δ (YBCO) superconducting thin films. J. Mater. Sci. Mater. Electron. 30, 4137–4143 (2019). 45. Huang, J. et al. Enhanced flux pinning properties of YBCO thin films with various pinning landscapes. IEEE Trans. Appl. Supercond. 27, 1–5 (2017). 46. Blatter, G., Feigel’Man, M. V., Geshkenbein, V. B., Larkin, A. I. & Vinokur, V. M. Vortices in high-temperature superconductors. Rev. Mod. Phys. 66, 1125–1388 (1994). 47. Crabtree, G. W. & Nelson, D. R. Vortex physics in high-temperature superconductors. Phys. Today 50, 38–45 (1997). 48. Figueras, J. et al. The loss of vortex line tension sets an upper limit to the irreversibility line in YBa 2 Cu 3 O 7. Nat. Phys. 2, 402–407 (2006). 49. Eley, S., Miura, M., Maiorov, B. & Civale, L. Universal lower limit on vortex creep in superconductors. Nat. Mater. 16, 409–413 (2017). 50. Liu, S. L., Shao, H. M. & Wang, Z. H. The activation energy U(T, H) in Er-doped melted textured growth YBa 2Cu3O7-δ crystal. Supercond. Sci. Technol. 20, 444–447 (2007). 51. Lundqvist, B., Rydh, A., Eltsev, Y., Rapp, A. & Andersson, M. Empirical scaling of the vortex glass line above 1 T for high super- conductors of varying anisotropy. Phys. Rev. B. Condens. Matter Mater. Phys. 57, 14064 (1998). 52. Rydh, A., Rapp, A. & Andersson, M. Consistent description of the vortex glass resistivity in high- Tcsuperconductors. Phys. Rev. Lett. 83, 1850–1853 (1999). 53. Andersson, M., Rydh, A. & Rapp, A. Scaling of thevortex-liquid resistivity in optimally doped and oxygen-deficient YBa2Cu3O7− δsingle crystals. Phys Rev. B - Condens. Matter Mater. Phys. 63, 1–9 (2001). Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 10 Vol:.(1234567890) www.nature.com/scientificreports/ 54. Loram, J. W., Luo, J., Cooper, J. R., Liang, W. Y. & Tallon, J. L. Evidence on the pseudogap and condensate from the electronic specific heat. J. Phys. Chem. Solids 62, 59–64 (2001). 55. Junod, A., Roulin, M., Revaz, B. & Erb, A. Experimental survey of critical fluctuations in the specific heat of high-temperature superconductors. Phys. B Condens. Matter. 280, 214–219 (2000). 56. Deutscher, G. Impact of pseudo-gap states on the pinning energy and irreversibility field of high temperature superconductors. APL Mater. 2, 096108 (2014). 57. Norman, M. R., Pines, D. & Kallin, C. The pseudogap: Friend or foe of high Tc?. Adv. Phys. 54, 715–733 (2005). 58. Daou, R. et al. Broken rotational symmetry in the pseudogap phase of a high-Tc superconductor. Nature 463, 519–522 (2010). 59. Kurth, F. et al. Unusually high critical current of clean P-doped BaFe2As2 single crystalline thin film. Appl. Phys. Lett. 106, 2 (2015). 60. Putzke, C. et al. Anomalous critical fields in quantum critical superconductors. Nat. Commun. 5, 1–6 (2014). 61. Ishida, J., Iimura, S. & Hosono, H. Effects of disorder on the intrinsically hole-doped iron-based superconductor KC a2 F e4 A s4 F2 by cobalt substitution. Phys. Rev. B 96, 1–7 (2017). 62. Hashimoto, K. et al. A sharp peak of the zero-temperature penetration depth at optimal composition in BaFe2(As1-xPx)2. Science 336, 1554–1557 (2012). 63. Talantsev, E., Crump, W. P. & Tallon, J. L. Thermodynamic parameters of single- or multi-band superconductors derived from self-field critical currents. Ann. Phys. 529, 1–18 (2017). 64. Talantsev, E. F. & Tallon, J. L. Universal self-field critical current for thin-film superconductors. Nat. Commun. 6, 2 (2015). 65. Ando, Y. & Segawa, K. Magnetoresistance of untwinned YBa2Cu3Oy single crystals in a wide range of doping: Anomalous hole- doping dependence of the coherence length. Phys. Rev. Lett. 88, 4 (2002). 66. Tallon, L., Loram, W., Cooper, R., Panagopoulos, C. & Bernhard, C. Superfluid density in cuprate high-Tc superconductors: A new paradigm. Phys. Rev. B Condens. Matter Mater. Phys. 68, 1–4 (2003). 67. Wen, H. H. et al. Hole doping dependence of the coherence length in La2-xSr xCuO4 thin films. Europhys. Lett. 64, 790–796 (2003). 68. Gordeev, S. N. et al. Resistivity of clean YBa2Cu3O7−δ in the pinned-liquid state. Phys. Rev. B 60, 10477–10483 (1999). 69. Ivlev, B. I. & Kopnin, N. B. Flux flow in layered high-Tc superconductors. Phys. Rev. B 42, 10052–10059 (1990). 70. Bartolomé, E. et al. Intrinsic anisotropy versus effective pinning anisotropy in YBa2Cu3O7 thin films and nanocomposites. Phys. Rev. B 100, 054502 (2019). 71. Blatter, G., Geshkenbein, V. B. & Koopmann, J. A. G. Weak to strong pinning crossover. Phys. Rev. Lett. 92, 3–6 (2004). 72. Xu, A., Jaroszynski, J., Kametani, F. & Larbalestier, D. Broad temperature range study of Jc and Hirr anisotropy in YBa2Cu3Ox thin films containing either Y2O3 nanoparticles or stacking faults. Appl. Phys. Lett. 106, 052603 (2015). 73. Obradors, X. et al. Progress towards all-chemical superconducting YBa2Cu3O7-coated conductors. Supercond. Sci. Technol. 19, 2 (2006). 74. Obradors, X. et al. Growth, nanostructure and vortex pinning in superconducting YBa 2 Cu 3 O 7 thin films based on trifluoro- acetate solutions. Supercond. Sci. Technol. 25, 123001 (2012). 75. Sieger, M. et al. Tailoring microstructure and superconducting properties in thick BaHfO3 and Ba2Y(Nb/Ta)O6 doped YBCO films on technical templates. IEEE Trans. Appl. Supercond. 27, 2 (2017). 76. Cayado, P. et al. Untangling surface oxygen exchange effects in YBa2Cu3O6+x thin films by electrical conductivity relaxation. Phys. Chem. Chem. Phys. 19, 14129–14140 (2017). 77. Liang, R., Bonn, D. A. & Hardy, W. N. Evaluation of CuO2 plane hole doping in YBa2 Cu3 O6+x single crystals. Phys. Rev. B Condens. Matter Mater. Phys. 73, 1–4 (2006). 78. Goyal, A. & Wee, S. H. Optimal, nanodefect configurations via strain-mediated assembly for optimized vortex-pinning in super- conducting wires from 5K–77K. J. Phys. Conf. Ser. 871, 012039 (2017). 79. Majkic, G. et al. Engineering current density over 5 kA mm2 at 42 K, 14 T in thick film REBCO tapes. Supercond. Sci. Technol. 31, 10 (2018). Acknowledgements The authors acknowledge financial support from Spanish Ministry of Economy and Competitiveness through the “Severo Ochoa” Programme for Centres of Excellence in R&D (Grant No. SEV-2015-0496), the EU ULTRASU- PERTAPE (ERC-2014-ADG-669504) and EUROTAPES project (FP7-NMP-Large-2011-280432), and the Spanish CONSOLIDER Excellence Network (Grant No. MAT2015-68994-REDC), COACHSUPENERGY project (Grant No. MAT2014-56063-C2-1-R) and SuMaTe project (RTI2018-095853-B-C21), both cofinanced by the European Regional Development Fund, and from the Catalan Government with Grant No. 2014-SGR-753 and 2017-SGR- 1519. Authors also thank the network collaboration of EU COST action NANOCOHYBRI CA16218.We also acknowledge the Scientific Services at ICMAB. We thank Juri Banchewski for some of the last transport experi- ments. A.S. is grateful for illuminating discussions on the manuscript with Stephan Steinhauer (KTH, Sweden). Author contributions A.S., A.P. and T.P designed the experimental study. A.S. performed the experimental work and analysed the data. A.S. prepared the manuscript with contributions from co-authors. All authors contributed to the scientific discussion. Competing interests The authors declare no competing interests. Additional information Supplementary Information The online version contains supplementary material available at https://​doi.​org/​ 10.​1038/​s41598-​021-​87639-4. Correspondence and requests for materials should be addressed to A.S. 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. Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 11 Vol.:(0123456789) www.nature.com/scientificreports/ 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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://​creat​iveco​mmons.​org/​licen​ses/​by/4.​0/. © The Author(s) 2021 Scientific Reports | (2021) 11:8176 | https://doi.org/10.1038/s41598-021-87639-4 12 Vol:.(1234567890)