Abstract

The main focus of this article is to investigate the theoretical interplay of magnetism and superconductivity in a two-band model for the iron-based superconductor Ba1−xKxFe2As2. On the basis of experimental results, the two-band model Hamiltonian was considered. We obtained mathematical statements for the superconductor Ba1−xKxFe2As2 superconducting (SC) transition temperature, spin-density-wave (SDW), transition temperature, superconductivity order parameter, and SDW order parameter from the Bogoliubov transformation formalism and the model Hamiltonian. Furthermore, an expression for the dependence of the SDW transition temperature on the SDW order parameter and the dependence of the SC transition temperature on the SDW order parameter was obtained for Ba1−xKxFe2As2. By substituting the experimental and theoretical values of the parameters in the derived statements, phase diagrams of the SC transition temperature versus the SDW order parameter and the SDW transition temperature versus the SDW order parameters have been plotted to demonstrate the dependence of the SDW order parameter on transition temperatures. By combining the two-phase diagrams, we depicted the possible coexistence of superconductivity and magnetism for the Ba1−xKxFe2As2 superconductor. The phase diagrams of temperature versus SC order parameter and temperature versus SDW order parameter were also plotted to show the dependence of order parameters on temperature for the Ba1−xKxFe2As2 superconductor.

1. Introduction

Fe-based superconductors (FebSc), which greatly build up the family of unconventional superconductors, have many different systems. The major FebSc systems are 1111, 122, 111, 11, 245, 42622, 1144, and 12442 (e.g., LaFeAsO, Ba1−xKxFe2As2, LiFeAs, FeSe1−xTex, Rb2Fe4Se5, Sr4V2O6Fe2As2, RbEuFe4As4, and KCa2Fe4As4F2) [15]. Iron-based square-planar sheets, which are essential for the unconventional superconductivity of all FebSc systems, represent their shared structural component [6]. The parent compounds of FebSc systems are antiferromagnets, spin-density-wave (SDW) metals, and they are magnetic bad metals. Increasing doping in 122 systems can destroy antiferromagnetism and lead to superconductivity [7]. Superconductivity occurs in specific planes, Fe-As plane [8]. The iron-containing plane of the superconductor, in which pnictogen or chalcogen atoms protrude above and below the plane, is not flat because the pnictogen and chalcogen atoms are much larger than iron atoms [9]. The pnictogen and chalcogen atoms are tightly grouped in a tetrahedral structure [10]. Due to the close packing of the Fe atoms, the five Fe 3d orbitals serve as charge carriers in superconductivity [11]. The sign-changing s-wave symmetry characterizes the FebSc symmetry [12].

Most of the iron-pnictides become superconductors only when doped (added an impurity) with holes or electrons. depends on doping concentration [13]. For the 122 iron-based superconductors family, AFe2An2, where (A indicates alkaline earth metals such as Ba, Sr or Ca and Eu) and An is a pnictide (As, P), superconductivity can be achieved with the introduction of dopants [14]. There are several ways to introduce dopants [15]. These are (1) hole doping is achieved with substituting A for monovalent B+ (B Cs, K, Na) atoms partially in the blocking layer, and this substitution should add an excess hole into the system, for example, Ba1−xKxFe2As2 [16, 17] (2) partially substitute Fe for transition metals (Co, Ni, Pd, Rh) into FeAs layers and yields electrons into the system. In this method, dopants are directly doped into the Fe layer, which can stabilize the system; for example, Co (A(Fe1−xCox)2As2), Rh (A(Fe2−xRhx)As2), Ni (A(Fe1−xNix)2As2) [1820] and we get electron-doped pnictides that form an abundant phase diagram where superconductivity and magnetism exist together, and (3) replacing arsenic partially with phosphorus, and phosphorus generates a chemical pressure effect that suppress SDW and emerges superconductivity at the corresponding unit cell volume [21, 22].

The parent compound BaFe2As2 is an iron arsenide compound and crystallizes with the tetragonal ThCr2Si2 type structure, space group I4/mmm [23]. Its crystal structure and SDW transition are strongly affected by hole doping with potassium [24]. The compound BaFe2As2 changes structural and magnetic transitions simultaneously at [25]. Magnetic transition is strongly coupled with structural distortion. Potassium substitutions with hole doping suppress the SDW state of BaFe2As2, resulting in the befall of superconductivity [26]. The suppression of long-range magnetic ordering and the emergence of the superconducting (SC) state at the same time are correlated, suggesting that the spin fluctuation of the Fe moments is crucial in creating the superconductive state [27]. The coupling between lattice and spins results from the proximity of the structural and magnetic transitions [28]. The crystal distortion is driven by magnetisms because lower symmetry allows the spins to order and thus relieve magnetic frustrations [6]. Superconductivity, which is closely related to electron pair instability on the Fermi surface, SDW, which has been deduced from the nesting Fermi surface and observed directly in neutron diffraction experiments in the BaFe2As2 (122) class of materials. The same electronic state would participate in both the SC and SDW electron-pairs instability of the Fermi surface, and indeed the hole and electron-like Fermi sheets in Ba1−xKxFe2As2 have been observed [29].

