Next Article in Journal
The Effect of 10% Carbamide Peroxide Dental Bleaching on the Physical Properties of Invisalign Aligners: An In Vitro Study
Previous Article in Journal
Electrical Resistance Performance of Cable Accessory Interface Considering Thermal Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Role of Lanthanides and Bilayer Fe2As2 in the Electronic Properties of RbLn2Fe4As4O2 (Ln = Gd, Tb, and Dy) Superconductors

1
Department of Physics, School of Science, Zhejiang University of Science and Technology, Hangzhou 310023, China
2
Department of Physics, School of Sciences, Nantong University, Nantong 226019, China
3
Henan Provincial Key Laboratory of Smart Lighting, School of Information Engineering, Huanghuai University, Zhumadian 463000, China
4
Qilu Institute of Technology, School of Chemical and Biological Engineering, Jinan 250200, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(11), 4123; https://doi.org/10.3390/ma16114123
Submission received: 8 May 2023 / Revised: 28 May 2023 / Accepted: 29 May 2023 / Published: 1 June 2023

Abstract

:
The superconducting transition temperatures (T c ) of RbGd 2 Fe 4 As 4 O 2 , RbTb 2 Fe 4 As 4 O 2 , and RbDy 2 Fe 4 As 4 O 2 are 35 K, 34.7 K, and 34.3 K without doping, respectively. For the first time, we have studied the high-temperature nonmagnetic state and the low-temperature magnetic ground state of 12442 materials, RbTb 2 Fe 4 As 4 O 2 and RbDy 2 Fe 4 As 4 O 2 , using first principles calculations and comparing them with RbGd 2 Fe 4 As 4 O 2 . We also performed a detailed study of the effects of lanthanides and bilayer Fe 2 As 2 . We predict that the ground state of RbLn 2 Fe 4 As 4 O 2 (Ln = Gd, Tb, and Dy) is spin-density-wave-type, in-plane, striped antiferromagnets, and the magnetic moment around each Fe atom is about 2 μ B . We also found that the structural differences caused by the simple ionic radius have little effect on the properties of these three materials. Different lanthanide elements themselves play a major role in the electronic properties of the materials. It can be confirmed that the effect of Gd on RbLn 2 Fe 4 As 4 O 2 is indeed different from that of Tb and Dy, and the presence of Gd is more conducive to interlayer electron transfer. This means that Gd can transfer more electrons from the GdO layer to the FeAs layer compared to Tb and Dy. Therefore, RbGd 2 Fe 4 As 4 O 2 has a stronger internal coupling strength of the bilayer Fe 2 As 2 layer. This can explain why the T c of RbGd 2 Fe 4 As 4 O 2 is slightly higher than that of RbTb 2 Fe 4 As 4 O 2 and RbDy 2 Fe 4 As 4 O 2 .

1. Introduction

