Emergent electronic structure of CaFe2As2

TL;DR

This study uses density functional theory to analyze how pressure-induced structural changes in CaFe2As2 affect its electronic properties, including Fermi surface evolution and magnetic behavior, shedding light on its exotic phenomena.

Contribution

It provides detailed insights into the electronic structure and stability of CaFe2As2's collapsed tetragonal phase under pressure, highlighting the role of hybridization and Fermi surface changes.

Findings

cT phase is more stable with enhanced hybridization

Fermi surface around Γ point vanishes with pressure

Magnetic moment decreases in cT phase

Abstract

CaFe2As2 exhibits collapsed tetragonal (cT) structure and varied exotic behavior under pressure at low temperatures that led to debate on linking the structural changes to its exceptional electronic properties like superconductivity, magnetism, etc. Here, we investigate the electronic structure of CaFe2As2 forming in different structures employing density functional theory. The results indicate better stability of the cT phase with enhancement in hybridization induced effects and shift of the energy bands towards lower energies. The Fermi surface centered around $\Gamma$ point gradually vanishes with the increase in pressure. Consequently, the nesting between the hole and electron Fermi surfaces associated to the spin density wave state disappears indicating a pathway to achieve the proximity to quantum fluctuations. The magnetic moment at the Fe sites diminishes in the cT phase consistent with the magnetic susceptibility results. Notably, the hybridization of Ca 4s states (Ca-layer may be treated as a charge reservoir layer akin to those in cuprate superconductors) is significantly enhanced in the cT phase revealing its relevance in its interesting electronic properties.

Tables

Table 1. Table 1: The magnetic moment of Fe and ground state energy at various structural configuration of CaFe 2 As 2 , the experimental latice parameters are taken from neutron-kreyssing-PRB08

Phase	$p$	$a$	$c$	$z_{As}$	$d_{FeAs}$	$d_{AsAs}$	${\mu }_{Fe}$	$E\left(Mag\right)$	$E\left(NM\right)$
	(GPa)	(Å)	(Å)		(Å)	(Å)	(${\mu }_B$)	(meV)	(meV)
T	0	3.8915	11.690	0.372	2.4145	2.9926	2.37	0	150
cTp1	0.47	3.9785	10.6178	0.3657	2.3380	2.8519	1.37	60	-198
cTp2	0.63	3.9780	10.6073	0.3663	2.3405	2.8364	1.95	100	-194

Full text

I Introduction

CaFe2As2 exhibits a rich temperature - pressure phase diagram and provides an ideal platform to study the interplay between magnetism and superconductivity in Fe-based compounds, where the superconductivity is believed to be driven by magnetic fluctuations. The crystal structure of CaFe2As2 at room temperature and ambient pressure is tetragonal possessing $I4/mmm$ space group with lattice parameters $a$ = 3.8915(2)Å, $c$ = 11.690(1)Å, and $z_{As}$ = 0.372(1)phasetransition_NiN_PRB08 ; structure_Tompsett_physicab2009 ; neutron-kreyssing-PRB08 . Upon cooling at ambient pressure, it undergoes a concomitant transition to spin density wave (SDW) state and orthorhombic (O) structure from its room temperature paramagnetic tetragonal (T) structure at about 170 K phasetransition_NiN_PRB08 . The crystal structures in ambient conditions are shown in Fig. 1. On application of pressure ($P$), both the structural and magnetic transitions get suppressed, and a collapsed tetragonal (cT) phase emerges with no magnetic long range order neutron-kreyssing-PRB08 ; ct-pratt-goldman-PRB09 ; ct-sc-MKumar-PRB08 ; neutron-Goldman-PRB08 ; sc-ct-torikachi-PRL08 . Superconductivity can be achieved in this material at an intermediate pressure range or via suitable electron or hole doping sc-ct-torikachi-PRL08 ; sc_cafe2as2_saha_PRB_2012 .

Historical Background

