Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 14 March 2024
Sec. Condensed Matter Physics

Theoretical investigation of the Co-occurrence of superconductivity and antiferromagnetism in iron-based high-temperature superconductors

  • 1Department of Physics, College of Natural and Computation Sciences, Wolkite University, Wolkite, Ethiopia
  • 2Department of Physics, College of Natural and Computation Sciences, Dilla University, Dilla, Ethiopia

In this research work, the Co-occurrence of superconductivity and antiferromagnetism in Ca0.741La0.261Fe1xCoxAs2 iron-based high-temperature superconductors have been investigated. Based on the multi-band nature of Ca0.741La0.261Fe1xCoxAs2 iron-based superconductors, a two-band model Hamiltonian that contains intra- and inter-bands is developed. By employing the matrix form of the temperature-dependent Green’s function formalism with the two-band Hamiltonian, the mathematical expressions of superconducting transition temperature as functions of the superconducting order parameter and antiferromagnetism translational temperature as a function of the magnetic order parameter are obtained, respectively. The plotted graph of the superconducting and magnetic temperature as a function of the magnetic order parameter indicates a clear possibility of Co-occurrence of superconductivity and the antiferromagnetism order parameter in Ca0.741La0.261Fe1xCoxAs2 iron-based superconductors in the range of magnetic order parameters between 2.3 meV and 6.07 meV, which is in good agreement with experimental observations. This research contributes to understanding the complex behavior of high-temperature superconductors and provides valuable technological applications for other fields.

1 Introduction

Current research in the field of condensed matter physics is dedicated to characterizing the properties of materials based on their structure and electronic behavior. While most metal systems have weaker Coulomb interaction energies compared to electron kinetic energies, the interaction between electrons, whether direct or indirect, significantly influences the physical properties of strongly correlated electron systems. These systems include unconventional superconductors, Mott insulators, and heavy fermions. On the other hand, superconductivity is a complex phenomenon that requires extensive experimental and theoretical efforts across a wide range of materials such as oxides, magnetic compounds, and organic compounds with their wide range of modern technological and industrial applications like MRI, Maglev train, and electronic device. This particular subfield remains one of the most challenging areas to understand fully.

The phenomenon of superconductivity, initially observed in high-purity mercury by H. Kamerlingh Onnes in April 1911, is characterized by a sudden drop in electrical resistance to zero at 4.2 K [1, 2]. This discovery led to the subsequent identification of elemental and binary superconductors with critical temperature (Tc) values up to 23.2 K [3, 4], which remained unbeaten for over 50 years. Extensive experimental and theoretical efforts were made to understand the underlying microscopic mechanisms and the perfect diamagnetism exhibited by superconductors. In 1933, W. Meissner and R. Ochsenfeld discovered that a superconductor, when cooled in the presence of a static magnetic field, expels the magnetic field from its interior [36]. The theoretical breakthrough came in 1957 with the formulation of the BCS theory by John Bardeen, Leon Neil Cooper, and John Robert Schrieffer, which provided a microscopic explanation for superconductivity [68]. In the late 1970s and early 1980s, superconductivity was discovered in heavy fermion and nearly magnetic systems despite their lower critical temperatures [8]. Researchers began seeking new pairing interactions to eventually achieve high-temperature superconductivity. The pivotal moment came in 1986 when J. G. Bednorz and K. A. Muller discovered La(1-x) Bax CuO4 with a Tc of 30 K, sparking intense and continuous research in the field of high-temperature superconductivity [4, 911].

The emergence of LaO(1-x) Fx FeAs, a novel superconductor with a transition temperature of 26 K, marked a significant milestone in superconductivity research [11]. This breakthrough led to a shift in focus from high-temperature cuprates (HTSC or CuSC) to iron-based superconductors (FeSC), as evidenced by the redirection of researchers and funding. Observed as the second family of unconventional High-Tc superconductors, the FeSC exhibits several similarities to cuprate superconductors. These similarities include a layered crystal structure, relatively high critical temperature (Tc), superconductivity occurring near a magnetic phase induced by parameter adjustments like chemical doping or pressure, and notably, a highly comparable temperature-concentration (T - x) phase diagram. Moreover, it is commonly accepted that electronic conduction in the FeAs/FeSe layers is linked, while cuprates have charge carriers that are delocalized in the basal copper oxide planes [11]. However, there are some significant differences: in the edge-sharing tetrahedron, the pnictogen anions are arranged above and below the Fe plane instead of at the same height as in the copper oxide plane [12], and the superconductivity of cuprates is derived from doping a Mott insulator. However, the parent compound of FeSC is a poor metal due to partial gapping around the Fermi surface (FS) at low temperatures [1315], which does not initiate an insulating state. Doping the two-dimensional (2D) copper oxide plane of cuprate superconductors causes a sharp decline in Tc. The cuprates exhibit a predominantly s-wave gap symmetry, whereas the FeSC are generally immune to this effect [14]. This change is reflected in the early reviews, and FeSC research is now a major area of condensed matter physics, as evidenced by the more than 3000 citations that show its active pursuit [1618].

The microscopic pairing process for superconductivity has been the subject of extensive research since the discovery of Fe-pnictides, a novel class of high-temperature superconductors [12, 1922]. However, still, the microscopic pairing process of superconductivity (SC) has to be the main issue in unconventional high-temperature superconductors. The main objective of this theoretical study of the co-occurrence of antiferromagnetism (AFM) and superconductivity in a Co-doped Ca0.73La0.27FeAs2 (Co-CaLa112) high-temperature superconductor is to provide insight on the superconductivity pairing process. FeSC holds great interest for several compelling reasons. Firstly, it offers the opportunity to explore fascinating physics arising from the Co-occurrence of superconductivity and magnetism. Secondly, the wide range of compounds available for study, coupled with the multi-band electronic structure, holds the potential for uncovering the elusive mechanism behind high-temperature superconductivity and discovering methods to enhance it (Tc). Lastly, FeSC shows promise for practical applications due to their higher critical field (Hc) compared to cuprates, as well as their strong and isotropic critical currents. Based on these reasons make them appealing for electrical power and magnetic applications, while the Co-occurrence of magnetism and superconductivity makes them interesting for spintronic [23, 24]. Therefore, FeSC offers a valuable opportunity to explore the influence of structural and electronic factors on the physical properties and pairing mechanism of high-temperature superconductivity (Tc). In the study of FeSC, there is a significant focus on investigating the interplay between antiferromagnetism (AFM) and superconductivity (SC). This attention is driven by the belief that spin fluctuations play a crucial role in the pairing mechanism [68, 25].

