The Depairing Current Density of a Fe(Se,Te) Crystal Evaluated in Presence of Demagnetizing Factors

: The effect of the demagnetizing factor, regarding the determination of the de-pairing current density J dep , has been studied in the case of a Fe(Se,Te) crystal, using DC magnetic measurements as a function of a magnetic ﬁeld ( H ) at different temperatures ( T ). First, the lower critical ﬁeld H c 1 ( T ) values were obtained, and the demagnetization effects acting on them were investigated after calculating the demagnetizing factor. The temperature behaviors of both the original H c 1 values and the ones obtained after considering the demagnetization effects ( H demagc 1 ) were analyzed, and the temperature dependence of the London penetration depth λ L ( T ) was obtained in both cases. In particular, the λ L ( T ) curves were ﬁtted with a power law dependence, indicating the presence of low-energy quasiparticle excitations. Furthermore, by plotting λ − 2 L as a function of T , we found that our sample behaves as a multigap superconductor, which is similar to other Fe-11 family iron-based compounds. After that, the coherence length ξ values were extracted, starting with the H c 2 (T) curve. The knowledge of λ L and ξ allowed us to determine the J dep values and to observe how they are inﬂuenced by the demagnetizing factor.


Introduction
The initial discovery of superconductivity in fluorine-doped LaFeAsO [1] attracted a lot of interest in the scientific community.Nevertheless, it was an iron-based superconductor, but with oxygen in its stoichiometric formula.This prompted us to correlate this compound with previously discovered cuprates superconductors [2].Later, with the important discovery of the Ba 1−x K x Fe 2 As 2 compound, it became clear that the importance of this new class of materials (which can present not-oxide superconductors) breaks the link with cuprates.This highlights the fact that, for fabricating high-temperature superconductors, oxygen does not necessary play a crucial role.The 11 family in the class of iron-based superconductors has been intensively investigated in recent years so that we may understand the basic mechanism governing its superconductivity [3].These materials have a layered structure, like cuprates, and they show interesting peculiarities such as high values in the upper critical field, critical current density, and irreversibility field [4][5][6][7].Moreover, they have lower anisotropy values than cuprates, with associated high pinning energy values [8][9][10][11].Among these properties, there are also non-monotonic responses to magnetization with applied magnetic fields; this has led to particular phenomena that have been deeply studied in recent years [12][13][14][15][16].These phenomena are strictly correlated with vortex dynamics [17][18][19][20][21]. On the other hand, the presence of vortices inside type-II superconductors can also be investigated by analyzing the lower critical field, H c1 , together with the London penetration depth, λ L .H c1 and λ L are useful parameters for gathering information regarding the bulk thermodynamic properties of a sample, and they have been investigated for iron-based compounds in the past [22][23][24].In particular, compared with other physical quantities, penetration depth is a useful parameter to study the superconductivity of a compound intrinsically; this is because it is not sensitive to the aspects related to surface conditions.In particular, the study of the behavior of λ −2 L , as a function of temperature, can provide information on the superconductivity typology characterizing the sample (e.g., single gap BCS theory, the two gap model, etc.) [25][26][27].Moreover, the λ L values indicate the strength of the interactions between vortices, thus allowing us to deduce the magnitude of the effective pinning energy of the sample.λ L , together with the coherence length, ξ, represent the fingerprints of a superconductor, and this study becomes crucial when first characterizing a new superconductor.Moreover, by combining them, it is possible to obtain the de-pairing current density, J dep , which fixes the upper limit of the presence of superconductivity inside a superconducting material.The de-pairing current density is of significant importance for understanding the existing limits for increasing J c [28,29], and since it directly provides data on the critical velocity of the superfluid, it is essential for the investigation of the superconducting mechanism and the symmetry of the superconducting gap [30].Using this framework, by following the Ginzburg-Landau theory, the de-pairing current density, J dep , depends on the characteristic critical parameters H c1 and H c2 , and more specifically, the London penetration depth, λ L , and the coherence length, ξ [30,31].In this work, the influence of the demagnetization effects on the de-pairing current density, J dep , has been analyzed by studying a Fe(Se,Te) iron-based superconductor.We started by measuring the first magnetization curve at different temperatures in order to obtain the lower critical field H c1 values.We have noted that the demagnetization effects acting on the sample were significant, and they resulted in an underestimation of the real H c1 values.From the H c1 values, the London penetration depth, λ L , as a function of the temperature, was obtained, and it was noted that it is not possible to fit the penetration depth with the typical exponential behavior that characterizes the s-wave superconductor.In this context, the plot of λ −2 L , as a function of T, confirmed that our sample shows peculiarities which can be ascribed to a multigap superconducting behavior.Finally, after determining the coherence length, ξ, from the upper critical field, H c2 , the J dep values were calculated as a function of the temperature, by considering the demagnetization effects and not considering them; this provided very high values in the framework of the iron-based superconductors.