Electronic structure calculations based on first principle approaches were employed extensively to address the ground state properties, Fermi surface and pressure induced effects structure_Tompsett_physicab2009 ; dft_yildrim_prl2009 ; ganesh-EuFe2As2 ; dft_Colonna_prb_2011 ; dft_Zhang_prb_2009 . While ligand $p$ states play major role in the electronic properties of cuprate superconductors (the most studied unconventional superconductors), the main contribution at the Fermi level, $\epsilon_{F}$ in the Fe-based materials arises from the Fe 3$d$ orbitals. These systems are more complex than cuprates due to the contribution of multiple bands at $\epsilon_{F}$ and drastic changes in the Fermi surface topology arising from magnetic and structural transitions. For example, the two dimensional hole pockets at $\Gamma$-point in the tetragonal phase exhibit a transition to three dimensional topology at low temperatures ARPES-Liu ; ganesh-CaFe2As2 and in cT phase, they disappear due to the shift of the energy bands. The disappearance of the Fermi surface in cT phase rules out the possibility of Fermi surface nesting necessary for the transition to the SDW state. The Fe-based systems, in general, are found to be moderately correlated and inclusion of electron correlation strength, $U$ in the calculations renormalizes the energy bands without drastic change in the ground state properties dft_Diehl_prb_2014 ; dft_Kurmaev_prb_2009 ; ganesh-FeTeSe ; silke-arxiv ; silke-PRB . Jean Diehl et al.dft_Diehl_prb_2014 showed that the Hund’s coupling, $J$ ($\sim$0.8 eV) influences the electronic structure more strongly than the changes found due to $U$ ($\sim$ 4 eV). They also showed that upon transition from T to cT phase, the $d_{xy}$ bands become the least renormalized orbital from the highest renormalized one, which was attributed to increased hybridization of Fe 3$d_{xy}$ orbital with As 4$p_{x}$ and 4$p_{y}$ orbitals. The appearance of the cT phase under pressure and/or on suitable doping dft-Sanna-PRB12 ; sc_cafe2as2_saha_PRB_2012 has been attributed to the increased hybridization among As 4$p$ orbitals along the compressed $c$ axis. Yildrim et al. dft_yildrim_prl2009 found that the Fe spin state strongly depends on the $c$ axis reduction in cT phase. Based on fixed moment calculation, he showed that reduction by about half of the Fe magnetic moment corresponds to the $c$ axis collapse to the experimental value. Thus, there is a finite magnetic moment present at the Fe sites in the cT phase when it is not superconducting. Evidently, understanding of the importance of Fe moment in superconductivity is still at its infancy and needs further investigation. These results indicate that the electronic structure is very sensitive to the distance of As sites from the Fe plane, which is termed as ‘pnictogen height’. In the antiferromagnetic (AFM) orthorhombic (O) phase, the spins show antiparallel coupling along the longer axis, $a$ and parallel coupling along shorter axis, $b$. This has been theoretically found by many groups dft_yildrim_prl2009 ; dft_Deepa_njp_2009 but a satisfactory resolution of the preferred spin coupling along $a$ and $b$ axis is not yet found. First principle calculations under hydrostaticdft_Zhang_prb_2009 ; dft_Widom_PRB_2013 and uniaxialdft_Tomi_prb_2012 pressure conditions reproduced the first order transition to cT phase from AFM O/T phase accompanied by quenching of Fe magnetic moments and also rule out the possibility of an intermediate tetragonal phase in CaFe2As2 in contrary to various experiments neutron-kreyssing-PRB08 ; neutron_goldman_prb ; sc-p-park-JPCM08 , which is different from the behavior in BaFe2As2 at intermediate pressure range dft_Tomi_prb_2012 .

In the experimental front, the pressure induced effects in AFe2As2 (A = Ca, Ba, Sr) has been studied by many groups resulting into conflicting conclusions as discussed below. Superconductivity in BaFe2As2 & SrFe2As2 containing alkaline earth elements larger than Ca was found at a reasonably high pressure ($P$ $\approx$ 4 GPa and 3.2 GPa, respectivelyAlireza-JPCM09 ) with $T_{c}$ of about 29 K and 27 K sc_p_Awadhesh_Mani_epl_2009 ; sc_p_Hisashi_Kotegawa_JPSJ_2009 ; sc_p_Colombier_prb_2009 ; 34K_sc_p_Hisashi_Kotegawa_JPSG_2009 ; sc_p_Kazumi_Igawa_JPSG_2009 ; ct-sc-MKumar-PRB08 . On the other hand, CaFe2As2 shows superconductivity with $T_{c}$ upto 12 K in a low pressure range from 0.2 to 0.9 GPa sc-ct-torikachi-PRL08 although the smaller atomic size of Ca is expected to employ less strain into the system. A typical temperature-pressure phase diagram shown in Fig. 2 exhibits the experimental results obtained by Torikachvili et al. sc-ct-torikachi-PRL08 , which is consistent with the results in other studies sc_p_Torikachvili_prb_2009 . In Fig. 2, we observe that the transition to orthorhombic AFM phase is suppressed upon increasing pressure and disappears above $P$ = 0.35 GPa with the emergence of the cT phase. The transition temperature from T to cT phase gradually increases with pressure. The superconductivity appears at an intermediate $P$ range from 0.3 GPa to 0.9 GPa. In fact, superconductivity with $T_{c}$ $\approx$ 10 K at 0.69 GPa pressure has been detected by several groups sc-p-park-JPCM08 ; sc_p_lee_prb_2009 . Subsequent experiments, however, did not find superconductivity although cT phase could be reached on application of pressure. It was suggested that the earlier measurements carried out using convectional liquid medium clamp pressure cell might have made the pressure quasi hydrostatic due to freezing of the medium. These later results obtained using helium as a pressure medium question the relation of superconductivity to cT phase and in line with other views suggesting non-magnetic behavior in cT phase that rules out spin-fluctuation. Evidently, the observance/non-observance of superconductivity under pressure, non-observance of other structural phases under pressure, no superconductivity in ambient conditions, etc. are curious and remains to be an outstanding puzzle in this field.

Objective

A study of the reported resistivity data exhibiting superconductivity reveal an interesting scenario; the transitions in the resistivity data at $P$ = 0.23 GPa are significantly broad. On the other hand, the results exhibiting absence of superconductivity show sharp transitions throughout the pressure range involving the transitions from T to O, O to cT and T to cT phases. Moreover, the superconductivity always found to appear at a pressure where the system is in the proximity of the transition to cT phase or some volume fraction is already in cT phase. Thus, the origin of superconductivity in this material appears to have relation to the multi-crystallographic phases within the intermediate pressure regime. Nuclear magnetic resonance NMR_baek_prb ; NMR_kawasaki_iop , $\mu$ spin rotation msr_goko_prb , and neutron diffraction studies neutron-kreyssing-PRB08 ; neutron_goldman_prb ; sc_cafe2as2_Proke_PRB_2010 support the possibilities of coexisting O, cT and T phases. Some studies indicated signature of an unidentified third phase appearing in the intermediate pressure regime sc_p_lee_prb_2009 . While the magnetic order in single or multiple phase scenario can be attributed to the O phase, it is difficult to delineate the phase responsible for superconductivity. Clearly, superconductivity does not appear to be related to pure cT or O phase under pressure ct_sc_absence_Yu_PRB_2009 . Then the question is; how the superconductivity sets in CaFe2As2 under pressure? How non hydrostatic pressure helps to get supeconductivity while hydrostatic pressure does not? Superconductivity also appears via doping of charge carriers in these material - how different is that mechanism from the pressure induced one?

