1 Introduction
The interplay of superconductivity and magnetism has been a research focus for decades. Usually, superconductivity and magnetic order will compete with each other, and thus superconductivity seldom occurs in strong magnetic materials [
1]. In conventional s-wave superconductors, local magnetic moments would break up spin singlet Cooper pairs and hence strongly suppress superconductivity [
2]. In the case of rare-earth superconductors, 4
f electrons may have a localized magnetic moment which could weakly couple with conduction electrons, and the coexistence of superconductivity and magnetism is possible to realize. There have been some examples of such materials, like
RMo
6S
8 [
3,
4],
RRh
4B
4 [
3,
4],
RBa
2Cu
3O
7−δ [
5–
8],
RNi
2B
2C [
9,
10],
RPdBi [
11] (
R = some of the rare earth elements), and iron-based superconductors containing rare earth elements [
12]. In ferromagnetic ErRh
4B
4 [
13] and HoMo
6S
8 [
14], and antiferromagntic (AFM) GdMo
6S
8 [
15], the competition between superconductivity and magnetic long-range order has resulted in the reentrant superconductivity behavior. In the
RNi
2B
2C and
RPdBi systems, the competitive relationship between AFM transition temperature
TN and superconducting (SC) transition temperature
Tc can be described by the Gennes scaling [
10,
11]. Moreover, the coexistence of magnetism and superconductivity may lead to some exotic phases. For example, UTe
2 where superconductivity coexists with ferromagnetic spin fluctuations is proposed to have a spin-triplet pairing state [
16], and the coexistence of magnetism and superconductivity in the
RPdBi systems can serve as a unique platform to investigate topological orders with multisymmetry breaking states [
11].
The
R2O
2Bi (
R = rare earth) families are a series of layered compounds with anti-ThCr
2Si
2 structure. Superconductivity has been found by excess oxygen incorporation in some of
R2O
2Bi members, irrespective of the presence of magnetic ordering [
17]. Superconductivity in non-magnetic Y
2O
2Bi was found dependent on the oxygen content [
18]. The oxygen incorporation will cause the expansion of inter-Bi-layer distance, which is proposed to be crucial to the occurrence of superconductivity. The obvious change in the carrier density companied with oxygen incorporation was not found. Meanwhile, the search for superconductivity in the magnetic
R2O
2Bi system has produced various outcomes. The Ce
2O
2Bi was found to be an AFM Kondo lattice compound without observation of superconductivity for temperature down to 0.3 K [
19]. Although Er
2O
2Bi, Tb
2O
2Bi and Dy
2O
2Bi have been reported to be superconducting, the detailed investigation of their properties is lacked [
17,
20]. Whether the 4
f electron magnetism has remarkable influence on superconductivity is an interesting issue and it is worthy of careful study. Up to now, among the magnetic
R2O
2Bi members, only the Er
2O
2Bi is found to be superconducting without the CaO incorporation which serves as an oxidant to obtain excess oxygen in
R2O
2Bi [
17]. However, the Ca ion impurity would also be incorporated in the samples and it may affect the superconducting properties. Thus it is better to tune the oxygen content without CaO incorporation in studying the interplay of magnetism and superconductivity in a series of
R2O
2Bi.
Here, we systematically studied the Er2O2Bi polycrystalline compounds with various oxygen content by performing magnetization, specific heat and resistivity measurements. The AFM transition associated with the layer occurs at TN of about 3 K, which indicates that there is a coexistence of superconductivity and AFM order in Er2O2Bi. The first-principle calculations suggest that the Fermi surface of Er2O2Bi is similar to that of Y2O2Bi, implying superconductivity in both compounds could share the same mechanism and the conducting monatomic Bi2−square net layer may be crucial to the occurrence of superconductivity. Above TN, there seems a kind of short-range magnetic correlation which enhances the specific heat significantly and reduces the magnetic entropy above TN. Furthermore, we have synthesized four samples with different nominal oxygen content x in the Er2OxBi compounds (x = 1.8, 1.9, 2.0, and 2.1). Upon increasing x, the actual oxygen content indeed increases, and TN is slightly suppressed. Meanwhile the hole-type carrier density and Tc both increase, indicating that the AFM order may compete with superconductivity. The variation of hole-type carrier density may account for the change of Tc and TN. This work provides a new example to study the interplay of rare earth magnetism and superconductivity in such an anti-ThCr2Si2-type structure.
2 Experimental details
The Er2OxBi (x = 1.8, 1.9, 2.0, 2.1) polycrystalline compounds were synthesized by solid-state reactions. High-purity Er2O3 (Alfa, 99.99%), Er (Alfa 99.8%) and Bi (Alfa 99.999%) were used as starting materials. The mixed powders were pressed into pellets and covered with Ta foils. An excess amount of Bi was added because of its high vapor pressure. The pellets was sealed in an evaluated silica tube, and heated at 1000◦C for 20 h in furnace. After it had cooled down to room temperature naturally, the pellets was reground and pressed again, then annealed at 1000◦C to improve the sample homogeneity. All the operations were taken in an Ar-filled glove box except the heating processes.
The powder X-ray diffraction (XRD) patterns were recorded on a PANalytical X-ray diffractometer with Cu K
α radiation at room temperature. We used the Rietveld refinements employing program GSAS with EXPGUI [
21] to obtain lattice parameters. The DC magnetization was measured on two magnetic property measurement systems (MPMS-XL5 for
T > 1.8 K and MPMS3 with
3He refrigerator for
T < 1.8 K ). Electrical resistivity and Hall effect were measured by a standard four-probe method on a physical property measurement system (PPMS-9). Specific heat measurements were also carried out on the same PPMS-9 by a heat-pulse relaxation method. We also used an Oxford superconducting magnet system equipped with a
3He cryostat to measure the superconducting transitions.
The first-principle calculations were performed with the Vienna Ab initio Simulation Package (VASP) code [
22,
23]. Throughout the calculations, the Perdew, Burke, and Ernzerhof (PBE) parameterization of generalized gradient approximation (GGA) to the exchange correlation function was used [
24]. The spin–orbit coupling (SOC) is included using the second variational method. The energy cutoff of plane-wave basis was chosen to be 600 eV, which was sufficient to converge the total energy to 1 meV/atom. A Γ-centered 18 × 18 × 18 Monkhorst–Pack [
25] k-point mesh was chosen to sample the Brillouin zone in the calculations. Considering that the Er-4
f states are fully localized at high temperature, the paramagnetic (PM) band structure was calculated by assuming the Er-4
f open-core electron. While for the AFM state, the states of
f electrons were treated as semi-core valence states by performing LDA+
U with
U = 6 eV [
26–
29].
3 Results and discussion
Figure 1(a) shows the XRD patterns of Er
2O
xBi with different nominal oxygen content (
x). It is obvious that all the four samples have the same main phase of anti-ThCr
2Si
2 type structure with a few minor impurity phases depending on the oxygen content. From the XRD patterns, the sample with
x = 2.0 has the best quality with least amount of impurity phases, therefore the most of superconducting property’ characterization is performed based on this sample. We used the Rietveld refinements employing program GSAS with EXPGUI [
21] to obtain lattice parameters of the four samples. Figure 1(b) shows the Rietveld refinement of the XRD pattern of the sample Er
2O
2Bi. All diffraction peaks can be indexed with the crystalline structure of Er
2O
2Bi, and their (
h k l) indices are also indicated for each peak. The fitting parameters are
a = 3.8459(1) Å and
c = 13.1585(5) Å with the parameters
Rwp = 7.65% and
χ2 = 1.93 which indicate a good convergence. The inset of Fig. 1(a) displays the variations of the lattice constants
a and
c with nominal oxygen content. As the oxygen content
x increases, there is a remarkable increase in the lattice constant
c, and the lattice constant
a remains almost constant. This
c-axis expansion due to excess oxygen incorporation is also observed in Y
2O
2Bi [
18]. The inset of Fig. 1(b) displays the anti-ThCr
2Si
2-type structure of Er
2O
2Bi which is composed of a Bi square net layer and two ErO layers.
Figure 2(a) shows the temperature-dependent resistivity for the sample Er
2O
2Bi from 300 K to 0.4 K. From the inset of Fig. 2(a), there are two obvious transitions in resistivity which should correspond to
Tc and
TN, respectively. The
χ–
T curves around
TN are also plotted in the inset for comparison, which indicates the good coincidence between the transitions in both resistivity and magnetic susceptibility at
TN. The increased magnetic field does not depress the transition in magnetic susceptibility around
T of 3 K, excluding a possible superconducting transition. Figure 2(b) shows the resistivity of Er
2O
2 Bi near the
Tc under different magnetic fields. The resistivity decreases to zero abruptly below the SC transition temperature
Tc = 1.23 K under zero magnetic filed. Here we define
Tc as the temperature at which the resistivity drops to 50% of normal-state resistivity just above the SC transition. The
(zero resistivity criterion) is 1.17 K. With the magnetic field increasing,
Tc is suppressed continuously. The reduced temperature (
T /
Tc) dependence of upper critical magnetic field
Hc2 is presented in Fig. 2(c). We tried to use the one-band Ginzburg–Landau theory,
Hc2 =
Hc2(0)(1−
t2)/(1 +
t2) (where
t =
T /
Tc), to fit the upper critical field and the fitting result is not good. It is found that there seems to be a positive curvature close to
Tc, which could be a characteristic of multi-band superconductor [
30–
33]. Therefore, we introduce the two-band Ginzburg–Landau theory [
34],
Here a0 = 2(λ11λ22–λ12λ21), a1 = 1+(λ11–λ22)/λ0, a2 = 1–(λ11–λ22)/λ0, , h =Hc2D1/2ϕ0T , η = D1/D2, U (x) = Ψ(x+1/2)–Ψ(x), Ψ(x) is the digamma function. λ11 and λ22 are the intraband BCS coupling constants, λ12 and λ21 are the interband BCS coupling constants, and D1 and D2 are the in-plane diffusivity of each band. The fitting parameters are listed in Table 1.
Based on this two-band model, we can obtain
Hc2(0) = 0.061 T which is lower than the BCS Pauli limit, i.e.,
= 1.84
Tc∼2.3 T. Compared to the upper critical field of Y
2O
2Bi which is about 0.0574 T [
18],
Hc2(0) of Er
2O
2Bi does not show much difference, in spite of the substitution of non-magnetic Y
3+ ion by the magnetic Er
3+ ion. It indicates that superconductivity in both compounds could share the same mechanism and the conducting monatomic Bi
2−square net layer may be crucial to the occurrence of superconductivity. Then, we estimate Ginzburg–Landau coherence length by
Hc2 = Φ
0/(2
πξ2) and
ξ(0) of about 730 Å is obtained. Apparently
ξ(0) is much larger than the lattice constants, while the length of magnetic unit cell is usually comparable to the lattice constants. The larger
ξ(0) value (than the length of magnetic unit cell) could account for the coexistence of superconductivity and antiferromagnetism [
10].
Figure 2(d) shows the temperature dependence of magnetic susceptibility χ (left axis) and inverse magnetic susceptibility, 1/χ (right axis) measured under the magnetic field H = 100 Oe. The χ(T ) satisfies the modified Curie–Weiss law χ = χ0 + C/(T−θW) where χ0 is a temperature independent susceptibility, C is the Curie constant and θW is the Weiss temperature. The red solid line shows the Curie–Weiss fit in the temperature range between 7–300 K and the derived parameters C = 11.77 emu/mol-Er K, χ0 = 0.00372 emu/mol-Er and θW = – 22 K. The negative Weiss temperature indicates an AFM exchange interaction between the magnetic moments. The calculated effective magnetic moment based on the Curie constant is 9.7 µB /Er which is very close to that of the free Er3+ ion (9.8 µB /Er). It implies that the 4f electrons of Er3+ ion are very localized in the temperature range between 7–300 K. Figure 2(e) displays the temperature dependence of magnetic susceptibility χ from T = 0.4 K to 5 K. It is clear that there exist two transitions which correspond to the AFM order and SC transition respectively. The red hollow circles and black solid circles stand for the field cooling (FC) and zero-field cooling (ZFC) processes under the magnetic field H = 10 Oe around Tc, respectively. The data represented by the blue solid squares is measured under H = 100 Oe to show the change around TN. To avoid demagnetization effect, the sample is shaped into a long and narrow cuboid (about 3 mm× 0.66 mm × 0.4 mm) and the magnetic field is applied along the long axis of the sample. The SC volume fraction is estimated to be about 133% based on the ZFC data, confirming the bulk nature of superconductivity in Er2O2Bi. Although the demagnetization factor is close to 0, the overestimated SC volume fraction which is beyond 100% may be caused by the error in measuring the sample volume. The inset of Fig. 2(e) is the field dependence of magnetization M (H) at T = 0.5 K which is a typical behavior of type-II superconductor. Figure 2(f) shows the magnetic susceptibility χ measured in the temperature T -sweep mode containing both T -up and T -down processes in the temperature range between 2–4 K. To avoid the error of T , the speed of T sweep is as slow as 0.15 K/min. There is a clear magnetic hysteresis around T of 3 K under 100 Oe, suggesting a possible first-order transition.
Figure 3(a) shows the temperature dependence of specific heat
C of Er
2O
2Bi. The peak at
TN is so sharp, indicating a possible first-order transition. Such a sharp peak in specific heat has been observed in the magnetic first order phase transitions in other compounds, i.e., DyCu
2Si
2 [
35]. With the magnetic field increasing, the peak is suppressed gradually, consistent with the AFM order. As shown in the inset of Fig. 3(a),
TN decreases monotonously with increasing magnetic field. Figure 3(b) displays the temperature dependence of specific heat
C of Er
2O
2Bi among the temperature range 0–2 K under magnetic fields
µ0H = 0, 0.01, 0.05 and 0.1 T, respectively. There is no obvious peak of specific heat observed at the SC transition temperature
Tc = 1.23 K, which is probably hidden in the shoulder of the large AFM peak. However, from the inset of Fig. 3(b), the different behaviors of
C/
T at the lowest temperature may be induced by SC transition. The
C/
T under zero magnetic field has a drop around 0.6 K, while the drop of
C/
T disappeared under
µ0H = 0.01 T, and
C/
T become coincident under
µ0H = 0.05 and 0.1 T. Concerning the
C/
T behavior in Y
2O
2Bi [
18], there is also an upturn below about 0.6 K which is proposed to be caused by a Schottky anomaly originating from the
209Bi nuclei (
I = 9/2). Therefore, we may ascribe the upturn below 0.6 K in
C/
T of Er
2O
2 Bi to the Bi nuclear contributions and the drop below 0.6 K under 0 T to the SC transition. With increasing magnetic field, the SC transition is suppressed and the upturn behavior gradually appears. When the magnetic field is beyond
Hc2(
T ), i.e.,
µ0H = 0.05 and 0.1 T, the drop in
C/
T induced by the SC transition disappears and the
C/
T would remain almost the same with further increasing magnetic field.
Figure 3(c) shows the specific heat divided by temperature, C/T , versus T 2 under the magnetic field µ0H = 0 and 9 T. The following equation is used to fit the C/T data,
where
γ and
β denote the coefficient of electron and phonon contributions for Er
2O
2Bi, respectively. The derived parameters are
γ = 324 (425) mJ/(mol·K
2) and
β = 3.1 (3.7) mJ/(mol·K
4) under the magnetic field
µ0H = 0 (9) T. The derived Debye temperatures Θ
D is 146 K and 138 K under
µ0H = 0 and 9 T respectively, by using the formula Θ
D = [12
π4N R/(5
β)]
1/3. Since the magnetic field usually should not change the phonon contribution to specific heat, the very close values of Debye temperatures Θ
D for
µ0H = 0 T and 9 T are reasonable. The Debye temperature of Er
2O
2Bi (146 K) is a little smaller than that of La
2O
2Bi, which is about 150 K [
19]. The difference between their Debye temperatures may result from a little larger mass of Er ions compared to the La ions. Based on Θ
D = 146 K, the electron–phonon coupling constant
λe−ph is estimated to be 0.53 from the McMillan equation [
36]:
where µ* is the Coulomb pseudopotential, and an empirical value of 0.13 is used. This value of λe−ph in Er2O2Bi is as the same as the value λe−ph = 0.53 in Y2O2Bi, both pointing to the typical weak coupling conventional BCS superconductors.
The Sommerfeld coefficient
γ∼324 mJ/(mol·K
2) (under zero magnetic field) is quite large. Usually, the large
γ value may result from Kondo interactions or strong magnetic fluctuations, which is difficult to distinguish them only based on the specific heat data. After carefully investigating the behaviors of resistivity and magnetic susceptibility, we argue that it should be caused by the magnetic fluctuations rather than the Kondo interactions. There is no obvious feature of Kondo effect in the temperature dependence of resistivity, such as a negative logarithmic temperature dependence of resistivity induced by Kondo scattering. The magnetic moments of Er
3+ are localized from 300 K to 7 K derived from magnetic susceptibility, which is not consistent with the Kondo effect model (magnetic moments will be partially screened by the conduction electrons at low temperature). Instead, it seems that the magnetic fluctuations could account for the large
γ, which is supported by the calculated magnetic entropy (to be shown below). Considering the contribution to specific heat from spin fluctuations in Er
2O
2Bi do not include a
T 3 ln
T term which is a character of ferromagnetic fluctuations, the spin fluctuations here may be antiferromagnetic type [
37].
Figure 3(d) shows the magnetic entropy calculated by integrating Ce/T over temperature. The Ce denotes the electronic contribution to specific heat, which is separated from the phonon contribution. The phonon contributions are calculated by the Debye model with ΘD = 146 K. The magnetic entropy does not recover R ln 2 at TN and keeps increasing, and finally saturates as temperature is close to 18 K. Assuming that the ground state of Er2O2Bi is Kramers two-fold degenerate, the free magnetic entropy should be R ln 2. The reduced magnetic entropy at TN may be interpreted by a kind of short-range magnetic correlation above the AFM long-range order temperature (TN). The saturation value of magnetic entropy at T = 18 K is a little smaller than R ln 2, which is probably introduced by the error in calculating the phonon contributions based on the Debye model.
Figure 4(a) shows the temperature dependence of Hall resistivity for the sample Er2O2Bi. The magnetic field dependence of Hall resistivity at high temperature is nearly linear which can be fitted by the usual one-band model. However, when the temperature decreases, the magnetic field dependence becomes nonlinear, which could be induced by magnetic skew-scattering and the effects of two types of charge carriers. Figure 4(b) displays the temperature dependence of Hall coefficient of Er2OxBi with x= 1.8, 1.9 and 2.0, respectively. With increasing nominal oxygen content (x), the positive Hall coefficient is decreasing, implying the increase in the hole-type carrier density. Since the samples with x = 2.1 and x = 2.0 have almost the same parameters of c-axis, and meanwhile their Tc and TN are also the same, we believe that both samples should have the same oxygen content, and the excess oxygen content in the sample with nominal x = 2.1 may lead to the formation of a few impurity phases like Er2O3. Thus the sample with x = 2.1 does not have good quality and its Hall coefficient is not represented for discussions. Figure 4(c) shows the x dependence of carrier density nH and mobility µ calculated based the Hall coefficient and resistivity at 100 K. The carrier density is calculated by a single band approximation. It is obvious that there is an increase in the carrier density as x increases from 1.8 to 2.0, which means the oxygen incorporation indeed induce more hole-type charge carriers.
Figure 5(a) shows temperature dependence of resistivity of the samples Er2OxBi with different oxygen content x. As x increases, Tc increases from 1 K to 1.23 K. We use the peak of T dχ/dT displayed in Fig. 5(b) to define the TN. In contrast to the SC transition, TN decreases with increasing x. The parameters, nH, Tc and TN of the samples are presented in the Fig. 5(c). The increase in the carrier density can account for the increase of Tc and the decrease of TN, which reveals that the AFM order may compete with superconductivity.
The primitive unit cell (with the simplest AFM configuration) and Brillouin zone of Er2O2Bi are shown in Figs. 6(a) and (b). The band structures of PM phase (black line) and AFM phase (red line) are shown in Fig. 6(c). In the AFM phase, the Er-4f orbitals are away from the Fermi level, located around 6 eV (1 eV) below (above) Fermi energy EF. Thus the electronic structure of AFM phase resembles that of PM phase near Fermi surface with two bands derived from Bi-p orbits crossing EF, which indicates that the magnetic moments of Er3+ ions probably have little effect on the electronic structure near the Fermi level. This result may account for the similar electron-phonon coupling of Er2O2Bi and Y2O2Bi. By integrating the 4f contribution below EF, each Er ion has a spin polarization of 2.9 µB in the z direction, corresponding to a spin configuration of S = 3/2. Therefore, the effective magnetic moment can be obtained with µeff = . In this calculation, it is found gj = 1.2, resulting in µef f = 9.58 µB , in good agreement with above experimental data, implying the fully localized nature of 4f electrons in Er2O2Bi in the AFM ground state.
4 Conclusion
In summary, we have systematically investigated the superconducting properties and AFM order of the Er2O2Bi compounds with different oxygen content by performing resistivity, magnetization, specific heat measurements and first-principle calculations. There may exist strong magnetic fluctuations or short-range AFM correlation above TN and they contributes to the significantly enhanced Sommerfeld coefficient of specific heat. The magnitude of Hc2(0) and electron-phonon coupling constants in both Er2O2Bi and Y2O2Bi are very close, implying that superconductivity in both magnetic or nonmagnetic R2O2Bi (R= Y and Er) may share the same origin. The band structures from the first-principle calculations support the view that the charge carriers in Er2O2Bi mainly come from the Bi square net layer. Therefore, the conducting monatomic Bi2−square net layer may be crucial to the occurrence of superconductivity in R2O2Bi. By adjusting the nominal oxygen content in Er2OxBi, we obtained several samples with different carrier density. The Er2O2Bi compound is a typical example with the coexistence of superconductivity and AFM order in the series of R2O2Bi (R = rare earth or Y) materials. This study suggested that there may be a competition between superconductivity and AFM order in Er2O2Bi. It is worth further studying to understand the interplay of rare earth magnetism and superconductivity in such an anti-ThCr2Si2-type structure.