Results and Discussion
In order to study the lower critical field, H c1 behavior, as a function of temperature T, the virgin magnetic moment vs. magnetic field was measured at different temperatures.In Figure 1, the first branch of the superconducting hysteresis loop was reported as reaching up to 0.3 T in the temperature range of 2.5 K to 10 K.The initial linear decrease of the magnetic moment m was visible due to the Meissner state, and the reduction of the superconducting signal was visible due to the increase in temperature.To determine H c1 at a fixed temperature, the first value that deviated from the linear trend of the Meissner state (black dashed line) meaning that the vortices have penetrated the sample building up a critical state [32][33][34].Considering all the temperatures, the H c1 (T) behavior was obtained and then fitted using the following equation [35]: where H c1 (0) is the value of H c1 at T = 0 K, T c = 14.5 K [4], and n is the exponent.As per the fitting procedure, H c1 (0) = 143 Oe and n = 1.54 were obtained.The H c1 (T) curve, together with its fit, is presented in Figure 2a (black squares and black solid line).It is worth noting that the H c1 (T) curve shows an upward trend with a negative curvature; this was not predicted in the single-band gap description of the mean-field theory, therefore, it is evident that two energy gaps exist [36].Generally, a superconductor immersed in a magnetic field is subject to demagnetization effects at low fields due to its finite dimensions [37,38].
Considering that the demagnetization effects are stronger when the sample is in a perpendicular field configuration, these effects cannot be overlooked in our case.In particular, the H c1 values are underestimated because, due to the demagnetization effects, the sample experiences a magnetic field higher than the applied one.Therefore, in order to take into account the demagnetization effects acting on the sample, the H c1 values must be properly scaled using the so-called demagnetizing factor.More specifically, the demagnetized , as a function of temperature, were calculated using the formula reported by Yeshurun et al. [39].
where N is the demagnetizing factor which can assume values ranging between 0 and 1.
The N value can be determined using the following formula [37,38]: where M is the magnetization in the Meissner state and H is the applied field.In our case, N has been already estimated to be approximately 0.76 in Ref. [4], which is quite a high value; this aligns with the perpendicular field configuration.Therefore, using Equation ( 2) with N = 0.76, the demagnetized values for H c1 were calculated and fitted with the following equation: where (T) curve, together with its fit, is presented in Figure 2a (red circles and red dashed line).where  (0) is the value of  at T = 0 K,  = 14.5 K [4], and n is the exponent.As per the fitting procedure,  (0) = 143 Oe and n = 1.54 were obtained.The  (T) curve, together with its fit, is presented in Figure 2a (black squares and black solid line).It is worth noting that the  (T) curve shows an upward trend with a negative curvature; this was not predicted in the single-band gap description of the mean-field theory, therefore, it is evident that two energy gaps exist [36].Generally, a superconductor immersed in a magnetic field is subject to demagnetization effects at low fields due to its finite dimensions [37,38].Considering that the demagnetization effects are stronger when the sample is in a perpendicular field configuration, these effects cannot be overlooked in our case.In particular, the  values are underestimated because, due to the demagnetization effects, the sample experiences a magnetic field higher than the applied one.Therefore, in order to take into account the demagnetization effects acting on the sample, the  values must be properly scaled using the so-called demagnetizing factor.More specifically, the demagnetized  values  , as a function of temperature, were calculated using the formula reported by Yeshurun et al. [39].
where N is the demagnetizing factor which can assume values ranging between 0 and 1.
The N value can be determined using the following formula [37,38]: where M is the magnetization in the Meissner state and H is the applied field.In our case, N has been already estimated to be approximately 0.76 in Ref. [4], which is quite a high value; this aligns with the perpendicular field configuration.Therefore, using Equation ( 2) with N = 0.76, the demagnetized values for  were calculated and fitted with the following equation: where  (0) is the value of  at T = 0 K,  = 14.5 K [4], and n is the exponent.From the fitting procedure,  (0) = 597 Oe and n = 1.54 were obtained.The  () curve, together with its fit, is presented in Figure 2a (red circles and red dashed line).As previously mentioned, the  values are higher than the  ones, as reported in Figure 2a.In Figure 2b, the ratio of the  () and  () values was reported to show a value equal to approximately four.On the other hand, the  (0)/ (0) value was also reported to be a red star in Figure 2b, and it is still equal to approximately four, thus confirming the viability of the fitting procedures that were previously performed.From the determination of the lower critical field values, the temperature dependence of the London penetration depth,  , can be obtained using the following formula [25]: where  = 2.07 × 10 Oe cm 2 is the magnetic flux quantum and k is the Ginzburg-Landau parameter.Using the Ginzburg-Landau theory, and by following the approach reported in Ref. [40], k can be calculated using the relation k =  2 (0)/(2 1/2  (0)), where  (0) is the thermodynamic critical field. (0) is calculated using  ,  , and  values at a temperature of zero (i.e.,  (0) = ( (0) ×  (0)) / ≈ 8 kOe and  (0) =  (0) ×  (0) / ≈ 16.5 kOe); therefore, obtaining k ≈ 40 and k demag ≈ 20 aligns with other Fe-chalcogenide superconductors [40][41][42][43].For the calculation of  , Equation ( 5), both  and  , together with the k and k demag values, were used, and the results are reported in Figure 3.Both the  () curves were fitted using the following equation: where  (0) is the London penetration depth at T = 0 K,  = 14.5 K [4], and n is the exponent.From the fitting procedure,  (0) = 208 nm, n = −0.72,and  (0) = 92 nm, and n = −0.72 were obtained, having considered the  , k, and  , k demag values for the  calculation, respectively.At low temperatures,  does not show the typical exponential behavior expected for a fully gapped clean s-wave superconductor [44].In As previously mentioned, the H demag c1 values are higher than the H c1 ones, as reported in Figure 2a.In Figure 2b, the ratio of the H demag c1 (T) and H c1 (T) values was reported to show a value equal to approximately four.On the other hand, the H demag c1 (0)/H c1 (0) value was also reported to be a red star in Figure 2b, and it is still equal to approximately four, thus confirming the viability of the fitting procedures that were previously performed.From the determination of the lower critical field values, the temperature dependence of the London penetration depth, λ L , can be obtained using the following formula [25]: where φ 0 = 2.07 × 10 −7 Oe cm 2 is the magnetic flux quantum and k is the Ginzburg-Landau parameter.Using the Ginzburg-Landau theory, and by following the approach reported in Ref. [40], k can be calculated using the relation , and H c2 values at a temperature of zero (i.e., H c (0 ≈ 16.5 kOe); therefore, obtaining k ≈ 40 and k demag ≈ 20 aligns with other Fe-chalcogenide superconductors [40][41][42][43].For the calculation of λ L , Equation ( 5), both H c1 and H demag c1 , together with the k and k demag values, were used, and the results are reported in Figure 3.Both the λ L (T) curves were fitted using the following equation: where λ L (0) is the London penetration depth at T = 0 K, T c = 14., k demag values for the λ L calculation, respectively.At low temperatures, λ L does not show the typical exponential behavior expected for a fully gapped clean s-wave superconductor [44].In general, a power law temperature dependence of λ L implies the presence of low-energy quasiparticle excitations [45].In this framework, when plotting  as a function of T (see Figure 4), the curve behavior was noted as being completely different from the behavior predicted using the single-gap BCS theory, which should have an opposite concavity; this indicates a more probable multigap superconductivity in our sample [25,27,46,47].This aligns with the results reported in the literature for Fe-based pnictides [48][49][50][51].On the other hand, it should be noted that the  (T) curvature may be also a sign that the sample has an anisotropic single gap nature [52].Finally, it is worth noting that the  values obtained from  appear low in respect to other Fe(Se,Te) samples reported in the literature [45,46,53,54], indicating a low vortex-vortex interaction that, in the framework of vortex dynamics, usually characterizes a single vortex state.This vortex lattice configuration, due to the strength of the effective pinning energy, typically allows the sample to carry high transport currents in high fields that are suitable for power applications of this class of superconducting materials.After the calculation of the  values, the coherence length ξ values were extracted, starting with the  2 (T) curve reported in Figure 12 of Ref. [4], using the following equation [55,56]: where  = 2.068 × 10 −15 Tm 2 is the magnetic flux quantum.In this framework, when plotting λ −2 L as a function of T (see Figure 4), the curve behavior was noted as being completely different from the behavior predicted using the single-gap BCS theory, which should have an opposite concavity; this indicates a more probable multigap superconductivity in our sample [25,27,46,47].This aligns with the results reported in the literature for Fe-based pnictides [48][49][50][51].On the other hand, it should be noted that the λ −2 L (T) curvature may be also a sign that the sample has an anisotropic single gap nature [52].Finally, it is worth noting that the λ demag L values obtained from H demag c1 appear low in respect to other Fe(Se,Te) samples reported in the literature [45,46,53,54], indicating a low vortex-vortex interaction that, in the framework of vortex dynamics, usually characterizes a single vortex state.This vortex lattice configuration, due to the strength of the effective pinning energy, typically allows the sample to carry high transport currents in high fields that are suitable for power applications of this class of superconducting materials.After the calculation of the λ L values, the coherence length ξ values were extracted, starting with the H c2 (T) curve reported in Figure 12 of Ref. [4], using the following equation [55,56]: where φ 0 = 2.068 × 10 −15 Tm 2 is the magnetic flux quantum.The ξ values, as a function of temperature, are shown in Figure 5 together with the fit with the equation [57]: where ξ(0) is the coherence length at 0 K, T c = 14.5 K [4], and n is the exponent.From the fitting procedure, we obtained ξ(0) ≈ 3 nm, which coheres with the value, ξ(0) ≈ 2.7 nm, which was calculated using H c2 (0) = 46.5 T from Ref. [4], and n ≈ −0.64, from Ref. [57].
The determination of the penetration depth, and the coherence length values, allows us to estimate the de-pairing current density, J dep .The de-pairing current density is the current value above which the superconductivity of the superconductor is completely broken; this is because it is the value wherein the kinetic energy of the superconducting carriers equals the binding energy of the Cooper pair.The de-pairing current density can be expressed, thanks to the Ginzburg-Landau theory, as follows [30,31]: where µ 0 is the vacuum permeability.It is worth underlining that this equation is usually only considered valid for temperatures approaching T c .Nevertheless, by taking into account the calculations performed by Bardeen [58], the de-pairing current density deviates from the GL theory by a maximum factor of 1.5 at T = 0 K; this indicates that it is possible to study its temperature behavior with reasonable error.Moreover, Equation ( 9) is usually valid for the single band superconductor, whereas Fe(Se,Te) shows peculiarities which can be ascribed to a multigap superconducting behavior, as previously mentioned.In this context, another important parameter to consider is γ H = H c2,ab /H c2,c , and in particular, its temperature dependence.In fact, in a single gap superconductor, the γ H parameter is temperature independent, while in a multiband superconductor, the contributions of different bands lead to a non-constant γ H , with each band contributing differently as T is changed.For Fe(Se,Te), it is evident in several works [59][60][61] that γ H is slightly dependent on temperature, therefore, making reasonable the use of the GL formula, although it gives overestimated J dep values, especially at low temperatures.The J dep values.obtained using Equation ( 9), are reported in Figure 6.It is worth noting the presence of two J dep (T) curves.In particular, the black curve was obtained using λ L without the demagnetization correction, whereas the red curve was obtained by considering λ demag L in conjunction with the demagnetization correction.Both sets of J dep values align with the highest values reported for the different iron-based families; this demonstrates the very good quality of this crystal [62][63][64][65][66].It is evident that the J demag dep values that take the demagnetizing factor into account are five times higher than the J dep values obtained without considering the demagnetizing factor.This helps us understand how important the demagnetizing factor is when estimating different important superconducting parameters such as H c1 , λ L , and J dep .It is important to note that the role of λ L is crucial to Equation (9) since it is squared, and therefore, a small λ L change generates a large J dep variation.In conclusion, in light of the fact that H c1 and λ L can strongly depend on stoichiometry, the possibility of tuning it by modifying the fabrication process and parameters could be exploited to enhance the de-pairing current density.The ξ values, as a function of temperature, are shown in Figure 5 together with the fit with the equation [57]: where ξ(0) is the coherence length at 0 K,  = 14.5 K [4], and n is the exponent.From the fitting procedure, we obtained ξ(0) ≈ 3 nm, which coheres with the value, ξ(0) ≈ 2.7 nm, which was calculated using  (0) = 46.5 T from Ref. [4], and  ≈ −0.64, from Ref. [57].values obtained without considering the demagnetizing factor.This helps us understand how important the demagnetizing factor is when estimating different important superconducting parameters such as  ,  , and  .It is important to note that the role of  is crucial to Equation ( 9) since it is squared, and therefore, a small  change generates a large  variation.In conclusion, in light of the fact that  and  can strongly depend on stoichiometry, the possibility of tuning it by modifying the fabrication process and parameters could be exploited to enhance the de-pairing current density.

Materials and Methods
We analyzed a FeSe0.5Te0.5 (nominal composition) crystal with the following dimensions: 3 × 3 × 0.2 mm 3 .The crystal was created using the Bridgman technique, and  = 14.5 K. Details concerning the creation of the crystal are reported elsewhere [4].A SEM-EDX analysis was performed on the sample, which showed the presence of twin boundaries and a slight deviation from the nominal composition in terms of stoichiometry (Fe0.96Te0.59Se0.45)[67].This is probably due to micro inhomogeneity and the phase values obtained without considering the demagnetizing factor.This helps us understand how important the demagnetizing factor is when estimating different important superconducting parameters such as  ,  , and  .It is important to note that the role of  is crucial to Equation ( 9) since it is squared, and therefore, a small  change generates a large  variation.In conclusion, in light of the fact that  and  can strongly depend on stoichiometry, the possibility of tuning it by modifying the fabrication process and parameters could be exploited to enhance the de-pairing current density.

Materials and Methods
We analyzed a FeSe0.5Te0.5 (nominal composition) crystal with the following dimensions: 3 × 3 × 0.2 mm 3 .The crystal was created using the Bridgman technique, and  = 14.5 K. Details concerning the creation of the crystal are reported elsewhere [4].A SEM-EDX analysis was performed on the sample, which showed the presence of twin boundaries and a slight deviation from the nominal composition in terms of stoichiometry (Fe0.96Te0.59Se0.45)[67].This is probably due to micro inhomogeneity and the phase

Materials and Methods
We analyzed a FeSe 0.5 Te 0.5 (nominal composition) crystal with the following dimensions: 3 × 3 × 0.2 mm 3 .The crystal was created using the Bridgman technique, and T c = 14.5 K. Details concerning the creation of the crystal are reported elsewhere [4].A SEM-EDX analysis was performed on the sample, which showed the presence of twin boundaries and a slight deviation from the nominal composition in terms of stoichiometry (Fe 0.96 Te 0.59 Se 0.45 ) [67].This is probably due to micro inhomogeneity and the phase separation of magnetic premises, which is typical for crystal growth and synthesis in Fe-SeTe [68][69][70] and its basic compound FeSe [71][72][73].The sample was characterized in a dc magnetic field that was applied perpendicularly to its largest face (H||c).In particular, the dc magnetic moment, as a function of the field, m(H), was measured using a Quantum Design PPMS-9T equipped with a VSM option.To avoid the effect on the sample response caused by the residual trapped field inside the PPMS dc magnet [74], this field was reduced below 1 × 10 −4 T [75].Regarding the m(H) measurements, the sample was first cooled down to the measurement temperature in the zero field and thermally stabilized for at least Condens.Matter 2023, 8, 91 8 of 11 20 min.Then, the field was ramped with the fixed sweep rate value to +9 T, then it was reduced to −9 T, and finally, to +9 T again in order to acquire the complete hysteresis loop.

Conclusions
By studying the DC magnetic moment as a function of the magnetic field (H) at different temperatures (T) on a Fe(Se,Te) crystal created using the Bridgman technique, the effect of the demagnetizing factor on the de-pairing current density J dep values was estimated.In order to achieve this, the London penetration depth and the coherence length were evaluated.Initially, the lower critical field, H c1 , was obtained as a function of temperature that focused on the impact that demagnetization effects have on its values.In particular, it was found that H c1 values were underestimated by a factor equal to four, even at T = 0 K. Starting with these results, the London penetration depth, λ L , was calculated as a function of temperature using both the original H c1 values and the values obtained after considering the demagnetization effects (H demag c1 ).The λ L (T) curves did not show the typical exponential behavior that was expected for a fully gapped clean s-wave superconductor, but rather, they exhibited a power law dependence that indicated the presence of low-energy quasiparticle excitations.In this framework, λ −2 L , as a function of T, was graphed, and multigap-like behavior was found in our sample, which aligns with the behavior of other iron-based samples reported in the literature.Additionally, we found that the λ demag L values obtained from H demag c1 are lower than the values obtained from other Fe(Se,Te) samples reported in the literature, which suggests a possible single vortex state.In this vortex lattice configuration, the effective pinning energy is high, therefore, reducing the dissipations inside the material, and improving the current transport properties is suitable for power applications of this class of superconducting materials.After that, the coherence length ξ values were extracted, starting with the H c2 (T) curve, as reported in our previous work.Combining the λ L and ξ values, the J dep values that took the demagnetization effects into account, as a function of temperature, were five times higher than the values that did not take the demagnetization effects into account.Tuning the stoichiometry of the compounds by modifying the fabrication procedures, and thus, changing the H c1 and λ L values by considering the demagnetizing factor, can push the limits of the de-pairing current density values of the materials.
0) is the value of H demag c1 at T = 0 K, T c = 14.5 K [4], and n is the exponent.From the fitting procedure, H demag c1 (0) = 597 Oe and n = 1.54 were obtained.The H demag c1 Condens.Matter 2023, 8, x FOR PEER REVIEW 3 of 12