Here, we studied the detailed electronic structure of CaFe2As2 employing density functional theoretical (DFT) methods at varied conditions with an emphasis on the tendency of the system to evolve to a certain electronic structure due to change in structural parameters induced by the pressure. We observe interesting spectral weight transfer with the application of pressure leading to complex change in Fermiology of the system. The covalency between Ca, As and Fe states appear to play the key role in deriving the electronic properties of this system.

II Results

The electronic structure of a material is very sensitive to the change in crystal structure as the hybridization of various electronic states depends on the bond length, distortion in the lattice structure (leading to a change in bond angle) etc. Often it is found that covalency induces key interactions leading to distortion in the crystal structure Andersen-PRL ; CaSrRuO-PRB ; BaIrO3-PRB that derives the ground state properties.

Density of States

In Fig. 3, we show the calculated density of states (DOS) of CaFe2As2 in tetragonal phase. There are two distinctly separated energy regions in the figure contributed by As 4$s$ shallow core levels (-14 to -10 eV) and the valence electronic states above -6 eV. In the valence band regime, we observe signatures of hybridization between Fe 3$d$ and As 4$p$ states. Major contribution from As 4$p$ states appear in the bonding energy bands between -6 to -2 eV energies. The contribution for antibonding states appear above -2 eV and possess dominant Fe 3$d$ character. The valence electrons in the vicinity of the Fermi level, $\epsilon_{F}$ consist of essentially Fe 3$d$ states. The Fe-As hybridization leads to finite contribution from As 4$p$ states ($\sim$ 4% of the total contribution at $\epsilon_{F}$). Interestingly, the contribution of the Ca 3$d$ states are also found to be significant ($\sim$ 2.4%). This suggests that although Fe 3$d$ states play the dominant role in the electronic properties of CaFe2As2 as also found in earlier studies, the role of As 4$p$ and Ca 3$d$ states are non-zero due to covalency between these states with Fe 3$d$ states. The contribution from Ca 4$s$ states, however, is almost negligible ($\sim$ 0.2%) in the vicinity of $\epsilon_{F}$.

At ambient pressure, CaFe2As2 undergoes a magneto-structural transition from paramagnetic tetragonal to antiferromagnetic orthorhombic phase at 170 K. In Fig. 4, we show the calculated DOS for tetragonal and orthorhombic structures to investigate the effect in the electronic structure due to the above mentioned transitions. While the DOS at $\epsilon_{F}$ in the tetragonal structure is large and has a decreasing slope with the increase in energy, a pseudogap like feature appears in the magnetically ordered phase as shown in Fig. 4(b). $\epsilon_{F}$ is situated at the middle of the pseudogap exhibiting particle-hole symmetry in its close proximity. The position and width of the pseudogap observed here matches well with the LDA+DMFT results at ambient pressureDMFT_Yin_naturematerials_2011 (here, LDA stands for local density approximation and DMFT for dynamical mean field theory). While there is no apparent change in overall Fe 3$d$ bandwidth due to the transitions, significant change is observed in the details, such as the shift of the Fe 3$d$ spectral weight away from $\epsilon_{F}$ leading to significant broadening of the sharp feature above $\epsilon_{F}$. A new feature appear around -2 eV energy. The center of mass of the As 4$p$ contributions shifts to lower energy along with an enhanced contribution near $\epsilon_{F}$ indicating a stronger mixing of the Fe 3$d$ and As 4$p$ electronic states. Such enhanced mixing/covalency between Fe 3$d$ and As 4$p$ states might be the driving force for the magnetic and structural transitions in the ground state.

In Fig. 5, we superimpose the Fe 3$d$ partial density of states (PDOS) for collapsed tetragonal and tetragonal phases to have a comparative study. Fe atoms in tetragonal structure experience tetrahedral crystal field created by the As atoms. A regular tetrahedral crystal field splits the Fe 3$d$ levels into a doubly degenerate $e_{g}$ band and a triply degenerate $t_{2g}$ band. The distortions in the crystal structure of CaFe2As2 leads to a deviation from the regular tetrahedral field and the degeneracy of the crystal field split bands are lifted with $d_{xy}$ band slightly shifted above the doubly degenerate $d_{xz}$ & $d_{yz}$ bands. The contributing energy range and width derived by the hopping interaction terms remain very similar for all these bands as manifested in the figure. The $t_{2g}$ band consisting of $d_{xz}$, $d_{yz}$ and $d_{xy}$ orbitals lies close to $\epsilon_{F}$ and are the major contributors to the Fermi surface. The orbitals $d_{x^{2}-y^{2}}$ and $d_{z^{2}}$ possessing $e_{g}$ symmetry are situated far below the Fermi level.

In the cT phase, the As-Fe-As angle (= 106.13$\degree$) deviates from the regular tetrahedral angle at ambient condition ($\approx$ $109.47\degree$) and gives rise to additional distortions in the FeAs4 tetrahedra that enhances further the crystal filed splitting tetrahedra_screp_sarapov2014 . In Fig. 5(a) and (b), it appears that all the energy bands in cT phase represented by dashed lines are shifted towards lower energies relative to those in T phase. The energy bandwidth of $d_{xz}$ + $d_{yz}$ bands increases slightly in cT phase compared to that in T phase while the other $t_{2g}$ and $e_{g}$ bandwidths remain unchanged. This implies that in cT phase, the hybridization between Fe-As along $z$ direction enhances as expected due to the compression along $c$-axis; the in-plane orbitals remain almost unaffected. In T phase, major contribution at $\epsilon_{F}$ comes from ($d_{xz}$ + $d_{yz}$) and $d_{xy}$ electronic states. The $d_{xy}$ band shifts below $\epsilon_{F}$ in cT phase and thus, contribution at $\epsilon_{F}$ in cT phase arises primarily from degenerate ($d_{xz}$ + $d_{yz}$) orbitals. Evidently, the carrier density at $\epsilon_{F}$ is significantly less in the cT phase due to the shift of the energy bands, which is in contrast to the transport data indicating a decrease in electrical resistivity in the cT phase compared to the tetragonal phase sc-ct-torikachi-PRL08 . While more studies are required to understand this phenomena, the decrease in resistivity may be attributed to the enhancement of the mobility of charge carriers; the increase in hybridization of Fe ($d_{xz}$ + $d_{yz}$) and As 4$p$ states facilitates higher degree of the hopping of charge carriers.