Experiment has been revealed that in potassium substitution barium iron arsenide superconductor (Ba1−xKxFe2As2), superconductivity appears for in both the orthorhombic and tetragonal phases and for , the SDW and SC orders interplay at low temperature [30]. The coexistence of magnetism and superconductivity strongly emphasizes that magnetic spin fluctuations are involved in the Cooper pairing in Fe-based superconductors [24]. As the doping level (x in Ba1−xKxFe2As2) increases, the electron Fermi pocket continuously shrinks, whereas the hole pocket becomes larger. SDW ordering coexisting with superconductivity at . When the doping concentration reaches 0.3, the SDW transition is completely suppressed, and increases to and at , . The electron pocket finally vanishes beyond , where superconductivity also disappears. In this article, we study the interaction (coexistence) between magnetism and superconductivity in a two-band model for the iron-based superconductor Ba1−xKxFe2As2 by using the Bogoliubov transformation formalism.

Based on the concepts of electronic structure, this article theoretically investigates the coexistence of magnetism and superconductivity with potassium substitution (or hole doping) in the blocking layer of barium iron arsenide (Ba1−xKxFe2As2) superconductor in the two-band model. By considering a two-band model of Hamiltonian and using the Bogoliubov transformation formalism methodology, It tried to find the mathematical expression for the SC critical temperature , SC order parameter , the SDW order parameter and SDW transition temperature ().

2. Mathematical Formulation of the Problem

A system of Hamiltonian conduction electron interacting with phonons for the interplay of superconductivity and magnetism in our compound (Ba1−xKxFe2As2) in a two-band model, band and band , can be expressed as follows [3135]:

Here, the first term is the kinetic energy. The second term is the kinetic energy of holes. The third term is the pair interaction of the SC order parameter, and the fourth term is the pair hopping between the fermion SC channel due to magnetic interactions with some coupling constant (). and are the wave vector and the spin projection on the z-axis. is the interaction potential of fermions, and energies of free electrons measured with respect to the chemical potential. and are creation (annihilation) operators of the electrons and holes, respectively.

To decouple the model Hamiltonian, one uses the mean-field approximation, and thus the mean-field Hamiltonian of both SC state and SDW state, expressed as follows:where the mean fields (order parameters) are as follows:where is the SC interaction potential.where is the magnetic interaction potential.

2.1. Pure SC State

A Bogoliubov transformation can be used to drive the gap equations by diagonalizing the mean-field Hamiltonian. That means a simple 2D rotation can diagonalize the mean-field Hamiltonian. Such a rotation is equivalent to introducing Bogoliubov operators, where the fermionic anticommutation relations remain valid. This kind of rotation is known as canonical or Bogoliubov transformation. The mean-field Hamiltonian in pure SC state Equation (2) can be written as follows:

The mean-field Hamiltonian in Equation (6) can be diagonalized using the relevant Bogoliubov transformation or quasi-particle transformation, and its diagonalized form is as follows:

For this transformation, the Bogoliubov quasiparticle matrix is defined as follows:

Their relation to fermion operators is defined as follows:where , , , and are coherence factors.

The diagonalized Hamiltonian has the form as follows:

The electron operators and can be expressed in terms of the quasiparticle operators as follows:where f is the Fermi–Dirac statistics, is the Energy at Fermi level.

Applying the properties of Fermi–Dirac statistics, the term becomes the following:

The coherence factors related as follows:

Inserting Equations (13) and (14) into Equation (12), one can get the following:

Substituting Equation (15) into Equation (3), the gap parameter becomes the following:where . Equation (16) is the SC order parameter in the s band.

Similarly, the b-band mean field value is found as follows:where is the Fermi–Dirac statistics, is the energy at the Fermi level.

Applying the properties of Fermi–Dirac statistics, the term also becomes the following:

The coherence factors related as follows:

Inserting Equations (18) and (19) into Equation (17), one gets the following:

Substituting Equation (13) into Equation (4), the gap parameter becomes the following:where . Equation (21) is the SC order parameter in the b band.

2.2. Pure SDW State