Finding new superconductors is an important and challenging task. After long-term efforts, especially in the past three decades, two basic routes for finding new superconducting materials have been discovered. The first is the chemical doping of potentially superconducting parent compounds, and the other is to synthesize new compounds containing two-dimensional superconducting active layers through structural design. Well-known examples are copper-based and iron-based high-temperature superconductors whose conductive active layers are CuO 2 planes and Fe 2 X 2 layers, respectively. Impressively, these far-reaching discoveries were made by chemical doping, and the structures of these high-temperature superconductors can be explained well by their crystal chemistry. Optimized doping levels and building block layers may set new records related to the superconducting transition temperature T c [1]. However, it is worth noting that superconductivity can also occur in undoped materials, where charge carriers are induced by intrinsic charge transfer through self-doping.
In recent years, the discovery of iron-based superconductors (FeSCs) has sparked a great deal of research activity. So far, dozens of FeSC crystal types have been discovered. Representative systems (with numerical abbreviations corresponding to their chemical compositions) that have been investigate in depth include (i) FeSe (11) [2], (ii) LiFeAs (111) [3], (iii) BaFe 2 As 2 (122) [4], (iv) K x Fe 2 y Se 2 (122*) [5], and (v) LaFeAsO (1111) [6,7]. All of these compounds consist of fluorite-like layers of Fe 2 X 2 (X = phosphorus group or chalcogenide); thus, the Fe 2 X 2 layer is considered a key building block for Fe-based superconductivity. This is further supported by the fact that some other compounds containing the same Fe 2 X 2 layer (but with a different barrier layer in between) also show superconductivity or become superconductors after proper chemical doping, with a relatively high T c . The 1111, 122, 111, and 11 iron-based high-temperature superconductors that have been discovered so far are all quasi-two-dimensional layered structures with Fe 2 X 2 layers. It was found that by adding calcium-state mineral-type layers between Fe 2 As 2 layers, the so-called 32522 and 42622 new systems can be obtained, and their typical representatives are Sr 3 Sc 2 O 5 Fe 2 As 2 [8] and Sr 4 V 2 O 6 Fe 2 As 2 [9]. It is also possible to insert a less symmetrical spacer layer between the Fe 2 As 2 layers.
So far, Cao Guanghan’s group [10] has designed several new structures containing Fe 2 X 2 layers. One of the candidate structures, “KLaFe 4 X 4 ”, was realized in the form of AeAFe 4 As 4 (Ae = K, Rb, and Cs; A = Ca and Sr) [11] and AkEuFe 4 As 4 (Ak=Rb and Cs) [12]. Another proposed structure, Ln 3 Fe 4 X 4 Z 2 [13], was realized in KCa 2 Fe 4 As 4 F 2 (hereafter referred to as 12442). At the same time, they also synthesized LnCa 2 Fe 4 As 4 F 2 (where Ln = Rb and Cs) [14], which are two sister compounds of KCa 2 Fe 4 As 4 F 2 , and they can be regarded as the interbreeding body of LnFe 2 As 2 (Ln = Rb and Cs) and CaFeAsF. This crystal structure has bilayer Fe 2 As 2 layers separated by insulating fluorite-type Ca 2 F 2 plates, similar to the case of double CuO 2 planes in cuprate superconductors. Its crystallization parameters were measured. Physical property measurements show that the bulk RbGd 2 Fe 4 As 4 O 2 , RbTb 2 Fe 4 As 4 O 2 , and RbDy 2 Fe 4 As 4 O 2 have T c values of 35 K, 34.7 K, and 34.3 K, respectively [15]. The development of novel superconductors whose bulk superconducting structures match the Fe 2 X 2 layers and satisfy thermodynamic or at least kinetic stability is important and challenging.
The 12442 material is synthesized from 1111 and 122 materials, AkLn 2 Fe 4 As 4 O 2 (Ak = K and Cs; Ln = Lanthanides) [15] iron-based superconducting material. This material exhibits superconductivity without extrinsic doping and without spin density wave (SDW) anomalies, which may imply a different superconducting mechanism than previous FeAs-based superconductors. Its crystal structure is unique in that it has bilayer Fe 2 As 2 layers [16]. At present, it has been found that a series of unique behaviors of 12442 may be related to these bilayer Fe 2 As 2 layers, such as the T c not being greatly affected by the purity [17] of the material; it not being greatly affected by the magnetic properties of lanthanide atoms; and the relationship of the T c to the bond angle of Fe-As and the height of the As atom relative to the Fe layer is inconsistent with the empirical rules of the previous numerical study [18]. At present, Ghosh et al. have discussed the electronic properties of KCa 2 Fe 4 As 4 F 2 [19], RbCa 2 Fe 4 As 4 F 2 , and CsCa 2 Fe 4 As 4 F 2 , and studies have shown that 12442 compounds have mixed multi-orbital and multi-band properties, and the contributions of these compounds are mainly from Fe-3D orbitals [20]. Unlike other hybrid iron-based superconductors, the contribution of As-4p z orbitals is negligible in these compounds. The substitution of base atoms with progressively larger atomic radii exerts chemical pressure inside the compound and induces orbital-selective evolution of the band structure and density of states [20]. However, there is no report on the electronic properties of RbTb 2 Fe 4 As 4 O 2 and RbDy 2 Fe 4 As 4 O 2 . At present, there are few articles on the electronic properties of RbGd 2 Fe 4 As 4 O 2 , and only the case of the nonmagnetic material without U is considered; the calculation results under magnetic conditions have not been considered [21].
In this paper, first principles methods, described in Section II, are applied to study the electronic structures and magnetic ground states of RbLn 2 Fe 4 As 4 O 2 . Section III contains the main results. It is confirmed that the ground states of RbLn 2 Fe 4 As 4 O 2 are SDW-type, in-plane, striped antiferromagnets, and the magnetic moment around each Fe atom is about 2 μ B . Several other differences between the Gd, Tb, and Dy compounds were analyzed, including the effect of the chemical pressure brought about by their introduction on the FeAs conductive layer and the role of differential charge transfer in the structurally unique double FeAs layer of 12442. A brief summary is provided in Section 4.

2. Materials and Methods

2.1. Crystal Structure

For all the calculations, the following experimental lattice parameters [22] are used in the unit cell: a = 3.901 Å and c = 31.343 Å for RbGd 2 Fe 4 As 4 O 2 , a = 3.890 Å and c = 31.277 Å for RbTb 2 Fe 4 As 4 O 2 , and a = 3.879 Å and c = 31.265 Å for RbDy 2 Fe 4 As 4 O 2 in nonmagnetic (NM) phase, ferromagnetic (FM) phases, and antiferromagnetic (AFM) phases. The shorter the Fe-As bond length, the stronger the metallicity, the stronger the ability to lose electrons, and the weaker the ability of the metal to control the outermost electrons. Therefore, the interaction between cations and free electrons becomes weaker, and the electron correlation decreases. The space group of NM, FM1, and FM2 are I4/mmm (No. 139), while the space group of AFM1, AFM2, AFM3 (AFM4), and AFM5 are Cmma (No. 67), P42/mmc (No. 131), I4/mm (No. 107), and P4/mmm (No. 123), respectively. Among them, AFM1 represents stripe antiferromagnetism (SAFM), that is, each Fe atom is aligned antiparallel to its second nearest neighbors along the a and b axes. AFM2 represents N e ´ el antiferromagnetism (NAFM); that is, the spin directions of the nearest neighbor atoms in the Fe layer are opposite. AFM3, AFM4, and AFM5 are different interlayer antiferromagnets. The space group belongs to P4/nmm, where FeAs layers and LnO layers are alternately arranged along the c-axis, forming a quasi-binary layered structure. The atoms in each layer are arranged in a four-coordinated tetrahedral structure and connected collinearly; in the FeAs layer, the 3D transition ferromagnetic atoms Fe form a two-dimensional square lattice, and charges are transferred from the LnO layer to the FeAs layer. The FeAs layer balances the chemical bonds between the (LnO) 1 + layer and the (FeAs) 1 layer.
The structures of RbLn 2 Fe 4 As 4 O 2 (Ln = Gd, Tb, and Dy) in the NM phase are shown in Figure 1. The As atoms in the FeAs bulk are no longer crystallographically equivalent. Here, the atomic coordinates are as follows: Rb (0, 0, 0); Ln (0.5, 0.5, z); Fe (0.5, 0, z); As1 (0.5, 0.5, z); As2 (0, 0, z); O (0.5, 0, 0.25). The antiferromagnetism of two Ln atoms is obtained by breaking the symmetry of the NM structure. As shown in Figure 1, we divide the FeAs layers in the electronic structure into an Fe 2 As 2 intralayer and an Fe 2 As 2 interlayer. It can be seen from the experimental results [22] that from Gd, Tb, to right Dy, both the Fe 2 As 2 interlayer (Fe1As1) distance and the Fe 2 As 2 intralayer (Fe2As2) distance become smaller, and these changes affect the superconducting transition temperature of the material.