The compression along $c$ axis in the cT phase influences the hybridization among As 4$p$ states too dft_yildrim_prl2009 leading to notable changes in their density of states. In Fig. 6(a) and (b), we show the PDOS corresponding to As 4$p_{z}$ and the degenerate As (4$p_{x}$+4$p_{y}$) energy bands. The results corresponding to the orthorhombic phase shown in the figure were calculated without considering magnetic order to study the effect on electronic structure due to the change in structural parameters only. From the figure, it is clear that all the As 4$p$ bands in O and T phases are essentially identical while the $p$ bands in cT phase exhibit a shift towards lower energies relative to those in T and O phases - a trend similar to the findings in Fe 3$d$ cases. In addition, there are several subtle changes observed for various $p$ orbitals. For example, the overall width of the As 4$p_{z}$ is somewhat larger in the cT phase compared to T phase due to enhancement of DOS near the band edges suggesting an effective increase in hybridization among the $p_{z}$ orbitals. The energy distribution of the density of states corresponding to the $p_{x}$ and $p_{y}$ electronic states exhibit spectral weight shift near the top of the As 4$p$ bands (see around -2 eV in the figure). The unoccupied part appears to be close to each other in all the three cases studied. The changes in the bonding bands between -6 to -2 eV indicate signature of enhancement of covalency in the system in the cT phase.

The results for the Ca 3$d$ and 4$s$ PDOS is shown in the Fig. 6(c) and (d), respectively. Ca 3$d$ states contribute primarily above $\epsilon_{F}$ (energy range 3 to 6 eV) with small contributions at and below $\epsilon_{F}$. The intensity of Ca 4$s$ states is distributed over larger energy range. While the major contribution of 4$s$ states appear between 3 to 8 eV above $\epsilon_{F}$, there is finite contribution below -3 eV energy. Interestingly, the contributions below $\epsilon_{F}$ exhibit a shift towards lower energies in the cT phase as found in Fe and As cases keeping the PDOS above $\epsilon_{F}$ within almost the same energy range in all the structural phases. In the crystal structure, the Ca layers are sandwiched between two As layers and the changes in Ca states below $\epsilon_{F}$ can be attributed to the hybridization of Ca 4$s$ states with the electronic states associated to As in the valence band. The compression ct_vanroekeghem_prb_2016 along $c$-axis due to the transition to the cT phase reduces the separation between the Ca and As atomic layers and the hybridization between Ca 4$s$ and As 4$p$ will be enhanced. It is clear that although the Ca contribution in the valence band regime is small, it is difficult to ignore the role of Ca in the structural changes and their electronic structure as well casrruo-core-PRB07 ; casrruo-core-SL14 . Almost all the high temperature (unconventional) superconductors exhibit effective two dimensional electronic structure, where the conduction sheet is separated by an insulating layer, often called ’charge reservoir’ layer. It is believed that the cations with $s$ electrons in the charge reservoir layer plays an important role in the pairing interactions, renormalizing disorder induced effects, etc charge-reservoir ; charge-reservoir_book . In the cuprate superconductors, elements like La, Hg, Pb, Bi form the charge reservoir layers and the Cu-O layers form the two dimensional conduction layers. In the present case, Ca layer sandwiched between conducting FeAs layers appears to play similar role and hence enhanced interaction of Ca 4$s$ states in the cT phase appears to have great implication for the onset of superconductivity in these systems.

The results for cT phase discussed so far correspond to the lattice constants found experimentally for the pressure of 0.63 GPa. Now, we study the density of states at different pressures utilizing the structural parameters obtained experimentally. The DOS calculated for tetragonal, collapsed tetragonal phases at $P$ = 0.47 GPa (cTp1) and $P$ = 0.63 GPa (cTp2) are superimposed in Fig. 7. The Fe 3$d$ PDOS in the non-magnetic phase shown in Fig. 7(a) exhibit significant intensity at $\epsilon_{F}$ in the tetragonal phase. With the application of pressure, Fe 3$d$ PDOS shift to lower energies reducing the intensity at $\epsilon_{F}$ to almost zero. Interestingly, the results for cTp1 and cTp2 appear identical indicating less dependence of the non magnetic phase on pressure once the pseudogap like feature at $\epsilon_{F}$ is present in the electronic structure. The overall bandwidth appears to be similar in every case.

The magnetic solution shown in Fig. 7(b) exhibits large separation between the up and down spin bands in the tetragonal phase with an exchange splitting of the order of 2 eV. With the application of pressure, up and down spin Fe 3$d$ PDOS shift towards each other reducing the exchange splitting with the increase in pressure. This change, however, is not monotonic; the data corresponding to 0.47 GPa shows much larger shift and smaller exchange splitting compared the case for $P$ = 0.63 GPa. Since the bottom of the valence band remains almost same in every case, the change at the top of the band also leads to an increase in bandwidth. Thus, the width of the valence band for $P$ = 0.47 GPa appears to be the largest. The transition to the magnetically ordered phase in tetragonal structure leads to a large dip ($pseudogap$) in DOS at $\epsilon_{F}$. Interestingly, the dip in the magnetically ordered phase reduces with the application of pressure in cT phase as shown in the inset of Fig. 7(b), where sum of the up and down spin Fe 3$d$ contributions are plotted. This indicates poorer degree of stabilization of the magnetic phase in cT structures and the electronic properties at $P$ = 0.47 GPa is closer to the proximity of instability in magnetic order relative to the other cases.