Figure 1 .
Figure 1.The field dependence of the initial magnetic moment curve is plotted for different temperatures.The black dashed line provides the linear fit for the low field m(H) curves.

Figure 1 .
Figure 1.The field dependence of the initial magnetic moment curve is plotted for different temperatures.The black dashed line provides the linear fit for the low field m(H) curves.

Figure 2 .
Figure 2. (a) Comparison between the temperature dependence of the lower critical field  (black squares) (where the black solid line is the fit of the data using the following equation,  () =  (0) 1 − ) and the temperature dependence of the lower critical field after considering the demagnetization effects  (red circles) (where the red dashed line is the fit of the data using the following equation  () =  (0) 1 − ).When performing the fit,  (0) = 143 Oe and n = 1.54 were obtained without considering the demagnetizing factor, and  (0) = 597 Oe and n = 1.54 were obtained by considering the demagnetizing factor.(b) Ratio of the  (T) and  (T) values.The red star indicates the ratio of  (0)/ (0).

Figure 2 .
Figure 2. (a) Comparison between the temperature dependence of the lower critical field H c1 (black squares) (where the black solid line is the fit of the data using the following equation, H c1 (T) = H c1 (0) 1 − T T c

Figure 3 .
Figure 3. Temperature dependence of the London penetration depth λ L obtained from H c1 (black squares) and H demag c1 (red circles).Both the curves were fitted with λ L (T) = λ L (0) 1 − T T c

12 Figure 4 .
Figure 4. Temperature dependence of the London penetration depth  , obtained from  (black squares) and  (red circles).The solid lines are a guide for the eyes.

Figure 4 .
Figure 4. Temperature dependence of the London penetration depth λ −2 L , obtained from H c1 (black squares) and H demag c1(red circles).The solid lines are a guide for the eyes.

Figure 6 .
Figure 6.Temperature dependence of the de-pairing current density,  with (red closed circles) and without (black closed squares) consideration of the demagnetizing factor.The solid lines are a guide for the eyes.

Figure 5 .
Figure 5. Temperature dependence of the coherence length, ξ (black squares), together with its fit and the following equation ξ(T) = ξ(0) 1 − T T c

Figure 6 .
Figure 6.Temperature dependence of the de-pairing current density,  with (red closed circles) and without (black closed squares) consideration of the demagnetizing factor.The solid lines are a guide for the eyes.

Figure 6 .
Figure 6.Temperature dependence of the de-pairing current density, J dep with (red closed circles) and without (black closed squares) consideration of the demagnetizing factor.The solid lines are a guide for the eyes.