Since the groundbreaking discovery of LaO(1-x) Fx FeAs (1111-family), a novel superconductor with a transition temperature of 26 K [11], significant progress has been made in the field of superconductivity research in Fe-pnictides. Includes the 1111-family REFeAs (O, F) (RE = rare-earth elements) [19, 20], the 122 families AeFe2 As2 (Ae = alkaline earth metals such as Ca, Sr, Ba) [21, 22], the 111-family AFeAs (A = alkali metals like Li, Na) [26, 27], the 11 family-FeSe [28, 29], and the 10-n-8 family Ca10 (Ptn As8) (Fe2 As2)5 (A = Pt, Pd, Ir; n = 3, 4) [27]. Obtaining a deeper understanding of the superconductivity mechanism in iron-based compounds and striving for advancements in Tc (critical temperature) necessitate the crucial discovery of new families of iron-based superconductors. Recently, the 112 family of iron-based superconductors, specifically the Ca(1-x) Ax FeAs2 (A = Rare Earth metal, such as La, Ce, Pr, Nd, etc.), has been experimentally confirmed [30, 31]. From this 112 family, the Ca(1-x) Lax FeAs2 (CaLa112) is one of the parent compounds with Tc up to 43K, crystalizes in a monoclinic lattice with the FeAsCa/LaAsCa/LaFeAs layer stacking [28]. CaLa112 is distinctive in several aspects due to the presence of zig-zag chains made of As layers alongside the prototypical FeAs layers consisting of edge-sharing FeAs4 tetrahedra. This combination of As chains and FeAs layers sets CaLa112 apart from other crystals. By introducing electron over doping in the parent compound of CaLa112 FeSC, the dual nature of moveable and local magnetism in FeSC is exemplified. This electron-doped CaLa112 FeSC undergoes a structural phase transition from monoclinic to triclinic at 58K, while a paramagnetic to stripe antiferromagnetic phase transition occurs at 54 K [32]. In addition to hole-like carriers introduced by Ca doping, electron-like carriers are incorporated through Co substitution on the iron sites, which contributes to the stabilization of superconductivity in Ca 0.73 La 0.27 FeAs2 [32]. The structure and magnetic phase transitions in Co-doped Ca0.73La0.27 FeAs2 (Co-CaLa112) are suppressed, leading to the emergence of bulk superconductivity with a critical temperature (Tc) of up to 20 K [30]. Doping experiments reveal the microscopic coexistence of AFM and SC in Co-doped samples with doping concentrations of x = 0.025 and 0.033 [30]. To explore the microscopic coexistence of AFM and SC in the 112 material, we utilize a simplified two-band model commonly used in the study of the relation between AFM and SC of FeSC [3335]. This model consists of a hole pocket located at the center of the Brillouin zone (BZ) and an electron pocket situated at the corner of the Brillouin zone as shown in Figure 1 Refs. [32, 33]. By studying the interplay between these phenomena, their connections can be uncovered, providing an empirical foundation for developing a comprehensive theoretical model.

FIGURE 1
www.frontiersin.org

FIGURE 1. Left: modified electronic band dispersion of the two-band model reflected in this paper, in the unfolded Brillouin zone. The circular hole-like FS is in the center, with SC order parameter 1, and the elliptical electron-like FSs are at =0,πand=π,0 with SC order parameter 2. The magnetic order with momentum Q=π,0 hybridizes hole and electron FSs separated by Q but leaves FSs at =±π,0 intact. Right: by doping one may adjust the size and shape of hole and electron bands, and also magnetic order parameter can be incommensurate, with momentum Q+q Ref. [33].

This study, which is based on an experimental perspective, focuses on the theoretical analysis of the interplay between superconductivity and antiferromagnetism in Ca (0.74 (1)) La (0.26 (1)) (Fe(1-x) Cox) As (2) iron-based high-temperature superconductors by considering intra and inter bands that the system incorporates two computing phenomena involving electron-hole-like pairing and electron-electron pairing. The mathematical expression for the superconducting order parameter, magnetic order parameter, superconducting transition temperature, and AFM transition temperature is found using the matrix form of the temperature-dependent Green’s function formalism with the two-band model Hamiltonian. Our research’s findings determine that AFM and SC are established over specific order parameter ranges, which will give clues for the mystery of high-temperature superconductors.

2 Formulation of the problem

Clarifying the origin of superconductivity in iron-based superconductors requires an understanding of their phase diagram. Spin density wave (SDW) with antiferromagnetism (AFM) is observed in parent compounds of iron-based superconductors below the Neel temperature. Superconductivity arises when hole or electron doping is applied, suppressing magnetization. The close relationship between AFM and SC phases suggests that spin fluctuation mediates the formation of a Cooper pair, leading to s+-wave order [34], where the gap function sign on hole Fermi surfaces (FSs) centered at momentum k = (0,0) is opposite to that on electron FSs at k = (0,π) and (π,0). Conversely, it is suggested that superconductivity arises from orbital fluctuation and takes the form of s++ wave order, where the two signs are the same [36]. s+-and s++ are possible candidates for SC symmetry in iron-based superconductors. Given the proximity of the AFM and SC phases, insights into the superconductivity of iron-based superconductors may be gained from examining their boundary.

Numerous experimental and computational investigations employing Density Functional Theory (DFT) have provided evidence that the Fermi surface of iron-based superconductors comprises two hole surfaces near the Γ point and two electron surfaces near the M point within the Brillouin zone of the Fe/cell [34, 35]. Additionally, these studies highlight the multi-band character of the band structure of FeScs. Calculations of the band structure have further revealed that the Fe 3d orbitals make significant contributions to the spectral weight in the vicinity of the Fermi energy. Neutron diffraction, X-ray diffraction, and NMR investigations all support the microscopically coexisting phase of AFM and SC orders in the Co-CaLa112 high-temperature superconductor [30]. An exploration of the coexistence phase has been conducted theoretically [34]. In the Co-CaLa112 high-temperature superconductor phase, it is generally believed that paramagnetic FSs are made up of a nearly circular hole pocket located at k= (0,0), along with an elliptical electron pocket at k=(π,0) and (0, π) as illustrated the model in Figure 1. The introduction of AFM order with a wave vector Q=(π,0) combines the hole and electron dispersions, resulting in the opening of an SDW gap [32]. There is an additional electron pocket at the edge of the Brillouin zone and a hole-like pocket at the center of the doped Brillouin zone, as observed experimentally in Fe pnictides, in addition to the appearance of an elliptical electron pocket at the corner and a circular hole pocket at the center. Since the fundamental properties of the SC and antiferromagnetic interactions and their interplay should not be significantly affected by the number of bands, we investigate a simple model with one circular hole and one elliptical electron band. In the case of pnictides, the twofold degeneracy of hole and electron states at the center and corners of the Brillouin zone is omitted since it seems to not affect magnetic order or superconducting [3741]. By utilizing 4 × 4 matrix representations of single particle normal and anomalous thermal Green’s functions within a two-band model Hamiltonian, along with appropriate mean-field approximations, the simultaneous occurrence of superconductivity and antiferromagnetism becomes a highly plausible scenario. This section focuses on investigating the coexistence of superconductivity and antiferromagnetism in the multi-band iron-based superconductor Ca (0.74 (1)) La (0.26 (1)) (Fe(1-x) Cox) As (2), as well as examining the influence of magnetic ordering on both the superconducting order parameter ¯ and the transition temperature (Tc) within the framework of the multi-band model Hamiltonian. Based on the multi-band nature of Ca0.741La0.261Fe1xCoxAs2 iron-based superconductors, a two-band model Hamiltonian which contains hole and electron bands is developed. The two-band model Hamiltonian includes the free fermion part H0, and the fermion interactions in superconducting and antiferromagnetic channels,

H=H0+H+HAFM

The free fermion channel of the Hamiltonian is

H0=kσε1kc1kσ+c1kσ+kσε2kc2kσ+c2kσ

Where ε1kandε2k measured from the Fermi energy μchemicalpotential,i.e,εik=εikμi=1,2, the k are momentum measured from center of BZ and k are deviations from Q0, we assume an inversion symmetry ε1,2k=ε1,2k.

The interaction that involves the exchange of electron pairs between the hole and electron pockets is the dominant term in the pairing interaction. However, there exist several other pair-scattering interactions that also contribute to the overall pairing mechanism [40, 42]. The pairing interactions in the superconducting channel for band 1 is