In order to study the difference between magnetic and non-magnetic solutions in more detail, we superimpose the results of magnetic and non-magnetic calculations for Fe 3$d$ and As 4$p$ PDOS in Fig. 8 for tetragonal structure. Fe 3$d_{z^{2}}$ spin-orbitals exhibit two features separated by about 2 eV in the non-magnetic phase. Fe 3$d_{x^{2}-y^{2}}$ exhibit similar structures with somewhat larger separation ($\sim$ 2.5 eV). Both these bands do not contribute at $\epsilon_{F}$. Fe 3$d_{xy}$ and (3$d_{xz}$ + 3$d_{yz}$) bands possess one intense feature in each case and appear in the vicinity of $\epsilon_{F}$ providing dominant contribution at $\epsilon_{F}$. In the magnetically ordered phase, each band splits into up and down spin bands due to the exchange interactions. Such splitting is found to be about 2 eV in 3$d_{z^{2}}$ and 3$d_{x^{2}-y^{2}}$ bands, while it is about 2.5 eV for all other cases. All the up spin bands appear below $\epsilon_{F}$ and the down spin states having 3$d_{xy}$, 3$d_{xz}$ and 3$d_{yz}$ symmetries primarily contribute at the Fermi surface indicating an essentially half-metallic scenario. As 4$p$ states exhibit interesting scenario in magnetically ordered phase. In addition to a small enhancement of As 4$p_{z}$ contributions at $\epsilon_{F}$, these states move towards lower energies in the magnetically ordered phase similar to the scenario in Fe 3$d$ case and exhibit small exchange splitting induced by Fe 3$d$ moments. In sharp contrast, the contribution of As (4$p_{x}$ + 4$p_{y}$) at $\epsilon_{F}$ in the magnetic phase is almost vanished, the energy bands move in opposite direction (towards higher energies as shown in Fig. 8(f)) and does not show exchange splitting. Absence of exchange splitting may be apprehended due to the absence of hybridization of these states with the Fe 3$d$ states (As layers are situated above and below the Fe layers). However, energy shift in opposite direction is curious.

The scenario is somewhat different in cT phase - the results for cTp2 phase is shown in Fig. 9. While the overall structure of the PDOS appear quite similar to the tetragonal case, there are significant differences. For examples, Fe 3$d_{z^{2}}$ and 3$d_{xy}$ PDOS exhibit exchange splitting of about 1 eV, which is much smaller than that in T phase. The behavior of 3$d_{x^{2}-y^{2}}$ and 3$d_{xz}$ + 3$d_{yz}$ PDOS, however, appear similar to the tetragonal case. Interestingly, the contrasting behavior of As 4$p_{z}$ and (4$p_{x}$ + 4$p_{y}$) bands is not present in this case - $both$ the bands show small exchange splitting and an energy shift towards lower energies. A compression along the $c$ axis in the cT phase brings Fe and As atomic layers closer to each other and hence, the electronic states within As layer ($xy$-plane) also experiences Fe 3$d$ moment induced interactions. Here, the contribution of up spin (4$p_{x}$ + 4$p_{y}$) bands at $\epsilon_{F}$ becomes finite, while it was absent in the non-magnetic case. It is to note here that the down spin contributions corresponding to 3$d_{xy}$ and 3$d_{xz}$ + 3$d_{yz}$ states are essentially unoccupied with $\epsilon_{F}$ pinning at the bottom of the conduction band. However, these bands are significantly occupied in cT phase with $\epsilon_{F}$ appearing close to the peak of the bands.

II.1 Fermi surface

In order to understand the Fermiology and its link to different crystal structures, the Fermi surfaces (FS) of CaFe2As2 are investigated below. Fig.10 depicts the Fermi surface in the tetragonal structure exhibiting quasi 2$D$ symmetry with the cylindrical axis as $k_{z}$-axis. In these calculations, $z_{As}$ is to optimize the structural parameters. It consists of four distinct Fermi sheets - two of them are wrapped around the $k_{z}$ axis at the brillouin zone center, $\Gamma$ and represent the hole states. Two Fermi sheets centered at the corner of the Brillouin zone are the electron pockets. The spin density wave type magnetic order observed in CaFe2As2 is associated to the Fermi surface nesting; a hole pocket centered around $\Gamma$-point is nested with an electron pocket around $M$-point. The dimensionality of the Fermi surface plays an important role in FS nesting; the 2$D$ behavior enhances the probability of nesting because parallel hole and electron sheets results in particle hole excitation and manifests as a Fermi surface nesting. The nesting wave vector $Q$ along $\Gamma-M$ direction in the primitive tetragonal unit cell is shown in Fig. 11. Since, the nesting vector $Q$ is almost equal for both the hole sheet, both contribute in magnetic ordering.

The dimensionality and the topology of the Fermi surface depends on $z_{As}$ that is linked to the pnictogen height. The Fermi surface calculated with the experimental parameters (without employing relaxation in the structure) is shown in Fig. 12. In this case, the Fermi surface consists of five sheets; the additional hole sheet appears at the zone center due to lifting of the degeneracy of one of the hole Fermi sheets in the relaxed structure. One of the electron sheet spreads around the perimeter of the Brillouin zone, the other remains unchanged. The change in the shape of the Fermi sheets at the zone corners and the zone center leads to the disappearance of the nesting in one of the hole sheets.

