Emergent electronic structure of CaFe2As2

CaFe2As2 exhibits collapsed tetragonal (cT) structure and varied exotic behaviour 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 that the stability of the cT phase under pressure arises from the enhancement in hybridization induced effects and shift of the energy bands towards lower energies. The Fermi surface centered around Γ 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.


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 2,[9][10][11][12] . These systems are more complex than cuprates due to the contribution of multiple bands at ε F and drastic changes in the Fermi surface topology arising from magnetic and structural transitions. For example, the two dimensional hole pockets at Γ-point in the tetragonal phase exhibit a transition to three dimensional topology at low temperatures 13,14 and in the cT phase, they disappear due to the shift of the energy bands. The disappearance of the Fermi surface in the 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 [15][16][17][18][19] . Jean Diehl et al. 15 showed that the Hund's coupling, J (~0.8 eV) influences the electronic structure more strongly than the changes found due to U (~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 3d xy orbital with As 4p x and 4p y orbitals. The appearance of the cT phase under pressure and/or on suitable doping 8,20 has been attributed to the increased hybridization among As 4p orbitals along the compressed c axis. Yildrim et al. 9 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 9, 21 but a satisfactory resolution of the preferred spin coupling along a and b axis is not yet found. First principle calculations under hydrostatic 12,22 and uniaxial 23 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 CaFe 2 As 2 in contrary to various experiments 3,24,25 , which is different from the behaviour in BaFe 2 As 2 at intermediate pressure range 23 .
On the experimental front, the pressure induced effects in AFe 2 As 2 (A = Ca, Ba, Sr) has been studied by many groups resulting into conflicting conclusions as discussed below. Superconductivity in BaFe 2 As 2 & SrFe 2 As 2 containing alkaline earth elements larger than Ca was found at a reasonably high pressure (P ≈ 4 GPa and 3.2 GPa, respectively 26 ) with T c of about 29 K and 27 K 5, 27-31 . On the other hand, CaFe 2 As 2 shows superconductivity with T c upto 12 K in a low pressure range from 0.2 to 0.9 GPa 7 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. 7 , which is consistent with the results in other studies 32 . 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 ≈ 10 K at 0.69 GPa pressure has been detected by several groups 24,33 . 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 citing the fact that non-magnetic behaviour in cT phase 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 remain to be outstanding puzzles 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 34,35 , muon spin rotation (μSR) 36 , and neutron diffraction studies 3, 25, 37 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 33 . 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 38 . Then the question is; how the superconductivity sets in CaFe 2 As 2 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 CaFe 2 As 2 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.

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 39-41 that derives the ground state properties. Below, we discuss the calculated results for different structural configurations of CaFe 2 As 2 . It is to note here that the calculated results in various cases are shown with respect to their Fermi energy. Since the Coulomb potential can be different in different calculations, such alignment of the energy levels may be questionable. In the present case, the ground states in all the configurations of CaFe 2 As 2 are metallic with no change in electron count of the system. Moreover, all the experimental (photoemission) data are procured after aligning the Fermi energy of the material with the Fermi energy of the instrument via establishing an electrical contact. Therefore, such alignment in the present case appears reasonable and represents results corresponding to their experimental counterparts. Fig. 3, we show the calculated density of states (DOS) of CaFe 2 As 2 in tetragonal phase.