H¯1k=kkV11SCk,k`c1k+c1k+c1kc1kkkV12SCk,kc1k+c1k+c2kc2k+c2k+c2k+c1kc1k

For band 2 is

H¯2k=kkV22SCk,k`c2k+c2k+c2kc2k12kkV21SCk,kc2k+c2k+c1kc1k+c+1kc1k+c2kc2k

We introduce a mean-field equation with order parameter in a self-consistent manner, in the superconducting state for band 1 and 2 as stated by;

11=kkV11SCk,k`c1kc1kand12=kkV12SCk,k`c2kc2k
22=kkV22SCk,k`c2kc2kand21=kkV21SCk,k`c1kc1k

where the 11 and 12 are intra and inter-band superconducting order parameter for band 1, and 22 and 21 are intra and inter-band superconducting order parameter for band 2, respectively

To consider the AFM state, we assume the following intra and interband Coulomb interactions:

HAFM=12kkα1,2AFMk,kc1k+Q+c1kc1k+c1k+Q+c2k+c2kc2k+c2k+Q12kkα12AFMk,kc2k+Q+c1kc+1kc2k+Q+c1k+c2k+Qc2k+Q+c1k

where the first and second terms are intra-band and inter-band Coulomb repulsions, respectively. α1,2AFM and α12AFM are the parameters of intra and inter-band Coulomb interaction, respectively. We solve a mean-field equation with order parameter in a self-consistent manner, in the AFM state with ordering vector Q, as stated by;

intraAFM=kkα11AFMk,k`c1kc1k+Q=kkα22AFMk,k`c2kc2k+Q

And

interAFM=kkα12AFMk,k`c1kc2k+Q=kkα21AFMk,k`c2k+Qc1k

where intraAFM and interAFM are intra-band and inter-band AFM order parameters, respectively, and take wave vector Q=(π,0). Note that in this calculation the interband AFM order parameter is neglected because very small and insignificant [34, 43].

The two-band model mean field Hamiltonian can be rewritten as;

H=kσε1kc1kσ+c1kσ¯1kkkc1k+c1k++kσε2kc2kσ+c2kσ¯2kkkc2k+c2k+¯Mkkkc1k+Q+c1k++c1k+c2k+Q++h.c(1)

Where ε1k and ε2k are the quasiparticle energies for bands 1 and 2, respectively, operators cikσ+cikσ annihilation (creation) or operator for hole and electron bands, respectively, and spin-σ electron of momentum- k and band i=1,2, summation represents sum over all the -k and h.c is the hopping parameters which are insignificant for this superconducting pairing terms. Within the fermionic basis, there are numerous approaches to decoupling the interacting Hamiltonian, we use mean-field approximation to decouple our assumed interacting orbitals to develop the mean-field Hamiltonian described in Eq. 1 [44], and to define order parameters. The mean field approximation assumes that the average interaction can be approximated by their average values, ignoring the individual fluctuations beyond the mean field, which helps to investigate the stability of superconductivity and magnetism in HTSC specifically in FeBSC, and understand the interplay between these two phenomena. As is typically practiced in mean-field theory, one substitutes a specific operator with its average value multiplied by a small fluctuating term. To provide a concrete example, consider the bilinear term, cik+cik+cik`cik` and replaced by cik`cik`c+ikcik+. From Eq. 1, the following effective order parameters are obtained. Thus,

¯1k=11+12=kkV11SCk,k`c1kc1k+kkV12SCk,k`c2kc2k(2)
¯2k=22+21=kkV22SCk,k`c2kc2k+kkV21SCk,k`c1kc1k(3)
¯Mk=kkα11AFMk,k`c1kc1k+Q+kkα12AFMk,k`c1kc2k+Q(4)

Where ¯1k¯2k are the effective superconducting order parameters for bands 1 and 2, respectively, with their first and second terms are intra-band and inter-band Coulomb repulsions, respectively and ¯Mk is antiferromagnetic exchange interaction for intra and inter-band pair hopping. In the following calculation we use the terms as; V11SC = V11, V12SC = V12, V22SC = V22, V21SC = V21, and α11AFM = α11, α12AFM = α12. Each band has its proper pairing interaction. V11 and V22. While pair interchanges between the two bands are assured by V12 = V21 term.

By introducing the τ dependent quantum operators to study the system, we need to transform to the Heisenberg picture [33].

ckστ=eτHckσeτH(5)
ckσ+τ=eτHckσ+eτH(6)

This satisfies the Heisenberg equation of motion.

ckσττ=H,ckστ(7)
ckσ+ττ=H,ckσ+τ(8)

Inserting Eq. 1 into the right-hand side of the equation Eq. 7 and Eq. 8 we obtain the system of differential equations for the time-dependent quantum operator [45].

c1kττ=ε1kc1k+¯1kc1k+(9)
c2kττ=ε2kc2k+¯2kc2k++¯Mkc1k+(10)
c1k+ττ=ε1kc1k++¯1k+c1k++¯Mk+c2+Qk+(11)
c2k+ττ=ε2kc2k++¯2k+c2k+(12)

Using Eqs 912 can drive a system of differential equations for the thermal, Green function techniques and solution for the equation of motion for both operators. In the state of superconductors, the Nambu-Gorkov formalism characterizes it as a state where symmetry is broken [46, 47]. To grip these equations handily, Nambu proposed the following a four-component space spanned by the four-component operators;

akτ=c1kτc1k+τc2kτc2k+τ(13)

And its corresponding conjugate

ak+τ,=c1k+τ,c1kτ,c2k+τ,c2kτ,(14)

Where akτanda+kτ, are nowadays commonly called NambuGorkov operators.

2.1 Green’s functions

From the NambuGorkov operators, we define the single particle 4×4 matrix of Green’s functions in Nambu space.

GTSk,ττ,=Takτak+τ,(15)

Where T is the time ordering operator and the subscript TS shows the elements of the 4 × 4 matrices by substituting Nambu operators the explicit yields,

GTSk,ττ,=Tc1kτc1k+τ,Tc1kτc1kτ,Tc1kτc2k+τ,Tc1kτc2kτ,Tc1k+τc1k+τ,Tc1k+τc1kτ,Tc1k+τc2k+τ,Tc1k+τc2kτ,Tc2kτc1k+τ,Tc2kτc1kτ,Tc2kτc2k+τ,Tc2kτc2kτ,Tc2k+τc1k+τ,Tc2k+τc1kτ,Tc2k+τc2k+τ,Tc2k+τc2kτ,(16)

To study the physical properties, we must define the following thermal Green’s functions;

Gijk,ττ,=Tcikτcjk+τ,(17)
Fijk,ττ,=Tcikτcjkτ,(18)
Fij+k,ττ,=Tcik+τcjk+τ,(19)

And

GijTk,ττ,=Tcik+τcjkτ,(20)

The correlation functions Gijk,ττ, and Fijk,ττ, are the normal and anomalous one-particle Green functions, which are elements of the Green function matrix GTS. i,j indicates the band index 1 or 2. Fij+k,ττ, implies the complex conjugate of Fij. We will often suppress the spin index. The anomalous green functions, Fij+k,ττ,, and GijTk,ττ, assumed that singlet pairing. As the creation and annihilation operators of electrons are fermionic, and with these specifications, the Green function matrix GTSk,ττ, can be written in a more compact form;

GTSk,ττ,=G11k,ττ,F11k,ττ,G12k,ττ,F12k,ττ,F11+k,ττ,G11Tk,ττ,F21+k,ττ,G12Tk,ττ,G21k,ττ,F21k,ττ,G22k,ττ,F22k,ττ,F12+k,ττ,G21Tk,ττ,F22+k,ττ,G22Tk,ττ,(21)