From various calculations, we observe that the electronic structure corresponding to the experimentally found crystal structure is significantly different from the one obtained for the optimised structure within the model calculation via relaxing As positions. The relaxed structure favors smaller $z_{As}$ value and hence, smaller pnictogen height compared to the experimentally observed values and larger degree of Fermi surface nesting, thereby enhancing the degree of magnetic order. The real material, however, forms with relatively larger pnictogen height and one nesting vector as also observed in ARPES measurements of various materials in similar class. One reason for such difference could be related to the underestimation of correlation induced effects in the density functional calculations - the system prefers to reduce the hybridization and hence itineracy via increasing the pnictogen height. The other possibility could be a proximity of the ground state crystal structure of CaFe2As2 to a metastable state and thus, the system gets arrested to different strained crystal structure during sample preparation instead of going to the ground state configuration. The finding of the change of volume fraction of different structures and stains in SrFe2As2 with preparation conditions Saha-PRL-SrFe2As2 could be a realization of this conjecture. The experimental findings of contrasting behavior in CaFe2As2 is also in line with this view.

In Fig. 13, we investigate the Fermi surfaces in all the structural phases realized in CaFe2As2 by varying parameters like temperature, pressure etc. In the O phase, without any magnetic order, the Fermi surface looks very similar to the T phase. The difference in the definition of $k$ vectors appears due to their Brillouin zone orientation. It is evident from the comparison of Figs. 13(a) and 13(b) that the Fermi sheets in the orthorhombic phase are shrinked compared to those in the tetragonal phase as evident from the larger separation between electron and hole sheets. The orthorhombic phase emerges at low temperature presumably due to thermal compression leading to a small distortion of the lattice in the $xy$ plane and small ($\sim$ 5%) compression in the $c$ axis. Thus, the pnictogen height becomes smaller as reflected in the Fermiology.

Fermi sheets in the O phase with AFM ordering are shown in Fig. 13(d). Here, the Fermi sheets look three dimensional and very different from all other cases due to the opening of energy gap at $\epsilon_{F}$ induced by the nesting of the electron sheets at the zone corner with the hole sheets at zone center.

Notable difference appears in the cT structure; the Fermi surface shown here correspond to the structural parameters at $P$ = 0.67 GPa with $z_{As}$ value slightly smaller than the case shown in Fig. 10. We chose these parameters corresponding to slightly higher pressure than the previous cases to demonstrate the changes in Fermi surface topology with greater clarity. The 2$D$ nature of the hole Fermi surfaces around the $z$-axis is lost completely in the results for the cT phase. Here, all the three hole pockets around the $\Gamma$-point disappeared. Only one hole sheet survives around the $Z$ point making it a closed one with larger radius. The electron sheets are not affected significantly. Such drastic change in the hole Fermi surface around $k_{z}$-axis rules out the electron-hole nesting responsible for magnetic long range order and is consistent with the experimental finding of absence of magnetic order in this phase.

II.2 Energy bands

We now turn to the investigation of the properties of the energy bands. While the density of states discussed above provide information about the contribution of various electronic states as a function of energy, these results will help to understand the $k$-resolved information, which are often helpful to understand the electronic properties and can also be observed in angle resolved photoemission measurements. In Fig. 14, we show the energy bands of CaFe2As2 for the tetragonal structure; the thickness of the line represents the orbital contributions. It is evident that the Fe 3$d$ symmetry adapted bonding bands appearing between -6 to -2 eV possess dominant As 4$p$ character. The bands close to the Fermi level are primarily constituted by Fe 3$d$ states with notable contributions from As 4$p$ states. Interestingly, while the bands around $\Gamma$ point forming the hole pocket in the Fermi surface possess essentially Fe 3$d$ character, one of the electron pocket around $X$-point possess large As 4$p$ character.

The detailed orbital symmetries near $\epsilon_{F}$ are shown in Fig. 15. At $\Gamma$, three energy bands are observed to cross $\epsilon_{F}$ - these bands are denoted by $\alpha$, $\beta$ and $\gamma$ as shown separately in Fig. 15(g) for better clarity. These three energy bands create three hole pocket around $\Gamma$ ARPES-Liu . There are two electron pockets around the high symmetry point, $X$. It is clear from Fig. 15(c) and Fig. 15(d) that the $\gamma$ band possesses mainly $d_{xy}$ character and the $\alpha$ & $\beta$ bands have contribution from degenerate $d_{xz}$ and $d_{yz}$ orbitals ganesh-CaFe2As2 . These three bands possess $t_{2g}$ symmetry. The $e_{g}$ bands constituted by $d_{x^{2}-y^{2}}$ and $d_{z^{2}}$ spin-orbitals appear far away from the Fermi level and hence do not play significant role in the electronic properties of the system. Overall, the bands near the Fermi level has very less As 4$p$ character except one electron pocket at $X$, which has As 4$p_{z}$ character. This is in line with the expectation from the fact that As layers appear above and below Fe layers and hence $p_{z}$ orbital is expected to play a major role in the Fe 3$d$-As 4$p$ hybridization.

In Fig. 16, we show the change in band structure in different phases along $X-\Gamma-Z$ symmetry lines. In collapsed tetragonal phase, the bands at $\Gamma$ shifts below the Fermi level and the hole pockets around $\Gamma$ disappear. The nonmagnetic calculation of the orthorhombic structure gives rise to similar band structure like tetragonal phase, but if we consider the antiferromagnetic configuration then the band structure becomes complex. The hole pocket at $\Gamma$ nesting with an electron pocket at $X$ gives rise to an energy gap at $\epsilon_{F}$ via folding of the energy bands at the zone boundary created by the magnetic supercell.