2.2. Calculational Methods

The full-potential linearized augmented plane wave Wien2K package [23] was used for the electronic structure calculations. We used the Perdew, Burke, and Ernzerhof (PBE) [24] version of the generalized gradient approximation (GGA) within density functional theory. The sphere radii for Rb, Gd, Tb, Dy, and O were taken as 2.5 bohr, 2.33 bohr, 2.32 bohr, 2.32 bohr, and 1.9 bohr and as 2.3 and 2.19 bohr for Fe and As, respectively. The basis set cut-off parameter R m t · K m a x = 7.0 was found to be sufficient, and the number of k points was 8000 for the tetragonal unit cell, 3000 for SAFM supercell, and 4500 for other supercells that were used for each of the calculations to provide an adequate sampling of small Fermi surfaces. The supercell shapes and spin distributions studied under all magnetisms are shown in Figure 2. Considering the strong correlation of the 4f electrons in Gd, Tb, and Dy, we applied the GGA+U (U = 9 eV) code for RbGd 2 Fe 4 As 4 O 2 , RbTb 2 Fe 4 As 4 O 2 , and RbDy 2 Fe 4 As 4 O 2 . We studied several references [25,26,27,28] and also calculated U = 6 eV and U = 12 eV at the same time. After comparison, it was found that U = 9 eV is more reasonable for these three materials.

3. Analysis and Discussion

We first present for the total energy difference of NM, FM, and several different SAFM configurations of RbTb 2 Fe 4 As 4 O 2 and RbDy 2 Fe 4 As 4 O 2 , as well as for Fe in each magnetic phase. Table 1 shows the magnetic moments. We found that the SAFM configuration within the layer has the lowest total energy, so all of the ground states of RbLn 2 Fe 4 As 4 O 2 are in the in-plane SAFM phase. The magnetic moment around each Fe atom is about 2 μ B , which is much smaller than the local density approximation (LDA) value in LaFeAsO [29]. As we will see later, DFT may be suitable for this compound due to the presence of bilayer Fe 2 As 2 layers, RbLn 2 Fe 4 As 4 O 2 , which have better metallic properties and smaller magnetic moments than LaOFeAs. In the following, we will discuss the basic characteristics of the electronic structures of NM phases and the ground state in-plane SAFM in RbGd 2 Fe 4 As 4 O 2 , RbTb 2 Fe 4 As 4 O 2 , and RbDy 2 Fe 4 As 4 O 2 and compare the differences between these three materials. We also analyze the effects of lanthanides and the role of bilayer Fe 2 As 2 in detail.

3.1. Nonmagnetic Phase

Figure 3 shows the project of bands (pbands) of the NM phases RbLn 2 Fe 4 As 4 O 2 . We mainly discuss the energy bands contributed by Fe, and the band structure consists of four hole-like energy bands around the Γ (0, 0) point and four electron-like energy bands around the X ( π , π ) point. The NM band structures of RbLn 2 Fe 4 As 4 O 2 are similar to those of 1111-type compounds such as LaFeAsO, which means that their ground states also tend toward an antiferromagnetic spin density wave state. The band dispersion of electrons along the Γ -Z direction of RbTb 2 Fe 4 As 4 O 2 and RbDy 2 Fe 4 As 4 O 2 is flat, showing considerable two-dimensional features, while the hole band structure of RbGd 2 Fe 4 As 4 O 2 along the Γ -Z direction shows dispersion. It can be seen that only RbGd 2 Fe 4 As 4 O 2 has an energy band contributed by Fe 3d z 2 near the Fermi level at the point closest to Γ . It shows that compared with Tb and Dy, Gd is more conducive to the occupation of electrons in the 3d z 2 orbital, which is beneficial to the interlayer charge transfer. This energy band also corresponds to the Fermi surface at the center of Figure 3d. Significantly different from the Fermi surfaces of conventional iron-based superconductors, such as the 1111 and 122 systems, there is an additional pocket of holes near the region center Γ (0, 0, 0). Compared with the three materials, only Gd has this Fermi surface. We further analyze the pDOS of Fe corresponding to five orbitals in the three materials, and we can find that in RbTb 2 Fe 4 As 4 O 2 , the three orbitals of 3d x 2 y 2 , 3d x z , and 3d y z are the main contributions at the Fermi energy and near the Fermi energy: 3d x 2 y 2 > 3d x z > 3d y z . Similar to RbDy 2 Fe 4 As 4 O 2 , we can see that the main source of contribution near the Fermi level is also the three orbitals of Fe 3d x 2 y 2 , 3d x z , and 3d y z , and the contributions of the three orbitals are in the same order, that is, 3d x 2 y 2 > 3d x z > 3d y z . However, the difference in the DOS of each orbital is larger than that of RbTb 2 Fe 4 As 4 O 2 . From the pDOS map of the Fe of RbGd 2 Fe 4 As 4 O 2 , it can be found that it is obviously different from the orbital contribution of RbDy 2 Fe 4 As 4 O 2 and RbTb 2 Fe 4 As 4 O 2 near the Fermi level. It can be found in Figure 3a–c that for RbGd 2 Fe 4 As 4 O 2 , five orbits near the Fermi level contribute, and the main contribution comes from the 3d x z of Fe, that is, 3d x z > 3d y z > 3d x 2 y 2 > 3d z 2 > 3d x y , in which the contribution of the 3d z 2 orbital is higher than the other two materials. It can thus be confirmed that the effect of Gd on 12442 is indeed different from that of Tb and Dy, and the presence of Gd is more conducive to the interlayer electron transfer.