The Hamiltonian Equation (2) can yield an SDW phase that is a pure magnetic phase for zero SC interaction. The mean-field Hamiltonian for the pure magnetic phase is as follows:

The mean-field Hamiltonian in Equation (22) can be diagonalized using the relevant Bogoliubov transformation or quasiparticle transformation, and its diagonalized form is as follows:

The Bogoliubov quasi-particle matrix for our transformation is as follows:

Operators in the Hamiltonian are expressed as follows:

The coherency factors are given by the following:

Therefore, the diagonalized Hamiltonian will have the form as follows:where and are the eigenvalues and given by an expression

and

Using the Bogoliubov transformation, we defined the SDW state in Equation (23); we can now calculate the expectation value of the operators in Equation (22) to find the SDW gap equation.

Substituting Equation (28) into Equation (5), it gives the SDW gap equation.

2.3. SC Order Parameter

There is clear information that show the coexistence of SDWs and superconductivity in both hole- and electron-doped families in most FebSc [22]. In our compound, Ba1−xKxFe2As2, Superconductivity can be obtained by (1) chemical doping [3638] and/or (2) the application of pressure [39]. The addition of a chemical dopant and/or the application of pressure on BaFe2As2 results in the eradication of the SDW and, therefore, in the reduction of magnetic order. Superconductivity is observed when the magnetic order is sufficiently weak; thus, superconductivity and weak magnetic order can still coexist [37, 38, 40]. The exact mechanism for the emergence of superconductivity is, however, still unknown. Therefore, we need to study the coupled system below the SC transition temperature. The mean-field Hamiltonian without the interaction part for both the SC state and the magnetic state is stated in Equation (2). Applying the second quasi-particle transformation to Equation (2), it gives the following:

Using Nambu notionwhere the Nambu spinors are as follows:

The energy matrix

Diagonalizing Equation (31) for each band gives the following:where and are the Bogoliubov quasi-particle matrix which is defined as the following transformations:

Their relation to fermion operators is defined as follows:where , , and are coherence factors that satisfy the following relation:andwhere and are the eigenvalues or excitations and are given by the following:

Where and

Therefore, the diagonalized Hamiltonian will have the form as follows:

Following the same procedure as for the SC state and the magnetic state, the gap equations for both the magnetic and SC states, the gap equations or order parameters are given by the following:where and .

Near the edge of the Fermi-surface, there is homogeneity, and the two order parameters can be taken as equal and opposite; with this assumption, we can take , and . With this in mind, the two SC gap equations of Equation (40) are simplified to the following:where and

For simplicity, the interaction potential and Equation (41) is simplified to the following:where is called the effective order parameter.

Changing summation to integration in the region and by introducing the density of state at the Fermi level of the interband interaction, that is as follows:

The density of state is equal to and Equation (42), becomes the following:where

If , and evaluating the integrals,

Let and , then Equation (56) written as follows:

From Laplacian transform with the Matsuber relation result, we can write Equation (58) as follows:

Apply the following equality.where and

From the above equality relations and

For low-temperature and denotes the Euler’s constant and its value is given by . and this means .so after some steps Equation (66) can be written as follows:

The second integral can be calculated using L’Hopitals Rule as follows:

After some steps, Equation (72) becomes the following:

Now Equation (54) yields,

Applying the same procedure can be written as follows:

Substituting Equations (74) and (75) into Equation (51), it gives the following:

For small , one can ignore the term. Thus, Equation (76) reduces to the following:

The expression for dependence of on in the inter band interaction becomes the following:

Equation (81) clearly shows the dependence of the SC transition temperature on the SDW order parameter. In the pure diamagnetism region, , and Equation (81) becomes the following:

This is like the well-known BCS expression for the SC transition temperature . From Equation (82), the SC coupling parameter can be found as follows:where , [21, 25, 30] and . Substituting these values in Equation (82.1), we will have

For perfect diamagnetism Equation (50) gives the following:

After some steps, Equation (83) becomes the following:

From the BCS theory,

Using logarithmic series

Ignoring the high-order terms, it reduces to the following:

2.4. Magnetic Order Parameter

Using the Bogoliubov transformation, we defined for SDW state in Equation (23); we are now able to calculate the expectation value of the operators in Equation (22) with SC state, the SDW gap equation or SDW order parameter is given by the following:

Applying the same procedure as for the pure SDW region, Equation (94) is simplified.where

Equation (95) can be simplified as follows:where , which is the effective order parameter. .

Changing the summation to integration in the region and by introducing the density of state at the Fermi level is that is as follows:

Inserting Equation (99) into Equation (98), we get the following:

Say , which is called the magnetic coupling parameter.

After some steps, one gets,

For small values of , we ignore term. Thus, Equation (107) reduces to the following:where is the SDW coupling parameter. Equation (110) shows that the SDW order parameter increases as the SDW transition temperature increases. For , Equation (110) is simplified to the following:where , [21, 22, 36] and . Substituting these values in Equation (110.1), we will have .

For pure magnetic region and Equation (106) becomes the following:

Let say and

From Laplacian transform with the Matsuber relation result, we can write the equation as follows:

Now let us apply the following equality:where and

From the above equality relations and

For low-temperature and denotes the Euler’s constant and its value is given by . and this means . so Equation (123) can be written as follows:

Therefore, it gives the following:

For a very small value of , goes to zero and . The magnetic coupling parameter becomes as follows:

Substituting Equation (127) into Equation (126), we get the following:

After some steps, we get the following:

The expression for the temperature dependence of SDW in the temperature range of [36] becomes the following:

This equation tells that if the temperature increases, the SDW order parameter decreases.

3. Result and Discussion

In this chapter, the effect of temperature on the SC order parameter and SDW order parameter on the SC transition temperature and on SDW transition temperature. It also describes the analysis we made for the conceptual analysis of the interplay of superconductivity and SDW in a two-band model for Ba1−xKxFe2As2. Using a two-band model Hamiltonian and using Bogoliubov transformation formalism, one obtained the mathematical expressions for the SC transition temperature , the SDW transition temperature, the SC order parameter and the SDW order parameter .

From Equations (82) and (93) and using the MATLAB script, we plot the phase diagram versus as Figure 1(a) illustrates this. Based on Equation (130), the phase diagram showing how the magnetic order parameter () depends on the temperature in the pure magnetic region shown in Figure 1(b) was plotted.

As seen in Figure 1(a), when the temperature increases, the SC order parameter decreases and gets lost as the temperature is equal to the critical temperature. The maximum value of SC order parameter, , occurs at and it vanishes at the SC transition temperature . As Figure 1(b) shows, magnetism decreases as the temperature increases and becomes lost at the SDW transition temperature . The maximum value of the SDW order parameter, , occurs at . This finding is also in agreement with experimental observations [21, 23, 25, 30, 37].

Based on Equation (81), we plotted the phase diagrams of versus . As can be observed in Figure 2(a), the SC transition temperature is suppressed as the value of the the SDW order parameter increases for Ba1−xKxFe2As2. From this figure, one can observe that the SDW order parameter suppresses superconductivity and increases the magnetic nature of the system.

Using Equation (110), one also plotted the phase diagram of the transition temperature versus SDW order parameter (Figure 2(b)). As can be observed from the figure, increases with increasing SDW order parameter of the Ba1−xKxFe2As2 superconductor.

Finally, by merging Figures 2(a) and 2(b), one depicted a region where both superconductivity and SDW coexist, as shown in Figure 3. This figure shows the possible interplay of superconductivity and SDW for Ba1−xKxFe2As2. As indicated in this figure, this finding is in agreement with experimental observations [21, 25, 30]. This figure also depicts there are regions that show SC state and antiferromagnetic state segregates, which indicates that there are regions where magnetic and SC phases are not mixed.

4. Conclusions

In this work, the possible coexistence of superconductivity and magnetism has been studied in the two-band model for the iron-based superconductor Ba1−xKxFe2As2. The superconductivity order parameter is weakened when the temperature increases and gets lost at the SC critical temperature. The SC transition temperature is weakened as the magnitude of the SDW order parameter increases. The SDW transition temperature increases with increasing SDW order parameter for Ba1−xKxFe2As2. This finding also proofs the dependence of the SDW order parameter on the temperature in the pure magnetic region of the iron-based superconductor Ba1−xKxFe2As2. The SDW order parameters is weakened when the temperature increases and vanishes at the SDW transition temperature. We have depicted the possible interplay of superconductivity and magnetism of the Ba1−xKxFe2As2 superconductor. The possible interplay of superconductivity and magnetism in Ba1−xKxFe2As2 initiates further study in the field of condensed matter physics for further studies. The results obtained in this article work are in agreement with previous studies [25, 30] that the coexistence implies that there is a weak exchange coupling between the itinerant SC electrons and the localized ordered spins in the Ba1−xKxFe2As2 superconductor.

Data Availability

The data used to support the findings of this study are included in the article.

Ethical Approval

The conducted research is not related to either human or animal use.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors contributed equally to this study. That is, they analyzed the data, drafted and edited the manuscript, and read and approved the final manuscript.