III Discussion

From the above results, it is clear that electronic structure of CaFe2AS2 is very sensitive to the crystal structure and pnictogen height is an important parameter in deriving the ground state properties. The system appears to favor smaller pnictogen height although the experiments found structures with somewhat higher values and significant dependence of the structure on preparation conditions. On application of pressure, the tetragonal symmetry gets arrested till the low temperatures and the electronic properties become significantly different from magnetically ordered phase observed in ambient conditions.

In cT phase, the $c$-axis collapse reduces Fe moment and Fe-As hybridization along $z$ direction enhances dft_yildrim_prl2009 . In fact, the change in Fe-As hybridization and consequent change in the magnetic moment of Fe turn out to be the key factors in deciding the magnetic and structural properties of CaFe2As2. To get a clear idea about this we have calculated the local magnetic moment at Fe-sites. The calculation are done for tetragonal phase and collapsed tetragonal phase. In the later case, we choose the lattice parameters corresponding to two different applied pressures. The phase realized under 0.47 GPa pressure is termed as cTp1 and the cT phase with $P$ = 0.63 GPa is called cTp2. Experimentally, the magnetization in the cT phase has been studied extensively and no signature of long range magnetism was found. There are controversy on the magnetic moment though - while Ma Long et al. MaLong-ChinesePhysB13 found signature of quenching of Fe moment on transition to cT phase induced by Pr substitution at Ca sites, magnetic susceptibility data by W. Wu et al. ct_sc_absence_Yu_PRB_2009 indicate finite moment of Fe in the cT phase though depleted from its value in ambient tetragonal phase. Moreover, the Fe magnetic moment appears to be sample dependent. We calculated Fe moments assuming parallel alignment of Fe-moments and compare the results with their values corresponding to the non magnetic phase in order to investigate the trend of the change of converged energies with the change in external parameters, which is often captured well employing DFT methodologies Andersen-PRL ; CaSrRuO-PRB ; BaIrO3-PRB . The table below gives the list of lattice parameters used for the calculations along with the calculated magnetic moment and ground state energies.

In the tetragonal structure, the converged energy for the magnetically ordered configuration is about 150 meV lower than that for the nonmagnetic (NM) configuration suggesting a preference for the magnetic order in the ground state. The magnetic moment of Fe atoms is found to be 2.37 $\mu_{B}$ in this structure. The cT structure manifests a contrasting scenario exhibiting significantly higher energy for magnetically ordered phase relative to the NM phase. Application of pressure leads to a gradual reduction of the lattice constant, $c$ and As-As bondlength ($d_{AsAs}$). The magnitude of the lattice constant, $a$ and $z_{As}$ exhibit anomalous scenario with the largest and smallest values at $P$ = 0.47 GPa, respectively. With the application of pressure, the nonmagnetic phase appears to be favored more. The ground state energy of the non-magnetic phase at cTp1 structure is found to be the lowest relative to all the calculations we have done.

The Fe moment appears to diminish with the application of pressure. Interestingly, the lowest moment of 1.95 $\mu_{B}$ is found for the pressure of 0.47 GPa (cTp1 phase) and then the Fe moment increases with further increase in pressure although the compression of $c$ is more at cTp2 phase compared to cTp1 phase. The value of $z_{As}$, however, is smaller in the cTp1 phase indicating a direct link of the Fe-moment to $z_{As}$ value. It is still an open question whether the magnetic moment at Fe sites actually vanishes or there is a small magnetic moment on Fe atom, but they don’t give any long range order dft_yildrim_prl2009 . All our calculations indicate that the magnetic configuration in the cT structure has higher energy compared to the non-magnetic phase. The magnetic moment decreases as the Fe-As bond length decreases. This can be understood considering the fact that decrease in Fe-As bondlength makes the conduction electrons more itinerant and hence less local moment. It is puzzling that the Fe-As bondlength is smallest in cTp1 structure and a small decrease in Fe-As bond length from 2.3405 Å in cTp2 structure to 2.338 Å in cTp1 structure changes the magnetic moment from $1.95\mu_{B}$ to $1.37\mu_{B}$.

The above results establish that the Fe spin state is highly sensitive to the Fe-As hybridization, which is linked to the pnictogen height too. The shorter As-As distance leads to increased hybridization within the As-plane and higher degree of As $p$ electron itineracy. However, the local character of the Fe $d$ electrons remain less affected. On the other hand, higher Fe-As hybridization leads to larger covalent splitting of the Fe-As bonding and antibonding bands leading to a shift of the valence bands. Some of the bands shifts below $\epsilon_{F}$ and become completely filled. The energy bands crossing the Fermi level will have higher itineracy and hence less local character, which is presumably one of the reasons for the reduction of Fe moment and the resistivity.

The calculated Fe-As bond length in the AFM ground state with orthorhombic structure is 2.3654 Å, which is larger than that in cT phase and shorter than the one in T phase. The Fe spin moment is found to be 1.87 $\mu_{B}$, which is much larger than the experimental value of 0.80 $\mu_{B}$. Various theoretical studies tried to capture the experimental moment by varying $z_{As}$ (Experimental value = 0.3664). It is observed that the spin moment becomes closer to the experimental value for the $z_{As}$ values of 0.3567 structure_Tompsett_physicab2009 . It is clear that $z_{As}$ and hence the pnictogen height is an important parameter leading to varied electronic properties in this system. Moreover, these results establish that the cT phase possesses finite magnetic moment as also found experimentally in magnetic susceptibility measurements for almost all kinds of phases prepared so far. The absence of long range order appears to be due to its proximity to the critical point of loosing Fermi surface nesting & no magnetic order ganesh-cr ; cao-casrruo ; maiti-casrruo .