Where T,s = 1,2,3,4, and G11k,ττ,, G22k,ττ, are the intraband Green’s functions for electrons of up spin, G11Tk,ττ,, G22Tk,ττ,, are the intraband Green’s functions for holes with spin down while G12k,ττ,, G21k,ττ,, are interpreted to be interband Green’s functions for electrons of up spin and G12Tk,ττ,,G21Tk,ττ, are taken to be interband Green’s functions for holes with spin down. F11k,ττ,, F22k,ττ, are the anomalous intraband Green’s functions while F12k,ττ,, F21k,ττ, represent the interband anomalous thermal Green. s functions involving electrons in different bands and F11+k,ττ,, F22+k,ττ,, are the complex conjugate of the anomalous intraband thermal Green’s function, while F12+k,ττ,,F21+k,ττ, are the complex conjugate of the anomalous interband thermal Green’s function.

In the four-component language, the sixteen equations of motion for the propagators lead to the energy matrix equation.

KOτGTSk,ττ,=δOττ,(22)

Where the operator KOτ is 4x4 energy matrix and δOττ, is 4x4 unit matrix given by

KOτ=τε1k¯1k00¯1k+τ+ε1k¯Mk+00¯Mkτε2k¯2k00¯2k+τ+ε2k(23)

And

δOτ=δOττ,0000δOττ,0000δOττ,0000δOττ,(24)

Eq. 22 is obtained from the results of Eqs. 912 by differentiating τ order products. We use the common way of Fourier transforming the correlation function, transitioning it from k-space to momentum space according to its established definition [48];

Gijk,ττ,=1βωnеiωnτGijp,iωn,Gijp,iωn=0βdτеiωnτGijk,ττ,(25)
Fijk,ττ,=1βωnеiωnτFijp,iωn,Fijp,iωn=0βdτеiωnτFijk,ττ,(26)

and

GijTk,ττ,=1βωnеiωnτGTijp,iωn,GijTp,iωn=0βdτеiωnτGijTk,ττ,(27)

where ωn=2n+1π/β, is Matsubara frequencies, β=1KBT with KB as the Boltzmann constant, T is the temperature, summation represent sum over all the ωn, and i,j show the band index. These transformations to momentum-dependent make the mathematics easily manageable. After performing the Fourier transformations of the correlation function and for the τ=τ,, the matrix product in Eq. 22 becomes;

iωnε1k¯1k00¯1k+iωn+ε1k¯Mk+00¯Mkiωnε2k¯2k00¯2k+iωn+ε2k×G11p,iωnF11p,iωnG12p,iωnF12p,iωnF11+p,iωnG11Tp,iωnF21+p,iωnG12Tp,iωnG21p,iωnF21p,iωnG22p,iωnF22p,iωnF12+p,iωnG21Tp,iωnF22+p,iωnG22Tp,iωn=1000010000100001(28)

From the matrix Eq. 28, we obtain the inverse of the Fourier transformed 4 × 4 Green’s function matrix in momentum space as

GT,S1p,iωn=iωnε1k¯1k00¯1k+iωn+ε1k¯Mk+00¯Mkiωnε2k¯2k00¯2k+iωn+ε2k(29)

and gives the components of the inverse Green function for ¯Mk=M where magnetic order parameter and suppressing k. To study the order parameters, therefore the equation of motions is;

F11p,iωn=c1kc1k=1Dωn¯1ωn2+ε22+¯22+M2¯2(30)
F12p,iωn=c1kc2k+Q=1DωnMEKε2iωn+ε1+¯1¯2MM3(31)

And

F22p,iωn=c2kc2k=1Dωn¯2ωn2+ε12+¯12M2¯1(32)
F21p,iωn=c2kc1k=1DωnMEKε1iωn+ε2+¯1¯2MM3(33)
Dωn=ωn2+ε12+¯12ωn2+ε22+¯22+2M2ωn2+2M2ε1ε2+M42M2¯1¯2(34)

The spectrum of the excitation energies of the quasiparticles in a superconductor is given by the poles of both the normal and anomalous thermal Green’s functions Gijp,iωn and Fijp,iωn at T=0. These poles are the denominators that vanish. Ek is the quasiparticle excitation energy given as in Eq. 35. For the generalized two-band model the Ek , which does not include the hopping parameters can be [49]; .

EK=±[12ε12+ε22+¯12+¯22+2M2±ε12+¯12ε22¯222+4M2ε12ε222+¯1+¯2212]12(35)

It is of interest to also evaluate the parameters ¯1,¯2, and ¯M for the case two-band model. They are related to the anomalous Green’s functions by self-consistency conditions Eqs 24.

¯1k=11+12=kkV11k,k`F11p,iωn+kkV12k,k`F22p,iωn(36)
¯2k=22+21=kkV22k,k`F22p,iωn+kkV21k,k`F11p,iωn(37)
¯Mk=M=kkα11k,k`c1kc1k+Q+kkα12k,k`c1kc2k+Q(38)

2.2 Physical properties

2.2.1 Superconductor order parameter

The gap parameter ¯ is a superconducting order parameter, which can be determined self-consistently from the gap equations Eq 36 and Eq 37. In matrix form, the order parameter for the superconducting state is given by;

¯i=jVijFijp,iωn¯j(39)

Where Vij is the pairing interaction constants and functions Fijp,iωn is anomalous green functions in a superconducting state are defined as above and Dωn=ωn2+ε12+¯12ωn2+ε22+¯22, ωn=2n+1π/β and ε12+¯12=E12 and ε22+¯22=E22 energy excitation for the band 1 and 2. Upon substitution, the equation of motion becomes

c1kc1k=ωn2+E22ωn2+E12ωn2+E22(40)
c2kc2k=ωn2+E12ωn2+E12ωn2+E22(41)

Let k,n(12n+12π2+x2)=tanhx22x where x=βE1, by substituting the result

F11p,iωn=β2tanhβE122βE1(42)

And

F22p,iωn=β2tanhβE222βE2(43)

By substituting these two equations, Eq. 42 and Eq. 43 in Eq. 36 and Eq. 37, respectively, the two gape equations in the superconducting state become;

¯1=11+12=1βkV11β2tanhβE122βE1¯1+1βkV12β2tanhβE222βE2¯2(44)
¯2=22+21=1βkV22β2tanhβE222βE2¯2+1βkV21β2tanhβE122βE1¯1(45)

Where V11 and V22 are pairing interactions for 1 and 2 bands, respectively, while the pair interchange between the two bands is assured by the V12 term. The quantity V12 has been supposed to be operative and constant in the energy interval for the higher band and the lower band, keeping in mind the integration range, the gap order parameter satisfies the system. If the intraband interactions are missing, i.e., V11=V22=0, the interband interaction solely induces the transition is V12=V21. Therefore, Eq. 44 and Eq. 45 Become

¯1=12=kV12tanhβE222E2¯2(46)
¯2=21=kV21tanhβE122E1¯1(47)

Converting the summation over k values into an integral with the cut-off energy from ±ћωD, ωD the boson cut-off frequency measured from the Fermi level and introducing the density of state at the Fermi level N10 and N20, and by applying identity 1=2N10 and 2=2N20, then Eq. 46 and Eq. 47 become;

¯1=12=2N10V12ћωDћωD¯2tanhβE222E2dϵ2(48)
¯2=21=2N20V21ћωDћωD¯1tanhβE122E1dϵ1(49)

Now, let’s study these two equations by considering different cases.