3.1.1. Minus LnO Layers

From the Fermi surface structures of RbLn 2 Fe 4 As 4 O 2 in the NM phase, it can be found that the Fermi surfaces of RbTb 2 Fe 4 As 4 O 2 and the 1111 system are very similar, while the structures of RbDy 2 Fe 4 As 4 O 2 are similar to the 122 system, so it is not difficult to discern the effect of the presence of Ln on these three materials. It is obviously different, so we calculated that if there is no rare earth element, we consider removing the Fermi surface of the LnO layer and the difference map of charge density to obtain the different effects of Ln series elements Gd, Tb, and Dy in the 12442 iron-based superconductors. If the LnO layer is removed, the electronic properties of RbLn 2 Fe 4 As 4 O 2 should be very similar, and the only difference is the structural parameters.
Compared with the energy band structure before removing the LnO layer, the obvious difference in RbLn 2 Fe 4 As 4 O 2 is that the characteristic energy band at X point is moved up, as shown in Figure 4. However, RbLn 2 Fe 4 As 4 O 2 looks very similar in general, indicating that the structural differences caused by the ionic radii have little effect on the properties of these three materials, and the different Ln elements themselves have a major impact on the properties of RbLn 2 Fe 4 As 4 O 2 .

3.1.2. Dy vs. Gd and Tb

Considering that Fe-based compounds and chalcogenides occupy Fe orbitals and that their relationship with properties is the focus of attention, to further study the effects of Gd, Tb, and Dy atoms on Fe orbitals, we calculated the charge densities of RbLn 2 Fe 4 As 4 O 2 in the NM phase. To facilitate the comparison, we have calculated the difference in density ρ (RbDy 2 Fe 4 As 4 O 2 )- ρ (RbGd 2 Fe 4 As 4 O 2 ) and ρ (RbDy 2 Fe 4 As 4 O 2 )- ρ (RbTb 2 Fe 4 As 4 O 2 ). To obtain the charge density of the two materials, both were calculated from the lattice constant of the previous compound. That is to say, the latter retains the lattice parameters and atomic positions of the former, but replaces the lanthanide elements to compare the effects of the two lanthanide elements on the electronic properties of the material.
Figure 5 shows a contour map of this differential density in two planes: the (110) plane containing the LnO layer Fe 2 As 2 intralayer (Fe2As2) and the Fe 2 As 2 interlayer (Fe1As1). The unit of electron density is electrons per unit cell. Positive values indicate areas where Dy attracts charges relative to Gd and Tb. As can be seen from Figure 5, Gd increases the charge in the Fe 3d z 2 -direction orbital relative to Dy, while the charge density in the z-direction of the LnO layer decreases. This indicates that Gd contributes to the transfer of charge from the LnO to the FeAs layer, confirming that previous Gd elements relative to Tb and Dy elements can play a role in enhancing the electron transfer between layers. While Tb reduces the charge in the Fe in-plane directional orbitals relative to Dy, the charge density in the z-direction of the LnO layer decreases. In addition, we can find that the Fe 2 As 2 intralayer and interlayer are different, especially as shown in ρ (RbDy 2 Fe 4 As 4 O 2 )- ρ (RbTb 2 Fe 4 As 4 O 2 ). These differences can be calculated by subtracting the orbital occupancy within the atomic ball and comparing the orbital occupancy numbers of the individual atoms. RbDy 2 Fe 4 As 4 O 2 minus RbTb 2 Fe 4 As 4 O 2 , within the atomic spheres. In units of 10 3 , the differences for Fe are [d z 2 ,d x 2 y 2 ,d x y ,d x z ,d y z ]:[−3.7,−5.1,−0.8,17.3,−16.3]. The net change for Fe is 8.6 × 10 3 electrons. For As, the analogous changes are p x , 2.9 ; p y , 2.9 ; and p z , 1.2 , for a small net change of 6.9 × 10 3 electrons. RbDy 2 Fe 4 As 4 O 2 minus RbGd 2 Fe 4 As 4 O 2 , within the atomic spheres. In units of 10 3 , the differences for Fe are [d z 2 ,d x 2 y 2 ,d x y ,d x z ,d y z ]:[−6.9,−4.2,−0.9,17.8,−15]. The net change for Fe is 9.2 × 10 3 electrons. For As, the analogous changes are p x , 3.6 ; p y , 3.6 ; p z , 1.4 , for a net change of 8.5 × 10 3 electrons. Thus, Tb induces 0.016 electrons into the Fe + As spheres compared to Dy, while Gd induces 0.018 electrons into the Fe + As spheres compared to Tb, with the decrease occurring in the interstitial region or the (Dy, Tb)O or (Dy, Gd)O regions, respectively. This shows that Gd can transfer more electrons from the GdO layer to the FeAs layer than Tb and Dy.

3.2. Ground State: In-Plane SAFM Phase