Interestingly, coupling of Ca 4$s$ states to the conduction layer becomes stronger in the cT phase as discussed in Fig. 6 and redrawn here to emphasize the scenario. In Fig. 17, we show the Ca 4$s$ PDOS along with Fe 3$d$ and As 4$p$ contributions for different structures. Ca 4$s$ contributions in orthorhombic and tetragonal phases look almost similar with major intensities above $\epsilon_{F}$, a reasonably intense feature between -6 to -3 eV and almost no intensity in the vicinity of $\epsilon_{F}$. While the intensity above $\epsilon_{F}$ remains dominant in the cT phase too like other cases, there is significant redistribution of PDOS in the lower energy regime. Most importantly the contribution in the vicinity of $\epsilon_{F}$ becomes almost double of the contributions in other two structures. The contributions in the energy range -3 to -6 eV appear due to the Fe-As bonding bands with dominant As 4$p$ character and -2.5 to 0 eV is contributed by anti-bonding bands with dominant Fe 3$d$ character. The intense feature between -2.5 to -4 eV in As 4$p$ PDOS are the non-bonding As 4$p$ states. The features in Ca 4$s$ PDOS resemble well with structures in the As 4$p$ PDOS as these states emerges due to the hybridization of Ca 4$s$ states with As 4$p$ states. Thus, the enhancement of 4$s$ contributions below $\epsilon_{F}$ can be attributed to enhanced hybridization of Ca and As states due to compression of $c$ axis in the cT phase.

Charge reservoir layers, generally insulating in nature, are believed to play an important role in superconductivity charge-reservoir ; charge-reservoir_book as these layers help to preserve two dimensionality of the system and can protect the conduction layers via screening of substitutional and/or other type of disorder induced effects. Effective two dimensional electronic structure seem to be essential to achieve high superconducting transition temperature and almost all the superconductors (conventional and unconventional) having high transition temperature possess two dimensional electronic structure. Thus, the higher degree of Ca-As hybridization in collapsed tetragonal phase and the observation of pressure induced superconductivity edoped_sc_cafe2as2_Lv_pnas_2011 ; sc_cafe2as2_Jeffries_PRB_2012 ; sc_cafe2as2_Neeraj_PRB_2009 ; NMR_kawasaki_iop appears to be linked.

IV Conclusion

In summary, we present here the detailed electronic structure of CaFe2As2 for different crystallographic structures using density functional theory. The calculated density of states indicate strong covalency induced effect in the electronic structure, which increases with the compression of $c$ axis observed in collapsed tetragonal phase. While the energy distribution of density of states appear similar in tetragonal and orthorhombic phase in ambient conditions, a significant shift towards lower energies is observed in collapsed tetragonal phase, which reduces the total energy of the system. The hole pockets at the $\Gamma$ point eventually vanishes at higher pressure ruling out the Fermi surface nesting related to the formation of spin density wave states. The electron pockets centered around the high symmetry points $M$/$X$/$Z$ evolves significantly in the cT phase. The enhancement of Fe-As hybridization in the cT phase presumably leads to higher degree of itineracy relevant for the changes in electrical resistivity in this system. Interestingly, the changes in Ca 4$s$ partial density of states suggest enhancement of Ca-As hybridization with the decrease in $c$, which is significant in the context of superconductivity as it is believed that such layers act as charge reservoir layer in various high temperature superconductors. We observe that the changes in the Fermiology and energy band shifts due to application of pressure depletes the magnetic moment at Fe-sites although they do not show signature of magnetic order indicating the proximity of such strained condition to the quantum fluctuations.

V Method

Electronic structure calculations were carried out using full potential linearized augmented plane wave method (FLAPW) captured in the Wien2k software wien2k . Convergence was achieved by fixing the energy convergence criteria to 0.0001 Rydberg (1 meV). For nonmagnetic calculations, we have used 10$\times$10$\times$10 $k$-points in the Brillouin zone and for Fermi surface calculations 39$\times$39$\times$10 $k$-points were used. We have used the Perdew-Burke-Ernzerhof generalized gradient approximation (GGA)dft_gga for our density functional theoretical calculation dft_ks . The Fermi surfaces were calculated using Xcrysdencrysden . The collapsed tetragonal phase of CaFe2As2 achieved via application of pressure/varying annealing conditions, possesses the same space group but a reduced $c$ axis and slightly increased $a$ axis. The volume of the unit cell reduces effectively. The lattice parameters at $P$ = 0.63 GPa are $a$ = 3.9780(1)Å, $c$ = 10.6073(7)Å, and $z_{As}$ = 0.372(1)neutron-kreyssing-PRB08 . Orthorhombic structure of CaFe2As2 appearing at low temperature and ambient pressure has $Fmmm$ space group with the lattice parameters $a$ = 5.506(2)Å, $b$ = 5.450(2)Å, $c$ = 11.664(6)Å, and $z_{As}$ = 0.36642(5)neutron-kreyssing-PRB08 . For Antiferromagnetic calculations, we used 10$\times$10$\times$10 $k$-points in the Brillouin zone and for Fermi surface 23$\times$23$\times$10 $k$-points.

Acknowledgements

K. M. acknowledges financial assistance from the Department of Science and Technology, Govt. of India under J.C. Bose Fellowship program and Department of Atomic Energy, Govt. of India.

Author Contributions

K.A. carried out all the calculations. K.M. initiated the study and supervised the project. K.A. and K.M. jointly analyzed the data and prepared the manuscript.

COMPETING FINANCIAL INTERESTS

The authors declare no competing financial interests.