Case 1; When T0β which implies, tanhβE2/21andtanhβE1/21, then integral becomes

¯1=12=N10V120ћωD¯2ε22+¯2212dε2(50)
¯2=21=N20V210ћωD¯1ε12+¯1212dε1(51)

Equation 50 and Eq. 51 can be

¯1=12=N10V12¯2ln2ћωD¯2(52)
¯2=21=N20V21¯1ln2ћωD¯1(53)

Now substitution Eq. 53 into Eq. 52

¯1=12=N10V12N20V21¯1ln2ћωD¯1ln2ћωD¯2(54)

For the integral 0ћωD1ε12+¯1212dε1=ln2ћωD¯1 and 0ћωD1ε22+¯2212dε2=ln2ћωD¯2 ,

If V12=V21¯1=12=21=¯2=¯. From this, Eq. 54 becomes,

¯=N10N20V212¯ln2ћωD¯2(55)

By rearranging, 1/V21N10N20=ln2ћωD/¯ taking the exponent

¯=2ћωDexp1V12N10N20(56)

From the BCS theory we have,

3.5KBTC=2¯

Therefore, Eq. 56 for the interband coupling constant of λ12=V12N10N20 becomes;

KBTC=1.14ћωDexp1λ12(57)

This expression for the superconducting transition temperature Tc is like the well-known BCS [7].

Case 2: To obtain a temperature-dependent superconductivity energy gap, we use the expression from Eq. 48 at T=Tc and ¯=0 it gives

1V12N10N20=1λ12=0ћωD1ε2+¯2tanhβε2+¯22dϵ(58)

By using the same techniques as above, this equation becomes

1λ12=ln1.14ћωDKBT¯πKBT21.05(59)

From Eq. 57

1λ12=ln1.14ћωDKBTc(60)

By substituting Eq. 60 into Eq. 59 and by rearranging

lnTTc=¯πKBT21.05(61)

Using the relation ln1x=1xx22+

¯T=3.06KBTc1TTc12(62)

Equation 62 demonstrates how the superconducting order parameter ¯T varies with temperature when the magnetic order parameter is zero and is analogous to the BCS model. From Eq. 62 at T = 0, the superconducting order parameter ¯0=3.06KBTc using the experimental values Tc=20K for Ca0.741La0.261Fe1xCoxAs2 FeSC [30], the superconducting order parameter is ¯0=3.02mev.

2.2.2 The effect of magnetism on the transition temperature (Tc)

From Eq. 36 and Eq. 37, the interchange interaction between the two bands in the two-band model is assured by V12 term. The interband interaction can induce the transition temperature (TC) in both antiferromagnetism and superconducting states. Kristoffele et al. [43] have shown that interband pairing is very efficient in enhancing TC. This is the characteristic feature of the band model. Therefore, the band equations can be

¯1=12=kk`V12(k,k`)F22p,iωn(63)
¯2=21=kk`V21(k,k`)F11p,iωn(64)

From this, to study the band gap, one can write these two equations in simultaneous equations as

¯1+¯2=12+21=2¯=V12kF22p,iωn+kF11p,iωn(65)

Where 12=21, 12+21=2¯ for V12=V21, using the two Green’s functions from equations Eq. 3032 and with energy in Eq. 34 after some mathematical steps, By applying the identity 1=2N10 and 2=2N20 and simplify Eq. 65 gives;

1λ12=0ћωD1ϵ2+¯M2tanhβϵ2+¯M22dϵ0ћωDM¯ϵ2+¯M2tanhβϵ2+¯M22dϵ(66)

Where λ12=V12N1+N2couplingconstant.

At T=TC,¯=0 using Laplace’s Transform, ωωn Matsubara frequency and proceeding through all the necessary steps, the first integral of Eq. 66 becomes

I1=0ћωD1ϵ2+M2tanhβϵ2+M22dϵ=0ћωD2βn=1ωn2+ϵ2+M2dϵ(67)

Where 2n+12π2/β2=ωn2, the Matsubara frequency Eq. 67 becomes

I1=0ћωD2βn=12n+12πβ22+ϵ2+M2dϵ=0ћωD1ϵtanhβϵ2dϵ0ћωDM24βn=01a41+x22dϵ
I1=ln1.14ћωDKBTCMπKBTC21.05(68)

The second integral in Eq. 66, I2 also can be evaluated as

I2=0ћωDM¯ϵ2+¯M2tanhβϵ2+¯M22dϵ(69)

Applying the L’ HOPITAL rule for the ¯0 Eq. 69 can be written as

I2=0ћωDLim¯0M¯ϵ2+¯M2tanhβϵ2+¯M22dϵ(70)

This equation also gives the limit.

I2=Mβ20ћωD1ϵ2+M2sech2βϵ2+M22dϵ(71)

Using sech2x=1tanh2x

0ћωD1ϵ2+M2dϵ0ћωD1ϵ2+M2tanh2βϵ2+M22dϵ(72)

By applying part techniques of integration and using Laplace’s transformation, the Matsubara frequency Eq. 72 gives

12lnM+ћωDMћωD(73)

Substituting this result in Eq. 71 and then Eq. 66 for the value of the two integration results (I1andI2) become

1λ12=ln1.14ћωDKBTCMπKBTC21.05Mβ4lnM+ћωDMћωD(74)

Neglecting the M2 from the expression because it is very small, the transition temperature TC is

TC=1.14ћωDKBe1λ12+aM(75)

Where a=β/4lnM+ћωD/MћωD

2.2.3 Antiferromagnetic transition temperature

We solve a mean-field equation self-consistently with an order parameter in the antiferromagnetism (AFM) state with ordering vector Q, as defined above in Eq. 38 with M as the antiferromagnetic order parameter, α11, and α12k,k` intraband and interband coupling constant of the antiferromagnetic order parameter, respectively, and in the first Brillouin zone of the paramagnetic phase, we take the order vector Q=π,0. Note that to study Co-occurrence, we neglect the inter-orbital order parameter kk`α12k,k`c1kc2k+Q12. Because these values are almost zero for the present orbital models [34]. The anomalous green’s function for the mean-field equation of the intraband pairing interaction from the inverse matrices;

c1kc1k+Q=1βωnеiωnτFijp,iωn(76)

The green function for the antiferromagnetic order parameter

F11p,iωn=c1kc1k+Q=1Dωn¯1ωn2+ε22+¯22(77)

Where Dωn=ωn2+ε12+¯12ωn2+ε22+¯22. Substituting Green’s function in the mean-field equation, the antiferromagnetic order parameter as a function of the superconducting order parameter becomes;

M=kk`α11(k,k`)c1kc1k+Q=α11(k,k`)βωn¯1ωn2+ε12+¯12(78)

Changing the summation into integration and introducing the density of states, N (0), we get