Here, we discuss magnetic energies and the basic features of the electronic structure of RbLn 2 Fe 4 As 4 O 2 and draw some parallels with related pnictides. The energies of seven of the simplest magnetic configurations relative to the FM1 phase are presented in Table 1, where the Fe atomic sphere moments are also provided. Because the energy of the NM state is much higher than that of the magnetic state, the difference is too large to be conducive to comparison. Therefore, the FM1 ferromagnetic state is used as the reference state. The DFT methods [30] overestimate the Fe moment, and thus the magnetic energy will also be exaggerated, but the relative ordering is still meaningful.
It can be obtained from the table that the ground states of RbLn 2 Fe 4 As 4 O 2 are all in-plane SAFM (AFM1), and their corresponding structures are shown in the Figure 2. Among them, the magnetic moment of Fe is almost 2 μ B , which is very close to the 1.92 μ B magnetic moment of the CaFeAsH [31] ground-state antiferromagnetism calculated by Wien2K before, and the magnetic moment of 2.11 μ B of the CaFeAs 2 [32] ground-state antiferromagnetic force is also very close. The magnetic moments of Fe determined by the neutron scattering experiment for the underdoped samples of LaFeAsO and BaFe 2 As 2 are 2.38 μ B [29] and 1.6 μ B [33], which are greater than the ordered magnetic moment, indicating that the magnetic fluctuation is strong even in the matrix.
Figure 2 shows a magnetic structure diagram; from left to right are FM1, FM2, AFM1, AFM2, AFM3, AFM4, and AFM5 phases. Here, the coordinates of Fe1 and Fe2 are (0, 0, z) and (0.5, 0.5, z), and Ln1 (0.5, 0.5, z) and Ln2 (0, 0, z) are the Ln of different layers. As mentioned above, AFM1 is SAFM, AFM2 is NAFM, and AFM3, AFM4, and AFM5 are different interlayer antiferromagnets. The corresponding AFM1 energy in Table 1 is the lowest, which means that the magnetic ground states of these three materials are in the SAFM state.
Figure 6 shows the Fe pDOS of RbLn 2 Fe 4 As 4 O 2 in the SAFM state. These three materials exhibit distinct electronic structural features in the SAFM state. The five spin-up orbitals are almost fully occupied, with complete electron occupation, while the spin-down orbitals are only partially filled, with an equal distribution of electrons and holes. The Fermi level contribution is mainly concentrated in the spin-down orbitals, where the Fe-3d x z orbital has the largest contribution to the Fermi level, followed by Fe-3d z 2 . Compared to RbLn 2 Fe 4 As 4 O 2 (Ln = Gd and Tb), the contributions of the Fe-3d x z and Fe-3d z 2 orbitals in RbDy 2 Fe 4 As 4 O 2 are reduced by approximately 0.1 eV.
Figure 7 shows the Fe atoms pbands of RbGd 2 Fe 4 As 4 O 2 , RbTb 2 Fe 4 As 4 O 2 and RbDy 2 Fe 4 As 4 O 2 , respectively. The corresponding Fermi level changes in the three materials near the T and Z positions are more obvious. The projected band structure shows that two energy bands, Fe-3d x 2 y 2 and Fe-3d y z and Fe-3d z 2 contribute to the band crossing the Fermi level. From the left to the right of the figure, it can be clearly found that the energy bands near the Fermi level move up sequentially, among which the energy bands of RbGd 2 Fe 4 As 4 O 2 are mainly contributed by Fe-3d y z and Fe-3d z 2 orbitals, while the energy bands of RbDy 2 Fe 4 As 4 O 2 are mainly composed of Fe-3d x 2 y 2 . It can be shown that Gd is more conducive to the electron transfer between layers than Tb and Dy, which is consistent with the conclusion drawn under the nonmagnetic phase.

4. Discussion and Summary

For the first time, we have studied the high-temperature nonmagnetic state and the low-temperature magnetic ground state of three 12442 materials, RbTb 2 Fe 4 As 4 O 2 , RbDy 2 Fe 4 As 4 O 2 , compared them with RbGd 2 Fe 4 As 4 O 2 , and performed a detailed study of the effects of lanthanides and bilayer Fe 2 As 2 . It was found that the two 12442 materials RbTb 2 Fe 4 As 4 O 2 and RbDy 2 Fe 4 As 4 O 2 have the same antiferromagnetic ground state. The ground state of RbLn 2 Fe 4 As 4 O 2 are spin-density-wave-type, in-plane, striped antiferromagnets, and the magnetic moment around each Fe atom is about 2 μ B . We also found that the structural differences caused by the simple ionic radius have little effect on the properties of these three materials. The electronic properties of the materials are strongly influenced by the specific lanthanide element present. It can be confirmed that the effect of Gd on RbLn 2 Fe 4 As 4 O 2 is indeed different from that of Tb and Dy, and the presence of Gd is more conducive to interlayer electron transfer. This means that Gd can transfer more electrons from GdO layer to FeAs layer compared to Tb and Dy. Therefore, RbGd 2 Fe 4 As 4 O 2 has a stronger internal coupling strength of the bilayer Fe 2 As 2 layer. This can explain why the T c of RbGd 2 Fe 4 As 4 O 2 is slightly higher than that of RbTb 2 Fe 4 As 4 O 2 and RbDy 2 Fe 4 As 4 O 2 . We hope that these findings will help to study such materials and understand their relationship with superconducting properties, thus providing theoretical support for searching for superconducting materials with higher T c .

Author Contributions