Density of States. In
There are two distinctly separated energy regions in the figure contributed by As 4s 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 3d and As 4p states. Major contribution from As 4p 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 3d character. The valence electrons in the vicinity of the Fermi level, ε F consist of essentially Fe 3d states. The Fe-As hybridization leads to finite contribution from As 4p states (~4% of the total contribution at ε F ). Interestingly, the contribution of the Ca 3d states are also found to be significant (~2.4%). This suggests that although Fe 3d states play the dominant role in the electronic properties of CaFe 2 As 2 as also found in earlier studies, the role of As 4p and Ca 3d states are non-zero due to covalency between these states with Fe 3d states. The contribution from Ca 4s states, however, is almost negligible (~0.2%) in the vicinity of ε F . At ambient pressure, CaFe 2 As 2 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 ε 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). ε 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 pressure 42 (here, LDA stands for local density approximation and DMFT for dynamical mean field theory). While there is no apparent change in overall Fe 3d bandwidth due to the transitions, significant change is observed in the details, such as the shift of the Fe 3d spectral weight away from ε F leading to significant broadening of the sharp feature above ε F . A new feature appear around −2 eV energy. The center of mass of the As 4p contributions shifts to lower energy along with an enhanced contribution near ε F indicating a stronger mixing of the Fe 3d and As 4p electronic states. Such enhanced mixing/covalency between Fe 3d and As 4p states might be the driving force for the magnetic and structural transitions in the ground state.
In Fig. 5, we superimpose the Fe 3d 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 3d levels into a doubly degenerate e g band and a triply degenerate t 2g band. The distortions in the crystal structure of CaFe 2 As 2 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  In the cT phase, the As-Fe-As angle (=106.13°) deviates from the regular tetrahedral angle at ambient condition (≈109.47) and gives rise to additional distortions in the FeAs 4 tetrahedra that enhances further the crystal filed splitting 43 . 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 ε F comes from (d xz + d yz ) and d xy electronic states. The d xy band shifts below ε F in cT phase and thus, contribution at ε F in cT phase arises primarily from degenerate (d xz + d yz ) orbitals. The carrier density at ε 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 7 . 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 4p states facilitates higher degree of the hopping of charge carriers.
The compression along c axis in the cT phase influences the hybridization among As 4p states too 9 leading to notable changes in their density of states. In Fig. 6(a) and (b), we show the PDOS corresponding to As 4p z and the degenerate As (4p x + 4p 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 4p 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 3d cases. In addition, there are several subtle changes observed for various p orbitals. For example, the overall width of the As 4p 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 4p 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 3d and 4s PDOS is shown in the Fig. 6(c) and (d), respectively. Ca 3d states contribute primarily above ε F (energy range 3 to 6 eV) with small contributions at and below ε F . The intensity of Ca 4s states is distributed over larger energy range. While the major contribution of 4s states appear between 3 to 8 eV above ε F , there is finite contribution below −3 eV energy. Interestingly, the contributions below ε F exhibit a shift towards lower energies in the cT phase as found in Fe and As cases keeping the PDOS above ε 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 ε F can be attributed to the hybridization of Ca 4s states with the electronic states associated to As in the valence band. The compression 44 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 4s and As 4p 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 45,46 . 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 47,48 . 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 4s 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 (cT p1 ) and P = 0.63 GPa (cT p2 ) are superimposed in Fig. 7. The Fe 3d PDOS in the non-magnetic phase shown in Fig. 7(a) exhibit significant intensity at ε F in the tetragonal phase. With the application of pressure, Fe 3d PDOS shift to lower energies reducing the intensity at ε F to almost zero. Interestingly, the results for cT p1 and cT p2 appear identical indicating less dependence of the non magnetic phase on pressure once the pseudogap like feature at ε 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 3d 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 ε 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 3d 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 3d and As 4p PDOS in Fig. 8 for tetragonal structure.  the exchange interactions. Such splitting is found to be about 2 eV in d 3 z 2 and − d 3 x y 2 2 bands, while it is about 2.5 eV for all other cases. All the up spin bands appear below ε F and the down spin states having 3d xy , 3d xz and 3d yz symmetries primarily contribute at the Fermi surface indicating an essentially half-metallic scenario. As 4p states exhibit interesting scenario in magnetically ordered phase. In addition to a small enhancement of As 4p z contributions at ε F , these states move towards lower energies in the magnetically ordered phase similar to the scenario in Fe 3d case and exhibit small exchange splitting induced by Fe 3d moments. In sharp contrast, the contribution of As (4p x + 4p y ) at ε 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 3d 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 cT p2 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 d 3 z 2 and 3d xy PDOS exhibit exchange splitting of about 1 eV, which is much smaller than that in T phase. The behaviour of − d 3 x y 2 2 and 3d xz + 3d yz PDOS, however, appear similar to the tetragonal case. Interestingly, the contrasting behaviour of As 4p z and (4p x + 4p 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 3d moment induced interactions. Here, the contribution of up spin (4p x + 4p y ) bands at ε F becomes finite, while it was absent in the non-magnetic case. It is to note here that the down spin contributions corresponding to 3d xy and 3d xz + 3d yz states are essentially unoccupied with ε F pinning at the bottom of the conduction band. However, these bands are significantly occupied in cT phase with ε F appearing close to the peak of the bands.
Fermi surface. In order to understand the Fermiology and its link to different crystal structures, the Fermi surfaces (FS) of CaFe 2 As 2 are investigated below. Figure 10 depicts the Fermi surface in the tetragonal structure   exhibiting quasi 2D symmetry with the cylindrical axis as k z -axis. In these calculations, z As is varied 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, Γ 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 CaFe 2 As 2 is associated to the Fermi surface nesting; a hole pocket centered around Γ-point is nested with an electron pocket around M-point. The dimensionality of the Fermi surface plays an important role in FS nesting; the 2D behaviour 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 Γ − 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 favours 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 observed in ARPES measurements 13,14 . 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 in a correlated system. The other possibility could be the proximity of the ground state crystal structure of CaFe 2 As 2 to a metastable state and thus, the system gets arrested to different strained crystal structure during the sample preparation instead of going to the ground state configuration. The finding of the change of volume fraction of different structures and stains in SrFe 2 As 2 with preparation conditions 49 could be a realization of this conjecture. The experimental findings of the contrasting behaviour in CaFe 2 As 2 is also in line with this view.
In Fig. 13, we investigate the Fermi surfaces in all the structural phases realized in CaFe 2 As 2 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 Fig. 13(a) and (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 (~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 ε 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.63 GPa with z As value slightly smaller than the experimentally found value in tetragonal structure. We chose the parameters corresponding to this higher pressure case to demonstrate the changes in the Fermi surface topology with greater clarity. The 2D 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 Γ-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 the absence of magnetic order in this phase.

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 CaFe 2 As 2 for the tetragonal structure; the thickness of the line represents the orbital contributions. It is evident that the Fe 3d symmetry adapted bonding bands appearing between −6 to −2 eV possess dominant As 4p character. The bands close to the Fermi level are primarily constituted by Fe 3d states with notable contributions from As 4p states. Interestingly, while the bands around Γ point forming the hole pocket in the Fermi surface possess essentially Fe 3d character, one of the electron pocket around X-point possess large As 4p character.
The detailed orbital symmetries near ε F are shown in Fig. 15. At Γ, three energy bands are observed to cross ε F -these bands are denoted by α, β and γ as shown separately in Fig. 15(g) for better clarity. These three energy bands create three hole pocket around Γ 13 . There are two electron pockets around the high symmetry point, X. It is clear from Fig. 15(c) and (d) that the γ band possesses mainly d xy character and the α & β bands have contribution from degenerate d xz and d yz orbitals 14 . These three bands possess t 2g symmetry. The e g bands constituted by − d x y 2 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 4p character except one electron pocket at X, which has As 4p 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 3d-As 4p hybridization.
In Fig. 16, we show the change in band structure in different phases along X − Γ − Z symmetry lines -in the AFM phase, the corresponding k-vector is in M − Γ − Z-direction. In collapsed tetragonal phase, the bands at Γ shifts below the Fermi level and the hole pockets around Γ disappear consistent with the ARPES data 50,51 . 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 Γ nesting with an electron pocket at X gives rise to an energy gap at ε F via folding of the energy bands at the zone boundary created by the magnetic supercell.

Discussion
From the above results, it is clear that electronic structure of CaFe 2 As 2 is very sensitive to the crystal structure and pnictogen height is an important parameter in deriving the ground state properties. The system appears to favour 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 9 . 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 CaFe 2 As 2 . To get a clear idea about this we have calculated the local magnetic moment at Fe-sites. The calculations 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 cT p1 and the cT phase with P = 0.63 GPa is called cT p2 . 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. 52 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. 38 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 [39][40][41] . The list of lattice parameters used for the calculations along with the calculated magnetic moment and ground state energies are given in Table 1.
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μ 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 bond length (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 favoured more. The converged energy of the non-magnetic phase at cT p1 structure is found to be the lowest relative to all the calculations we have done for tetragonal CaFe 2 As 2 . The converged energies for the NM and AFM phases in the orthorhombic structure is −31 meV and −298 meV, respectively exhibiting a scenario consistent with the observed ground state phase of the system in ambient conditions. From the calculated energies it is clear that a perturbation energetically favours cT phase consistent with the experimental findings.
The Fe moment appears to diminish with the application of pressure. Interestingly, the lowest moment of 1.37μ B is found for the pressure of 0.47 GPa (cT p1 phase) and then the Fe moment increases with further increase in pressure although the compression of c is more at cT p2 phase compared to cT p1 phase. The value of z As , however, is smaller in the cT p1 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 9 . 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 bond length makes the conduction electrons more itinerant and hence less local moment. It is puzzling that the Fe-As bond length is smallest in cT p1 structure and a small decrease in Fe-As bond length from 2.3405 Å in cT p2 structure to 2.338 Å in cT p1 structure changes the magnetic moment from 1.95 μ B to 1.37 μ 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 ε 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 of 2.3654 Å in the AFM ground state with orthorhombic structure is larger than that in the cT phase and shorter than the one in T phase. The Fe spin moment is found to be 1.87 μ B , which is much larger than the experimental value of 0.80 μ 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 2 . 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 [53][54][55] .
Interestingly, coupling of Ca 4s 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 4s PDOS along with Fe 3d and As 4p contributions for different structures. Ca 4s contributions in orthorhombic and tetragonal phases look almost similar with major intensities above ε F , a reasonably intense feature between −6 to −3 eV and almost no intensity in the vicinity of ε F . While the intensity above ε 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 ε 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 4p character and −2.5 to 0 eV is contributed by anti-bonding bands with dominant Fe 3d character. The intense feature between −2.5 to −4 eV in As 4p PDOS are the non-bonding As 4p states. The features in Ca 4s PDOS resemble well with structures in the As 4p PDOS as these states emerges due to the hybridization of Ca 4s states with As 4p states. Thus, the enhancement of 4s contributions below ε 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 47,48 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. In CaFe 2 As 2 , the higher degree of Ca-As hybridization in collapsed tetragonal phase and the observation of pressure induced superconductivity 35, 56-58 appears to be linked. It is to note here that while the non-hydrostatic conditions and doped systems (unit cell becomes too big for doped systems of interests) could not be treated in our calculations, the results presented here provide an understanding of the scenarios involving those issues via providing evidence that application of pressure and/or structural distortions help the system to evolve towards the proximity of quantum criticality via destruction of magnetic long range order and enhancement of the influence of charge reservoir layer in the electronic structure. We hope that these results will provide an incentive to initiate more studies in this direction.

Conclusion
In summary, we present here the detailed electronic structure of CaFe 2 As 2 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 the collapsed tetragonal phase. While the energy distribution of the density of states appears similar in tetragonal and orthorhombic phases 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 Γ 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 4s 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.

Method
Electronic structure calculations were carried out using full potential linearized augmented plane wave method (FLAPW) captured in the Wien2k software 59 . Convergence was achieved by fixing the energy convergence criteria to 0.0001 Rydberg (1 meV). For nonmagnetic calculations, we have used 10 × 10 × 10 k-points in the Brillouin zone and for Fermi surface calculations 39 × 39 × 10 k-points were used. We have used the Perdew-Burke-Ernzerhof generalized gradient approximation (GGA) 60 for our density functional theoretical calculation 61 . The Fermi surfaces were calculated using Xcrysden 62 . The collapsed tetragonal phase of CaFe 2 As 2 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) Å 3 .