M=α11(kk`)2N10β0ћωD¯1ωn2+ε12+¯12dϵ(79)

Using the Matsubara frequency 2n+12π2/β2=ωn2 and 1/2n+12π2+x2=tanhx/2/2x where x=βE1 Eq. 79 become

M=λATM0ћωD¯1ε12+¯12tanhβε12+¯122dϵ(80)

Where λATM=α11k,k`N1 (0) coupling constant, and E12=ε12+¯12.

Now, let us first solve the integral.

0ћωD1ε12+¯12tanhβε12+¯122dϵ=0ћωD2βn=1ωn2+ε12+¯12dϵ(81)

Using the Laplace transformation and Matsubara frequency, Eq. 81 become

0ћωD2βn=1ωn2+ε12+¯12dϵ=0ћωD2βn=12n+12πβ22+ϵ2+¯12dϵ=0ћωD1ϵtanhβϵ2dϵ0ћωD¯124βn=01a41+x22dϵ(82)

The two integrations give

I1=0ћωD1ϵtanhβϵ2dϵ=ln1.14ћωDKBTATM
I2=0ћωD¯124βn=01a41+x22dϵ¯1πKBTATM21.05

Substituting the I1 and I2 in Eq. 82 and Eq. 80 gives

M=λATM¯ln1.14ћωDKBTATM+¯πKBTATM21.05(83)

Neglecting ¯2 because very small for ¯=¯1 the antiferromagnetic order parameter is

M=λATM¯ln1.14ћωDKBTATM(84)

This equation also gives the antiferromagnetic transition temperature (TATM)

TATM=1.14ћωDKBeMλATM¯0(85)

3 Result and discussion

In this section, we have examined the results obtained from analyzing the normal and anomalous one-particle thermal Green’s functions in a two-band model of superconductivity. The analysis takes into account the potential intraband and interband superconducting interaction terms, which decouple both bands in the mean field approximation. We have derived expressions for the superconducting order parameters (¯1 and ¯2) for the two bands, the superconducting transition temperature (TC) in the pure superconducting state, and the superconducting transition temperature with a magnetic order parameter. Additionally, we have obtained the antiferromagnetic transition temperature (TATM).

Using Eq. 57, we have calculated the theoretical value of TC to be 20.35K, while the experimental value is 20 K [30] for the iron-based superconductor Ca0.741La0.261Fe1xCoxAs2. Additionally, we have employed reasonable approximations for various parameters. For instance, we set the interband coupling constant (λ12) to 0.25 and assumed a substantial boson energy of 84meV within the multi-band model. The electron-phonon coupling constant in FeSC is estimated to be in the range of λ∼0.17–0.21 [50]. Our theoretical prediction is in agreement with the experimental results, providing further evidence that supports their agreement [44, 51]. This observation offers valuable insights into the underlying mechanism responsible for the pairing in superconductivity. Interestingly, it is possible to obtain TC even when all the intraband and interband interactions correspond to repulsion between carriers, as long as the relation λ11λ22λ12λ21>0 is satisfied [52]. This relation is often used for Fe-based superconductors. To visualize the relationship between TC and the interband coupling constant (λ12), we have plotted TC against λ12 in Figure 2A. From the figure, it can be observed that as λ12 increases, the superconducting temperature (TC) for the Ca0.741La0.261Fe1xCoxAs2 superconductor exponentially increases, and vice versa. As the interband coupling constant increases, Tc also increases, supporting the commonly accepted pairing scenario for iron pnictids, which involves SFs mediating pairing.

FIGURE 2
www.frontiersin.org

FIGURE 2. (A) Interband coupling constant (λ12) versus superconducting temperature TC for Ca0.741La0.261Fe1xCoxAs2 Superconductor for ћωD=84mev (B) Phase diagram of superconducting order parameter versus temperature for various values of the magnetic order parameter for Ca0.741La0.261Fe1xCoxAs2 Superconductor by fixing Tc

In the study, the expression for the superconducting order parameter as a function of temperature (Eq. 62) and magnetic ordering was obtained (Eq. 85), and the phase diagram of the superconducting order parameter versus temperature for different values of the magnetic order parameter was plotted (Figure 2B). The results showed that the superconducting order parameter is suppressed when magnetic ordering is present, and this suppression becomes more significant as the value of the magnetic order parameter increases. The impact of magnetic ordering on superconductivity depends on the details of the magnetic structure and the electron bands. The phase diagram of the transition temperature (Tc) versus magnetic ordering was also plotted (Figure 3A) using Eq. 85, and it was observed that magnetic ordering suppresses the superconducting transition temperature. This suppression is likely due to the coupling between localized and conduction electrons, which is strong enough to break up the Cooper pairs. The effect of magnetic ordering on superconductivity depends on the details of the magnetic structure and the electron bands [34].

FIGURE 3
www.frontiersin.org

FIGURE 3. (A) Phase diagram of superconducting transition temperature versus magnetic order Parameter for Ca0.741La0.261Fe1xCoxAs2 Superconductor for ћωD=84mev and λ12=0.25. (B) Phase diagram of antiferromagnetic order temperature versus magnetic order Parameter (Eq. 85) for Ca0.741La0.261Fe1xCoxAs2 Superconductor for λ12=0.25, ћωD=84mev and ¯ (0) = 3.02mev.

The quantum mechanical interaction between the spins of localized electrons and the atomic magnetic moments of the system is the underlying cause of magnetic ordering, suppressing superconductivity. Below the transition temperature (TC), the exchange interaction attempts to align the Cooper pairs, imposing strict limits on the existence of superconductivity. Additionally, we have plotted the phase diagram of TATM (transition temperature to antiferromagnetism) versus M (Figure 3B). It can be observed that magnetic ordering enhances the Neel temperature (TATM). This implies that the antiferromagnetic moment lies in the basal plane for all values of M. By combining Figures 3A, B, we have identified the region between (2.3 < M < 6.07) for M, the magnetic order parameter where superconductivity and antiferromagnetic coexist, as depicted in Figure 4. Our findings are broadly consistent with experimental observations [30]. This also gives additional information for the physics of superconductivity. Based on our discovery, we have found that incorporating the interaction between intra-band and inter-band terms in the two-band model Hamiltonian leads to a multi-band dispersion at the Fermi surface. This is achieved through the utilization of the matrix form of thermal Green function formalism, which holds significance for further comprehensive investigations. By employing mean-field decoupling, it becomes possible to establish the order parameters for superconductivity and magnetism. These order parameters, along with the critical temperatures associated with superconductivity and magnetism, allow us to construct a phase diagram. In mean-field theory, the absence of fluctuations can lead to an overestimation of phase stability, allowing phases that would normally exist in competition to coexist. It is crucial to recognize that mean-field theory has its limits and might not precisely capture the phase boundaries or quantitative characteristics. Still, by encapsulating the fundamental qualitative features of the phase diagram, mean-field theory can offer insightful information. In short, superconducting (SC) and antiferromagnetism (AFM) phases coexist within a certain range of increasing magnetic order parameter values. At lower temperatures, the width of the coexistence region expands as the magnetic order parameter increases, indicating that the coexistence of these phases is influenced by the ellipticity of the electron bands.

FIGURE 4
www.frontiersin.org

FIGURE 4. Phase diagram of SC transition temperature and antiferromagnetic transition temperature versus magnetic order parameter. The figure demonstrates the Co-occurrence of superconductivity and antiferromagnetism (SC + AFM) region between 2.3<M<6.07 for M the magnetic order parameter.

4 Conclusion

In the present work, the possible interplay of superconductivity and antiferromagnetic in Ca0.741La0.261Fe1xCoxAs2 superconductor was studied. Based on the electronic band structure of the iron-based superconductor, we have developed the two-band model Hamiltonian for the system, which consists of intra and interband pairing interactions, and inter-orbital pair hopping. Using the normal and anomalous thermal Green function technique with the two-band model Hamiltonian, the self-consistent gap equations and the expressions for the transition temperatures and order parameters have been obtained. With this mathematical expression and relevant parameters numerically solved the results have been presented in the figure. Figure 2A shows that the interband coupling constant increases as the transition temperature increases for the Ca0.741La0.261Fe1xCoxAs2 iron base superconductor. The superconducting order parameter (¯) gets to zero at transition temperature (Tc), and it is suppressed when magnetic order parameter (M) sets in and the suppression becomes more important when the antiferromagnetic correlations grow. This is verified by plotting the phase diagrams of the superconducting order parameter (¯) versus temperature (T) by varying the value of M. Furthermore, Figure 3A, which is plotted Tc as the function of the magnetic order parameter of Ca0.741La0.261Fe1xCoxAs2 system indicates that Tc is decreasing when M is increased. On the other hand, in the triclinic state, magnetic ordering enhances the antiferromagnetism transition temperature (TAFM) as indicated in the phase diagram of antiferromagnetism order temperature (TAFM) versus antiferromagnetic order parameter M which is shown in Figure 3B. Lastly, by merging Figures 3A, B, we have found the intersection region of superconductivity and antiferromagnetic. The region under the two merged Figures shows the Co-occurrence of the two states established in the magnetic order parameter range of 2.3<M<6.07 for the system Ca0.741La0.261Fe1xCoxAs2 as shown in Figure 4 which is seen to the broad experimental agreement.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

GS: Conceptualization, Data curation, Investigation, Methodology, Software, Supervision, Writing–original draft, Writing–review and editing. DS: Conceptualization, Investigation, Methodology, Supervision, Writing–review and editing.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Kamerlingh-Onnes H. Disappearance of the electrical resistance of mercury at helium temperatures. Proc Koninklijke Akademie van Wetenschappen te Amsterdam (1911) 14:113.

Google Scholar

2. Van Delft D, Kes P. The discovery of superconductivity. Phys Today (2010) 63(9):38–43. doi:10.1063/1.3490499

CrossRef Full Text | Google Scholar

3. Gavaler JR. Superconductivity in Nb–Ge films above 22 K. Appl Phys Lett (1973) 23(8):480–2. doi:10.1063/1.1654966

CrossRef Full Text | Google Scholar

4. Kasap S, Capper P. Springer handbook of electronic and photonic materials. Berlin, Germany: Springer (2017).

Google Scholar

5. Meissner W, Ochsenfeld R. Ein neuer effekt bei eintritt der supraleitfähigkeit. Naturwissenschaften (1933) 21(44):787–8.

CrossRef Full Text | Google Scholar

6. Buckel W, Kleiner R. Superconductivity: fundamentals and applications. Hoboken, New Jersey, United States: John Wiley and Sons (2008).

Google Scholar

7. Bardeen J, Cooper LN, Schrieffer JR. Theory of superconductivity. Phys Rev (1957) 108(5):1175–204. doi:10.1103/physrev.108.1175

CrossRef Full Text | Google Scholar

8. Ott HR, Rudigier H, Fisk Z, Smith JL. U be 13: an unconventional actinide superconductor. Phys Rev Lett (1983) 50(20):1595–8. doi:10.1103/physrevlett.50.1595

CrossRef Full Text | Google Scholar

9. Sharma RG. A review of theories of superconductivity. Supercond Basics Appl Magnets (2021) 123–60.

CrossRef Full Text | Google Scholar

10. Hasan MS, Ali SS. Properties and types of superconductors. Mater Res Found (2022) 132:17–48.

CrossRef Full Text | Google Scholar

11. Kamihara Y, Watanabe T, Hirano M, Hosono H. Iron-based layered superconductor La [O1-x F x] FeAs (x= 0.05− 0.12) with T c= 26 K. J Am Chem Soc (2008) 130(11):3296–7. doi:10.1021/ja800073m

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Schafgans AA, Moon SJ, Pursley BC, LaForge AD, Qazilbash MM, Sefat AS, et al. Electronic correlations and unconventional spectral weight transfer in the high-temperature pnictide BaFe 2− x Co x as 2 superconductor using infrared spectroscopy. Phys Rev Lett (2012) 108(14):147002. doi:10.1103/physrevlett.108.147002

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Cooper LN. Origins of the theory of superconductivity. IEEE Trans Magn (1987) 23(2):376–9. doi:10.1109/tmag.1987.1065151

CrossRef Full Text | Google Scholar

14. Stewart GR. Superconductivity in iron compounds. Rev Mod Phys (2011) 83(4):1589–652. doi:10.1103/revmodphys.83.1589

CrossRef Full Text | Google Scholar

15. Kordyuk AA. Iron-based superconductors: magnetism, superconductivity, and electronic structure. Low Temp Phys (2012) 38(9):888–99.

CrossRef Full Text | Google Scholar

16. Lee C-H, et al. Effect of structural parameters on superconductivity in fluorine-free LnFeAsO1-y (Ln= La, Nd). J Phys Soc Jpn (2008) 77(8):83704.

CrossRef Full Text | Google Scholar

17. Nandi S, Kim MG, Kreyssig A, Fernandes RM, Pratt DK, Thaler A, et al. Anomalous suppression of the orthorhombic lattice distortion in superconducting Ba (Fe 1− x Co x) 2 as 2 single crystals. Phys Rev Lett (2010) 104(5):057006. doi:10.1103/physrevlett.104.057006

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tanui PK, Namwetako JS, Ingosy A, Khanna KM. Superconductivity and electron-hole superconductivity. Int J Recent Res Asp (2020) 7(3).

Google Scholar

19. Ren Z-A, Che GC, Dong XL, Yang J, Lu W, Yi W, et al. Superconductivity and phase diagram in iron-based arsenic-oxides ReFeAsO1− δ (Re= rare-earth metal) without fluorine doping. Europhys Lett (2008) 83(1):17002. doi:10.1209/0295-5075/83/17002

CrossRef Full Text | Google Scholar

20. Matsuishi S, Inoue Y, Nomura T, Yanagi H, Hirano M, Hosono H. Superconductivity induced by Co-doping in quaternary fluoroarsenide CaFeAsF. J Am Chem Soc (2008) 130(44):14428–9. doi:10.1021/ja806357j

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Rotter M, Tegel M, Johrendt D. Superconductivity at 38 K in the iron arsenide (Ba 1− x K x) Fe 2 as 2. Phys Rev Lett (2008) 101(10):107006. doi:10.1103/physrevlett.101.107006

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Kudo K, Iba K, Takasuga M, Kitahama Y, Matsumura J, Danura M, et al. Emergence of superconductivity at 45 K by lanthanum and phosphorus co-doping of CaFe2As2. Sci Rep (2013) 3(1):1478. doi:10.1038/srep01478

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Pitcher MJ, Parker DR, Adamson P, Herkelrath SJC, Boothroyd AT, Ibberson RM, et al. Structure and superconductivity of LiFeAs. Chem Commun (2008) 45:5918–20. doi:10.1039/b813153h

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Tapp JH, Tang Z, Lv B, Sasmal K, Lorenz B, Chu PCW, et al. LiFeAs: an intrinsic FeAs-based superconductor with T c= 18 K. Phys Rev B (2008) 78(6):060505. doi:10.1103/physrevb.78.060505

CrossRef Full Text | Google Scholar

25. Hsu F-C, Luo JY, Yeh KW, Chen TK, Huang TW, Wu PM, et al. Superconductivity in the PbO-type structure α-FeSe. Proc Natl Acad Sci (2008) 105(38):14262–4. doi:10.1073/pnas.0807325105

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Mizuguchi Y, Takano Y. Review of Fe chalcogenides as the simplest Fe-based superconductor. J Phys Soc Jpn (2010) 79(10):102001. doi:10.1143/jpsj.79.102001

CrossRef Full Text | Google Scholar

27. Löhnert C, Stürzer T, Tegel M, Frankovsky R, Friederichs G, Johrendt D. Superconductivity up to 35 K in the iron platinum arsenides (CaFe1− xPtxAs) 10Pt4− yAs8 with layered structures. Angew Chem Int. Ed. (2011) 50(39):9195–9. doi:10.1002/anie.201104436

CrossRef Full Text | Google Scholar

28. Katayama N, Kudo K, Onari S, Mizukami T, Sugawara K, Sugiyama Y, et al. Superconductivity in Ca1-<i>x</i>La<i>x</i>FeAs2: a novel 112-type iron pnictide with arsenic zigzag bonds. J Phys Soc Japan (2013) 82(12):123702. doi:10.7566/jpsj.82.123702

CrossRef Full Text | Google Scholar

29. Yakita H, Ogino H, Okada T, Yamamoto A, Kishio K, Tohei T, et al. A new layered iron arsenide superconductor:(Ca, Pr) FeAs2. J Am Chem Soc (2014) 136(3):846–9. doi:10.1021/ja410845b

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Jiang S, et al. Coexistence of superconductivity and antiferromagnetism in Ca $ _ {0.74 (1)} $ La $ _ {0.26 (1)} $(Fe $ _ {1-x} $ Co $ _ {x} $) as $ _ {2} $ single crystals (2016). https://arxiv.org/abs/1603.04899.

Google Scholar

31. Vavilov MG, V Chubukov A, Vorontsov AB. Jump in specific heat in the presence of a spin-density wave at the superconducting transition in iron pnictides. Phys Rev B (2011) 84(14):140502. doi:10.1103/physrevb.84.140502

CrossRef Full Text | Google Scholar

32. Jiang S, et al. Structural and magnetic phase transitions in Ca 0.73 La 0.27 FeAs 2 with electron-overdoped FeAs layers. Phys Rev B (2016) 93(5):54522.

CrossRef Full Text | Google Scholar

33. Vorontsov AB, Vavilov MG, V Chubukov A. Superconductivity and spin-density waves in multiband metals. Phys Rev B (2010) 81(17):174538. doi:10.1103/physrevb.81.174538

CrossRef Full Text | Google Scholar

34. Matsui Y, Morinari T, Tohyama T. Coexistence of antiferromagnetism and superconductivity in iron-based superconductors. J Phys Soc Jpn (2014) 83(9):094703. doi:10.7566/jpsj.83.094703

CrossRef Full Text | Google Scholar

35. Raghu S, Qi X-L, Liu C-X, Scalapino DJ, Zhang S-C. Minimal two-band model of the superconducting iron oxypnictides. Phys Rev B (2008) 77(22):220503. doi:10.1103/physrevb.77.220503

CrossRef Full Text | Google Scholar

36. Parker D, Vavilov MG, V Chubukov A, Mazin II. Coexistence of superconductivity and a spin-density wave in pnictide superconductors: gap symmetry and nodal lines. Phys Rev B (2009) 80(10):100508. doi:10.1103/physrevb.80.100508

CrossRef Full Text | Google Scholar

37. Lorenzana J, Seibold G, Ortix C, Grilli M. Competing orders in FeAs layers. Phys Rev Lett (2008) 101(18):186402. doi:10.1103/physrevlett.101.186402

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Izyumov Y, Kurmaev E. High-Tc superconductors based on FeAs compounds. Berlin, Germany: Springer Science and Business Media (2010).

Google Scholar

39. Thomale R, Platt C, Hu J, Honerkamp C, Bernevig BA. Functional renormalization-group study of the doping dependence of pairing symmetry in the iron pnictide superconductors. Phys Rev B (2009) 80(18):180505. doi:10.1103/physrevb.80.180505

CrossRef Full Text | Google Scholar

40. V Chubukov A, Vavilov MG, Vorontsov AB. Momentum dependence and nodes of the superconducting gap in the iron pnictides. Phys Rev B (2009) 80(14):140515. doi:10.1103/physrevb.80.140515

CrossRef Full Text | Google Scholar

41. Maier TA, Graser S, Scalapino DJ, Hirschfeld PJ. Origin of gap anisotropy in spin fluctuation models of the iron pnictides. Phys Rev B (2009) 79(22):224510. doi:10.1103/physrevb.79.224510

CrossRef Full Text | Google Scholar

42. V Chubukov A, V Efremov D, Eremin I. Magnetism, superconductivity, and pairing symmetry in iron-based superconductors. Phys Rev B (2008) 78(13):134512. doi:10.1103/physrevb.78.134512

CrossRef Full Text | Google Scholar

43. Nuwal A, Kakani S, Kakani SL. Two band model for the iron based superconductors. Indian J Pure Appl Phys (2015) 52(6):411–22.

Google Scholar

44. Wen XG. Quantum field theory of many-body systems: from the origin of sound to an origin of Light and electrons. Quan F. Theor Many-Body Syst. Orig. Sound Orig. Light Electrons (2010) 1–520. doi:10.1093/acprof:oso/9780199227259.001.0001

CrossRef Full Text | Google Scholar

45. Inkson JC. Many-body theory of solids: an introduction. Berlin, Germany: Springer Science and Business Media (2012).

Google Scholar

46. Nambu Y. Quasi-particles and gauge invariance in the theory of superconductivity. Phys Rev (1960) 117(3):648–63. doi:10.1103/physrev.117.648

CrossRef Full Text | Google Scholar

47. Rickayzen G. Green’s functions and condensed matter. Massachusetts, United States: Courier Corporation (2013).

Google Scholar

48. Mahan GD. Many-particle physics. Berlin, Germany: Springer Science and Business Media (2000).

Google Scholar

49. Kristoffel N, Örd T. On the possibility of HTSC in a two-band model with a semiconducting gap (A3C60). Phys Status Solidi (1993) 175(1):K9–K12. doi:10.1002/pssb.2221750129

CrossRef Full Text | Google Scholar

50. Brydon PMR, Daghofer M, Timm C. Magnetic order in orbital models of the iron pnictides. J Phys Condens Matter (2011) 23(24):246001. doi:10.1088/0953-8984/23/24/246001

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Fernandes RM, Pratt DK, Tian W, Zarestky J, Kreyssig A, Nandi S, et al. Unconventional pairing in the iron arsenide superconductors. Phys Rev B (2010) 81(14):140501. doi:10.1103/physrevb.81.140501

CrossRef Full Text | Google Scholar

52. Zhou W, Ke F, Xu X, Sankar R, Xing X, Xu CQ, et al. Correlation between non-Fermi-liquid behavior and superconductivity in (Ca, La)(Fe, Co) as 2 iron arsenides: a high-pressure study. Phys Rev B (2017) 96(18):184503. doi:10.1103/physrevb.96.184503

CrossRef Full Text | Google Scholar

Keywords: Co-occurrence, iron-based superconductors, green’s function, order parameter, intra-band and inter-band

Citation: Sherka GT and Shiferaw DA (2024) Theoretical investigation of the Co-occurrence of superconductivity and antiferromagnetism in iron-based high-temperature superconductors. Front. Phys. 12:1356768. doi: 10.3389/fphy.2024.1356768

Received: 16 December 2023; Accepted: 16 February 2024;
Published: 14 March 2024.

Edited by:

Jian-Xin Zhu, Los Alamos National Laboratory (DOE), United States

Reviewed by:

Christopher Lane, Los Alamos National Laboratory (DOE), United States
Jianbao Zhao, Canadian Light Source, Canada

Copyright © 2024 Sherka and Shiferaw. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Gebru Tesfaye Sherka, gebrut70k@gmail.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.