Conceptualization, Y.-N.H. and D.-Y.L.; Methodology, Y.-N.H., D.-Y.L. and H.-Y.M.; Formal analysis, Y.-N.H. and L.H.; Writing—original draft, Y.-N.H. and H.-P.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China under Grant No. 11904319 and No. 11974354, Zhejiang Provincial Natural Science Foundation of China, under Grant No. LQ18A040002.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Askerzade, I. Unconventional Superconductors; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar] [CrossRef]
  2. Huang, Y.N.; Pickett, W.E. Electronic Coupling between a FeSe Monolayer Film and SrTiO3 Substrate. Phys. Rev. B. 2017, 95, 165107. [Google Scholar] [CrossRef]
  3. Tapp, J.H.; Tang, Z.; Lv, B.; Sasmal, K.; Lorenz, B.; Chu, P.C.W.; Guloy, A.M. LiFeAs: An Intrinsic FeAs-based Superconductor withTc = 18 K. Phys. Rev. B Condens. Matter Mater. Phys. 2008, 78, 060505. [Google Scholar] [CrossRef]
  4. Alireza, P.L.; Ko, Y.T.C.; Gillett, J.; Petrone, C.M.; Cole, J.M.; Lonzarich, G.G.; Sebastian, S.E. Superconductivity up to 29 K in SrFe2As2 and BaFe2As2 at High Pressures. J. Phys. Condens. Matter 2009, 21, 012208. [Google Scholar] [CrossRef]
  5. Wang, F.; Yang, F.; Gao, M.; Lu, Z.Y.; Xiang, T.; Lee, D.H. The Electron Pairing of KxFe2-ySe2. EPL 2011, 93, 57003. [Google Scholar] [CrossRef]
  6. Boeri, L.; Dolgov, O.V.; Golubov, A.A. Is LaFeAsO1-xFx an Electron-Phonon Superconductor? Phys. Rev. Lett. 2008, 101, 026403. [Google Scholar] [CrossRef] [PubMed]
  7. Nitsche, F.; Jesche, A.; Hieckmann, E.; Doert, T.; Ruck, M. Structural Trends from a Consistent Set of Single-Crystal Data of RFeAsO (R = La, Ce, Pr, Nd, Sm, Gd, and Tb). Phys. Rev. B Condens. Matter Mater. Phys. 2010, 82, 134514. [Google Scholar] [CrossRef]
  8. Zhu, X.; Han, F.; Mu, G.; Zeng, B.; Cheng, P.; Shen, B.; Wen, H.H. Sr3Sc2Fe2As2O5 as a Possible Parent Compound for FeAs-based Superconductors. Phys. Rev. B Condens. Matter Mater. Phys. 2009, 79, 024516. [Google Scholar] [CrossRef]
  9. Sefat, A.S.; Singh, D.J.; Ovidiu Garlea, V.; Zuev, Y.L.; McGuire, M.A.; Sales, B.C. Variation of Physical Properties in the Nominal Sr4V2O6Fe2As2. Phys. C Supercond. 2011, 471, 143–149. [Google Scholar] [CrossRef]
  10. Zhi-Cheng, W.; Guang-Han, C. Self-Doped Iron-Based Superconductors with Intergrowth Structures. Acta Phys. Sin. 2018, 67, 207406. [Google Scholar] [CrossRef]
  11. Iyo, A.; Kawashima, K.; Kinjo, T.; Nishio, T.; Ishida, S.; Fujihisa, H.; Gotoh, Y.; Kihou, K.; Eisaki, H.; Yoshida, Y. New-Structure-Type Fe-based Superconductors: CaAFe4As4 (A = K, Rb, Cs) and SrAFe4As4 (A = Rb, Cs). J. Am. Chem. Soc. 2016, 138, 3410–3415. [Google Scholar] [CrossRef]
  12. Liu, Y.; Liu, Y.B.; Tang, Z.T.; Jiang, H.; Wang, Z.C.; Ablimit, A.; Jiao, W.H.; Tao, Q.; Feng, C.M.; Xu, Z.A.; et al. Superconductivity and Ferromagnetism in Hole-doped RbEuFe4As4. Phys. Rev. B. 2016, 93, 214503. [Google Scholar] [CrossRef]
  13. Jiang, H.; Sun, Y.L.; Xu, Z.A.; Cao, G.H. Crystal Chemistry and Structural Design of Iron-Based Superconductors. Chin. Phys. B 2013, 22, 087410. [Google Scholar] [CrossRef]
  14. Ghosh, A.; Ghosh, S.; Ghosh, H. Electron Correlation Induced Orbital Selective Lifshitz Transition in New Hybrid 12442 Iron Based Superconductors. Comput. Mater. Sci. 2020, 183, 109802. [Google Scholar] [CrossRef]
  15. Wu, S.Q.; Wang, Z.C.; He, C.Y.; Tang, Z.T.; Liu, Y.; Cao, G.H. Superconductivity at 33–37 K in ALn2Fe4As4O2(A = K and Cs; Ln = lanthanides). Phys. Rev. Mater. 2017, 1, 044804. [Google Scholar] [CrossRef]
  16. Kirschner, F.K.K.; Adroja, D.T.; Wang, Z.C.; Lang, F.; Smidman, M.; Baker, P.J.; Cao, G.H.; Blundell, S.J. Two-gap superconductivity with line nodes in CsCa2Fe4As4F2. Phys. Rev. B 2018, 97, 060506. [Google Scholar] [CrossRef]
  17. Komlev, A.S.; Gimaev, R.R.; Davydov, A.S.; Zverev, V.I. The influence of chemical impurities on the properties of heavy rare-earth metals (Tb, Dy, Ho): Experimental and theoretical approaches. Materialia 2021, 18, 101166. [Google Scholar] [CrossRef]
  18. Singh, B.; Kumar, P. Density functional study of ACa2Fe4As4F2 (A = K, Rb): Electronic structure, unconventional superconductors. In Proceedings of the 2018 National Conference on Advanced Materials and Nanotechnology (AMN-2018), Noida, India, 15–17 March 2018; AIP Conference Proceedings. Volume 2009, p. 020002. [Google Scholar] [CrossRef]
  19. Wang, Z.C.; He, C.Y.; Wu, S.Q.; Tang, Z.T.; Liu, Y.; Ablimit, A.; Feng, C.M.; Cao, G.H. Superconductivity in KCa2Fe4As4F2 with Separate Double Fe2As2 Layers. J. Am. Chem. Soc. 2016, 138, 7856–7859. [Google Scholar] [CrossRef]
  20. Pavlov, N.S.; Pervakov, K.S.; Nekrasov, I.A. Anomalous metallic oxygen band in the potential superconductor KCa2Fe4As4O2: A DFT study. Comput. Mater. Sci. 2023, 218, 111916. [Google Scholar] [CrossRef]
  21. Wang, Z.; Wang, G.; Tian, X. Electronic structure and magnetism of RbGd2Fe4As4O2. J. Alloys Compd. 2017, 708, 392–396. [Google Scholar] [CrossRef]
  22. Wang, Z.C.; He, C.Y.; Wu, S.Q.; Tang, Z.T.; Liu, Y.; Cao, G.H. Synthesis, Crystal Structure and Superconductivity in RbLn2Fe4As4O2 (Ln = Sm, Tb, Dy, and Ho). Chem. Mater. 2017, 29, 1805–1812. [Google Scholar] [CrossRef]
  23. Blaha, P.; Schwarz, K.; Tran, F.; Laskowski, R.; Madsen, G.K.H.; Marks, L.D. WIEN2k: An APW+lo Program for Calculating the Properties of Solids. J. Chem. Phys. 2020, 152, 074101. [Google Scholar] [CrossRef] [PubMed]
  24. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  25. Ameri, I.; Boularaf, A.; Drief, F.; Zaoui, A.; Kacimi, S. Ab-initio study of magnetic and electronic properties of the perovskites RFeO3: 4f-R valence electrons effects. J. Magn. Magn. Mater. 2021, 537, 168214. [Google Scholar] [CrossRef]
  26. Söderlind, P.; Turchi, P.; Landa, A.; Lordi, V. Ground-state properties of rare-earth metals: An evaluation of density-functional theory. J. Phys. Condens. Matter 2014, 26, 416001. [Google Scholar] [CrossRef]
  27. Kanoun, M.B.; Reshak, A.H.; Kanoun-Bouayed, N.; Goumri-Said, S. Evidence of Coulomb correction and spin-orbit coupling in rare-earth dioxides CeO2, PrO2 and TbO2: An ab initio study. J. Magn. Magn. Mater. 2012, 324, 1397–1405. [Google Scholar] [CrossRef]
  28. Jamnezhad, H.; Jafari, M. Structural, electronic, and optical properties of C-type Gd2O3: A density functional theory investigation. J. Comput. Electron. 2017, 16, 272–279. [Google Scholar] [CrossRef]
  29. Ma, F.; Lu, Z.Y. Iron-based layered compound LaFeAsO is an antiferromagnetic semimetal. Phys. Rev. B 2008, 78, 033111. [Google Scholar] [CrossRef]
  30. Ortenzi, L.; Gretarsson, H.; Kasahara, S.; Matsuda, Y.; Shibauchi, T.; Finkelstein, K.D.; Wu, W.; Julian, S.R.; Kim, Y.J.; Mazin, I.I.; et al. Structural Origin of the Anomalous Temperature Dependence of the Local Magnetic Moments in the CaFe2As2 Family of Materials. Phys. Rev. Lett. 2015, 114, 047001. [Google Scholar] [CrossRef]
  31. Huang, Y.N.; Liu, D.Y.; Zou, L.J.; Pickett, W.E. Role of Hydrogen in the Electronic Properties of CaFeAsH-based Superconductors. Phys. Rev. B. 2016, 93, 195148. [Google Scholar] [CrossRef]
  32. Huang, Y.N.; Yu, X.L.; Liu, D.Y.; Zou, L.J. Magnetism and Electronic Structures of Novel Layered CaFeAs2 and Ca0.75(Pr/La)0.25FeAs2. J. Appl. Phys. 2015, 117, 17E113. [Google Scholar] [CrossRef]
  33. Shimojima, T.; Ishizaka, K.; Ishida, Y.; Katayama, N.; Ohgushi, K.; Kiss, T.; Okawa, M.; Togashi, T.; Wang, X.Y.; Chen, C.T.; et al. Orbital-Dependent Modifications of Electronic Structure across the Magnetostructural Transition in BaFe2As2. Phys. Rev. Lett. 2010, 104, 057002. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (Color online) The structures of RbLn 2 Fe 4 As 4 O 2 (Ln = Gd, Tb, and Dy) in NM phase, including the structural diagrams of the lanthanide layer and the iron atomic layer.
Figure 1. (Color online) The structures of RbLn 2 Fe 4 As 4 O 2 (Ln = Gd, Tb, and Dy) in NM phase, including the structural diagrams of the lanthanide layer and the iron atomic layer.
Materials 16 04123 g001
Figure 2. (Color online) The structures of RbLn 2 Fe 4 As 4 O 2 (Ln = Gd, Tb, and Dy) in FM1, FM2, AFM1, AFM2, AFM3, AFM4, and AFM5 phases.
Figure 2. (Color online) The structures of RbLn 2 Fe 4 As 4 O 2 (Ln = Gd, Tb, and Dy) in FM1, FM2, AFM1, AFM2, AFM3, AFM4, and AFM5 phases.
Materials 16 04123 g002
Figure 3. (Color online) The Fe projection of bands, projected density of states, and the Fermi surfaces of (a,d) RbGd 2 Fe 4 As 4 O 2 , (b,e) RbTb 2 Fe 4 As 4 O 2 , and (c,f) RbDy 2 Fe 4 As 4 O 2 in NM phase.
Figure 3. (Color online) The Fe projection of bands, projected density of states, and the Fermi surfaces of (a,d) RbGd 2 Fe 4 As 4 O 2 , (b,e) RbTb 2 Fe 4 As 4 O 2 , and (c,f) RbDy 2 Fe 4 As 4 O 2 in NM phase.
Materials 16 04123 g003
Figure 4. Electronic band structure along high symmetry line Γ X M Γ Z of (a) RbGd 2 Fe 4 As 4 O 2 minus GdO layer, (b) RbTb 2 Fe 4 As 4 O 2 minus TbO layer, and (c) RbDy 2 Fe 4 As 4 O 2 minus DyO layer in NM phase.
Figure 4. Electronic band structure along high symmetry line Γ X M Γ Z of (a) RbGd 2 Fe 4 As 4 O 2 minus GdO layer, (b) RbTb 2 Fe 4 As 4 O 2 minus TbO layer, and (c) RbDy 2 Fe 4 As 4 O 2 minus DyO layer in NM phase.
Materials 16 04123 g004
Figure 5. Contour plots of the density difference ρ (RbDy 2 Fe 4 As 4 O 2 )- ρ (RbGd 2 Fe 4 As 4 O 2 ) and ρ (RbDy 2 Fe 4 As 4 O 2 )- ρ (RbTb 2 Fe 4 As 4 O 2 ). The large red region denotes where contours have been cut off due to the large and meaningless difference of the RbDy 2 Fe 4 As 4 O 2 and Rb(Gd/Tb) 2 Fe 4 As 4 O 2 densities in NM phase.
Figure 5. Contour plots of the density difference ρ (RbDy 2 Fe 4 As 4 O 2 )- ρ (RbGd 2 Fe 4 As 4 O 2 ) and ρ (RbDy 2 Fe 4 As 4 O 2 )- ρ (RbTb 2 Fe 4 As 4 O 2 ). The large red region denotes where contours have been cut off due to the large and meaningless difference of the RbDy 2 Fe 4 As 4 O 2 and Rb(Gd/Tb) 2 Fe 4 As 4 O 2 densities in NM phase.
Materials 16 04123 g005
Figure 6. For RbLn 2 Fe 4 As 4 O 2 with SAFM order: (a) RbGd 2 Fe 4 As 4 O 2 , (b) RbTb 2 Fe 4 As 4 O 2 , and (c) RbDy 2 Fe 4 As 4 O 2 projected DOS for the Fe atom. Note the small DOS at the Fermi energy (zero of energy).
Figure 6. For RbLn 2 Fe 4 As 4 O 2 with SAFM order: (a) RbGd 2 Fe 4 As 4 O 2 , (b) RbTb 2 Fe 4 As 4 O 2 , and (c) RbDy 2 Fe 4 As 4 O 2 projected DOS for the Fe atom. Note the small DOS at the Fermi energy (zero of energy).
Materials 16 04123 g006
Figure 7. For RbLn 2 Fe 4 As 4 O 2 with SAFM order: (a) RbGd 2 Fe 4 As 4 O 2 , (b) RbTb 2 Fe 4 As 4 O 2 , and (c) RbDy 2 Fe 4 As 4 O 2 projected band structure for the Fe atom. Upper row: spin up. Lower row: spin down.
Figure 7. For RbLn 2 Fe 4 As 4 O 2 with SAFM order: (a) RbGd 2 Fe 4 As 4 O 2 , (b) RbTb 2 Fe 4 As 4 O 2 , and (c) RbDy 2 Fe 4 As 4 O 2 projected band structure for the Fe atom. Upper row: spin up. Lower row: spin down.
Materials 16 04123 g007
Table 1. Total energy difference, for RbLn 2 Fe 4 As 4 O 2 (Ln = Gd, Tb, and Dy) of seven magnetic phases including FM1, FM2, AFM1, AFM2, AFM3, AFM4, and AFM5 phases. The reference is the energy of the FM1 phase ( Δ E = E ( S ) F M E F M 1 ). The corresponding magnetic moment in the Fe sphere is given.
Table 1. Total energy difference, for RbLn 2 Fe 4 As 4 O 2 (Ln = Gd, Tb, and Dy) of seven magnetic phases including FM1, FM2, AFM1, AFM2, AFM3, AFM4, and AFM5 phases. The reference is the energy of the FM1 phase ( Δ E = E ( S ) F M E F M 1 ). The corresponding magnetic moment in the Fe sphere is given.
Magnetic StructureFM1FM2AFM1AFM2AFM3AFM4AFM5
(a) Relative energy (meV)RbGd 2 Fe 4 As 4 O 2 0−11−429−239−6−27−8
RbTb 2 Fe 4 As 4 O 2 0−81−478−218156−44−6
RbDy 2 Fe 4 As 4 O 2 0−872−1236−1083−879−874−877
(b) Fe moment ( μ B )RbGd 2 Fe 4 As 4 O 2 0.800.832.022.070.850.810.84
RbTb 2 Fe 4 As 4 O 2 0.810.831.982.060.900.840.87
RbDy 2 Fe 4 As 4 O 2 0.910.841.912.030.860.840.90
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Huang, Y.-N.; Liu, D.-Y.; Mei, H.-Y.; Han, L.; Yang, H.-P. Role of Lanthanides and Bilayer Fe2As2 in the Electronic Properties of RbLn2Fe4As4O2 (Ln = Gd, Tb, and Dy) Superconductors. Materials 2023, 16, 4123. https://doi.org/10.3390/ma16114123

AMA Style

Huang Y-N, Liu D-Y, Mei H-Y, Han L, Yang H-P. Role of Lanthanides and Bilayer Fe2As2 in the Electronic Properties of RbLn2Fe4As4O2 (Ln = Gd, Tb, and Dy) Superconductors. Materials. 2023; 16(11):4123. https://doi.org/10.3390/ma16114123

Chicago/Turabian Style

Huang, Yi-Na, Da-Yong Liu, Hong-Ying Mei, Li Han, and Huan-Ping Yang. 2023. "Role of Lanthanides and Bilayer Fe2As2 in the Electronic Properties of RbLn2Fe4As4O2 (Ln = Gd, Tb, and Dy) Superconductors" Materials 16, no. 11: 4123. https://doi.org/10.3390/ma16114123

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop