Superlattice Symmetries Reveal Electronic Topological Transition in CaC6 with Pressure

1. Introduction

Graphite intercalant compounds (GICs) are among the many interesting carbon (C)-containing families of superconductors [1,2,3]. Superconductivity was discovered in the compound calcium intercalated graphite (CaC₆) in 2005 exhibiting a critical temperature (T_c) of 11.5K at ambient pressure [4]. In this compound, calcium (Ca) is intercalated between graphene layers. Initially, there was ambiguity regarding the number of graphene layers and Ca stacking in the primitive unit cells, as CaC₆ was first believed to adopt an hexagonal structure with double layers [5], similar to other metal-intercalated graphite compounds like YbC₆ (which has a T_c of 6.5K) [4]. This misconception was clarified through the synthesis of bulk CaC₆ [6], which confirmed a distinctive rhombohedral symmetry ($R\overline{3}m$) compared to other GICs. This discovery of T_c in CaC₆, coupled with aspiration to broaden the applications of graphene, ignited substantial research interest that yielded several pivotal findings [7,8,9].

Understanding of superconductivity in CaC₆ has been the subject of ongoing debate among researchers. In 2005, Csányi et al. investigated the influence of the free-electron-like interlayer states on the superconductivity of GICs, including YbC₆ and CaC₆, using density functional theory (DFT) [10]. They hypothesised that a weak coupling between the interlayer states and the graphene layers in GICs creates conditions favourable for soft charge fluctuations. Such conditions could, in turn, foster s-wave superconductivity via an excitonic pairing mechanism. Mazin [5], using DFT calculations, postulated a different origin for superconductivity in YbC₆ and CaC₆. Mazin proposed a different mechanism involving vibrations from the intercalating elements with the dominant role of the intercalant electronic states at the Fermi level (FL), leads to a pronounced coupling with soft intercalant modes. However, it should be highlighted that both Csányi and Mazin overlooked the correct symmetry of CaC₆ in their analyses.

Weller in 2005 [4] stressed the absence of a direct correlation between the amount of charge transferred and T_c. He speculated that the superconductivity evident in GICs could be associated with a resonant valence bond mechanism. A thorough first-principles electron-phonon coupling study was carried out on CaC₆ by Calandra et al. [11]. Employing the McMillan formula [12], the critical temperature was computed by Calandra et al. [11], aligning well with experimental data. Notably, the calculated isotope effect for Ca stood at 0.24K. This value contrasts with experimental observation of a large isotope effect, with about 0.5K shift in T_c for Ca reported by Hinks et al. [13].

Additionally, Hinks et al. [13] highlighted the differences in charge transfer between the two compounds CaC₆ and MgB₂. In MgB₂, complete charge transfer from magnesium to the boron network results in strong coupling between carriers and high-energy in-plane boron vibrations, leading to high-T_c superconductivity. However, in CaC₆, charge transfer is incomplete, leading to the presence of an interlayer band at the FL. This interlayer band couples with calcium phonons, and this interaction is believed to be responsible for the relatively high T_c observed in CaC₆. The precise details of this interaction and its contribution to the observed T_c are still not fully understood. The discrepancy for the Ca isotope effect suggests that the electron-phonon mechanism alone might not fully account for the superconducting attributes of CaC₆ [13].

Subsequent experimental investigations have yielded further evidence elucidating the superconductivity mechanism in CaC₆. A landmark study by G. Lamura [14] marked the first measurement of the magnetic penetration depth in bulk CaC₆, revealing a depth of 720$\pm 80Å$. Lamura emphasised that a single gap analysis aptly captures the data with minimal determination error. He attributed the pairing to the interaction between Ca s-band electrons and the phonon modes of both C out-of-plane and Ca in-plane. Further validation of the s-wave superconducting gap came from scanning tunnelling spectroscopy [15]. Both DFT calculations [16] and directional point-contact spectroscopy [17] identified an anisotropic gap, noting values of 1.35meV and 1.71meV, respectively. The inconsistency among experimental observations, as well as between theoretical predictions and experimental findings concerning the superconducting properties of CaC₆, highlights the necessity for further in-depth investigation.

Like isotopic effects, understanding pressure effects on the superconducting properties of CaC₆ can also provide insights into underlying mechanisms. Smith et al. [18] investigated the pressure dependence of T_c for CaC₆ and YbC₆ up to 1.2 GPa. CaC₆ exhibited a linear T_c increase with pressure, while YbC₆’s T_c initially rose, peaked, and then declined. In a subsequent experiment [19], the T_c of superconducting CaC₆ displayed a prominent linear rise with pressure, reaching 15.1K at 7.5 GPa. However, at 8 GPa, a shift to a new phase was noted, characterised by a diminished T_c and inferior metallic properties at room temperature. This phenomenon is believed to be due to pressure-induced phonon softening linked to an in-plane Ca phonon mode. Separately, Gauzzi et al. [20] analysed the room-temperature crystal structure of bulk CaC₆ under pressures up to 13 GPa. At 9 GPa, they identified an order-disorder transition coupled with lattice softening. Contrary to expectations of symmetry reduction, the CaC₆ structure favoured a disordered arrangement of intercalant Ca atoms in the a–b plane. While the $R\overline{3}m$ space group symmetry remained unchanged, there was a significant rise in isothermal compressibility and a pronounced broadening of Bragg peaks.

The Fermi surface (FS)-dependent superconducting gap and electron-phonon coupling in CaC₆ were studied by Sugawara et al. [21] and Valla and Pan [22] using Angle Resolved Photoemission (ARPES). They discovered that the superconducting gap is anisotropic and varies depending on the FS and provided insights into the electron-phonon coupling in the compound. Yang et al. [23] provided further insight into the superconducting mechanism of CaC₆ using ARPES, emphasizing the critical role of the interaction between the π* antibonding orbitals and interlayer bands. Their analysis of superconducting gaps and electron-phonon coupling strengths highlighted the importance of these parameters in achieving the superconducting phase transition. Additionally, they suggested the possibility of inducing superconductivity in a monolayer of graphene by creating an adatom superlattice. Theoretical investigations also extended beyond the superconductivity of CaC₆. Rahnejat et al. [24] observed a charge density wave (CDW) in CaC₆ using ARPES, which did not cause any distortion of the carbon nuclei, highlighting the remarkable rigidity of graphene. This finding suggested the possibility of inducing CDW and even superconductivity in a graphene-based field-effect transistor through electron-doping.

Geometrical and topological aspects of superconductivity, both in real and reciprocal spaces, are currently topical and gaining widespread interest from the research community. Topology helps to understand not only the details of the physical phenomena but also some general regularity connecting physical behaviours [25]. Topological properties are relevant for quantum materials, which include superconductors, graphene, topological insulators, Weyl semimetals, quantum spin liquids, and spintronic devices, among others [26]. Topology contributes to the scientific search for universality, assisting in identification of underlying organizing principles, and evaluation of concepts that can be developed to enable a deeper understanding [27].

In this work, using first-principles calculations, we systematically investigate the changes in topology of electronic band structures (EBSs) and FSs of CaC₆ under various pressures. We evaluate two software versions of DFT and benchmark specific computational parameters in order to determine, at meV-scale, details of the EBS and FS variation with pressure. The pressure-driven alterations to EBSs, FSs and electronic charge are analysed in terms of a double superlattice cell along the c axis using a hexagonal lattice with Space Group P1. The information derived in this article by using superlattices further validates their introduction and is complementary to the information provided in a companion article [28].

In CaC₆, bands at or near the FL show a cosine-modulated dependence. This dependence infers a bonding/antibonding (or in-phase/out-of-phase) modulation as noted in earlier work [28,29]. Interpretation of results from the superlattice perspective identifies an electronic topological transition (ETT), consistent with well-established understanding of ETT’s of other layered and non-layered superconductors. The findings provide fresh insight into the principles driving superconductivity in CaC₆, thereby enriching our understanding of superconducting properties.

2. Methods

We have conducted a comprehensive DFT analysis on the EBSs and FSs of CaC₆ within a pressure range from 0 GPa up to 16 GPa using Quantum ESPRESSO [30] and Materials Studio CASTEP [31] for comparison. Crystal structures are visualised with Crystal Maker V11.0.2 using geometry optimised cell parameters of experimentally determined values.

Due to the small superconducting gap for CaC₆ (~1.7 meV), achieving convergence that accurately predicts ground state properties is essential. Previous studies have indicated the feasibility of deducing critical temperatures from electronic band structures under higher convergence conditions [32,33,34]. Considering the necessity for high precision in determining the superconducting gap, to an accuracy of about ${10}^{-1}$ meV, we implemented a meticulous benchmarking process, with detailed results provided in the supporting file. This process entailed a thorough evaluation of parameters including plane wave cut-off energies, pseudopotentials, and k-point grids [32]. We provide details of this evaluation in the Supplemental section.

Our benchmarking tests varied cut-off energies from 40 Ry (=544 eV) to 200Ry (=2720 eV) and have shown that low cut off energies lead to large deviations from consistent results. As shown in Figure S2 (Supplemental section), we observed that the total energy of CaC₆ converges within a margin of 1 meV for cut-off energies higher than 120 Ry. Additionally, fluctuations in Fermi energies (FEs), with a consistent magnitude of ~ 0.2 meV, were only noted when cut-off energies exceeded 100 Ry. Given that Δ, the zero-temperature superconducting gap, is 1.79±0.08 meV [17,21], selecting cut-off energies above 120 Ry was deemed appropriate for accurate representation.

Further, we evaluated the total system energies and Fermi surface of CaC₆ using a variety of k-point grids. The convergence behaviour for both total energy and FE under different k-point grids is detailed in Figure S3 (Supplemental). Our findings suggest that the total energy stabilises using k-point grids of density higher than a 24×24×24 sampling grid, with energy fluctuations ranging from -1 meV to 1 meV beyond this point. The FE demonstrated similar trends, achieving convergence with the 24×24×24 k-point grid.

To select suitable pseudopotentials, we adopted a systematic approach, testing various options and comparing preliminary results with experimental data. In the context of CaC₆, a rhombohedral crystal structure has been experimentally determined [6]. The experimental X-ray diffraction data provide lattice constants: $a=5.17\mathrm{Å}$and α = ${49.55}^{\mathrm{o}}$. These lattice parameters are converted to an equivalent centrosymmetric hexagonal cell with D_3d Schoenflies point group symmetry for EBS calculations.

To assess and optimise the choice of pseudopotentials, we performed unit cell relaxation using ten distinct pseudopotentials. The pseudopotentials include GGA+PAW, GGA+US, LDA+PAW, LDA+US, GGA+US+rVV10 [35], GGA+US+vdW-DF [36,37], GGA+US+ vdW-DF2 [38], GGA+US+vdW-DF3-opt1 [39], GGA+US+vdW-DF3-opt2 [39] and GGA+US+vdW-DF-C6 [40]. These pseudopotentials encompass a range of exchange-correlation approximations and treatments of electron-electron interactions including van der Waal (vdW) interactions [41]. We evaluate these pseudopotentials in order to comprehensively investigate their influence on calculated properties.

By employing this comprehensive array of pseudopotentials and drawing comparisons with experimentally determined lattice constants and other pertinent properties, we aim to identify the pseudopotential that best aligns with observed behaviour of CaC₆, ultimately enhancing the accuracy and reliability of our computational predictions. Figure S4 (Supplemental) indicates the disparities between simulated lattice parameters and corresponding experimental data, for various pseudopotentials noted above. The chart offers insight into the accuracy of different pseudopotentials in reproducing experimental lattice parameters (assumed to be measured at low temperature, although the temperature conditions of the experimental XRD analysis are not disclosed [6]) and follows earlier investigations on pseudopotential choices for MB₂- and MB₆-type compounds [42,43,44].

As shown in Figure S4 (Supplemental), simulations using GGA+PAW, GGA+US and GGA+US+rVV10 pseudopotentials yield lattice parameter values that most closely approximate experimental values. While all three options display promising alignment with experimental data, considerations of computational efficiency became pivotal in the selection process. For simulations employing the GGA+US pseudopotential, the discrepancies in lattice parameters compared to experimental data stand at 0.19% for $a$ and 0.46% for $\alpha .$ These small deviations underscore the reliability of the chosen pseudopotential and capacity to represent the properties of CaC₆.

Following these benchmarking calculations, we chose to employ ultrasoft pseudopotentials [45] with a generalised gradient approximation (GGA) [46,47] for the exchange-correlation functional, based on the closest match between calculated and experimental reported lattice parameters. Eigenfunctions were expanded using a plane-wave basis set with a cut-off energy of 120 Ry. For the wavefunctions and density of states, we used a $24\times 24\times 24$ Monkhorst-Pack grid [48] and a 0.02 Ry Methfessel-Vanderbilt smearing [49]. We also explored norm-conserving pseudopotentials within the LDA and GGA approximations using CASTEP with cut-off energy 990 eV and Dk-grid 0.005 Å^-1 for additional comparisons.

For CaC₆, based on the atomic position symmetry, it is generally accepted that a rhombohedral structure applies with lattice parameters $a_R=0.517nm$and $\alpha ={49.55}^o$ and group symmetry $R\overline{3}m$ [28]. This structure can also be represented by an equivalent hexagonal structure with lattice parameters $a=b=0.433nm,\alpha =\beta ={90}^o,and\gamma ={120}^o.$ In a rhombohedral view of CaC₆, the calcium atoms are stacked along the c-axis in an AαAβAγ sequence. The AA stacking pattern represents two adjacent layers of graphene. In this structure, the intercalant metal atoms within each layer occupy one of three prismatic hexagonal sites, labelled as α, β, and γ. Following standard conventions [50] for high symmetry paths as depicted in Figure S5, key symmetry directions ΓZ and ΓA are related as follows (see also ref [28]):

ΓZ^R = 3ΓA^H

(1)

ΓZ^R = p/c.

(2)

c =13.572A

(3)

(a₁* + a₂* + a₃*)/3 = p/c

(4)

Where a₁*, a₂* and a₃* represent the reciprocal unit cell vectors in the rhombohedral lattice, while c denotes the hexagonal lattice spacing in the direction perpendicular to the graphene layers.

DFT calculations considered both a primitive rhombohedral unit cell and the equivalent hexagonal counterpart. The specific lattice constants used as input for calculations were taken from X-ray diffraction results of bulk CaC₆ [6]. Figure 1 shows schematics of the hexagonal structure, after geometry optimisation, for selected isostatic pressures. We have also explored several supercell calculations, which could lead to folding of reciprocal space [28,29] in selected high-symmetry directions. Figure 2 shows a schematic of the 2c double hexagonal supercell. The supercell periods effectively track the bonding/antibonding modulations (or alternating phase changes) in reciprocal space, indicated by an electronic band with cosine function near the FL [51].

Using an hexagonal cell with equivalent symmetry to the rhombohedral cell, the original cosine function, which runs along the direction perpendicular to the planes or layers, has a periodicity of the 2c-lattice parameter in real space [28]. By using a 2c-supercell, the periodicity of a 1c-lattice folded at the mid-point is equivalent in reciprocal space and defines a new BZ boundary where the non-bonding condition is defined.

For the EBS of the supercell, two cosine branches exist, one with lower energy for the bonding condition and another with higher energy for the antibonding condition. These two branches intersect the Γ point at the maximum bonding (lowest) and maximum antibonding (highest) energy values. The periodicities of these phase modulations may also originate from, and correspond to, atomic orbital symmetry as noted in earlier work on MgB₂ [32].

For clearer and deeper insight, we mapped the FSs under varying pressures, using FermiSurfer [52] to visualize the FSs from the Quantum ESPRESSO output data. Charge transfer data were calculated using the software Bader Charge Analysis [53,54,55,56]. CASTEP results on FSs were exported in bitmap format from their display in the software with a limited selection of FS branches and used for complementary illustration of information.

3. Results and Discussion

This section presents our findings from the DFT analysis of CaC₆, divided into three parts. We begin by looking at how the EBSs and FSs fold when double supercells are introduced. Then, we explore how EBSs and FSs change with different pressure levels and identify the ETT. Finally, we discuss the connection between the ETT, or Lifshitz transition, to the superconducting properties of CaC₆. Through this approach, we shed light onto factors that determine the superconducting behaviour of CaC₆ and identify further common features with ETTs for other superconductors.is section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

3.1. Electronic Band Structures, Fermi Surfaces and Superlattices

The DFT calculated EBS and reciprocal path directions the for the 2c double hexagonal unit cell are shown in Figure 3a,b, respectively. Equivalent results for the rhombohedral unit cell with high symmetry paths which follow the conventions described in reference [50] are given in Supplementary Figures S5 and S6, for convenience.

A nearly free electron band in the $\mathit{\Gamma }\mathit{Z}$ reciprocal space direction, with a folded cosine-modulated shape for the 2c double hexagonal cell, can be observed crossing the FL (see band and labelled regions inside the yellow box in Figure 3a). The mid-point of the cosine-modulated band is an inflection point at Z for the derivative of the band [29,57], separating folded regions with bonding and antibonding character, from Γ to $\mathit{Z}$ and back from $\mathit{Z}$to Γ, respectively. Constructing a double superlattice along 2c in real space folds the reciprocal space in two at the $\mathit{Z}$ point and folds the EBS; thus, highlighting clearly the bonding and antibonding regions [29,57].

For the hexagonal cell shown in Figures 1a and 1b, the EBS is triply folded in the c*-axis direction (see reference [28]) as expected from the relationship between reciprocal space directions for rhombohedral and hexagonal unit cells given in Equations (1)–(4). A double hexagonal supercell has the correct 2c-lattice dimension that locates the required folding at the midpoint of the original cosine function for distinct separation of bonding and antibonding regions. From Supplemental Figure S6, which displays the projected energy bands and DOS for the primitive rhombohedral unit cell, we can see that the majority of electrons at the Fermi level have Ca 4s and C with 2p_z-orbital character. There is a small fractional contribution of electrons from C with 2s-orbital character, particularly for the cosine-shaped band.

Figure 4a,b show cross sections of the FS of CaC₆ calculated for the 2c double hexagonal supercell at 0 GPa using CASTEP. In Figure 4a an additional FS band compared to Figure 4b is shown with partial transparency. However, as the number of bands increases, it becomes more difficult to distinguish overlapping sections. Figure 4b shows a magnified view of the FS shown in Figure 4a, but with one less FS band (notice that the circle arc in the equatorial, horizontal direction is missing in Figure 4b). On the other hand, the circular caps at the top and bottom of the FS, which were barely, if at all, visible in Figure 4a, are more clearly distinguishable in Figure 4b. This is one of the advantages of QE over CASTEP for displaying multiple FSs, because cross sections can be precisely cut at chosen planes, while CASTEP shows overlapping volumetric contours.

Figure 4c shows the procedure of how the rhombohedral FS surface can be folded six times to obtain the 2c double hexagonal FS. Notice that for the FSs directly calculated for the 2c double superlattices, the interesctions of the FSs with the BZ boundaries are perpendicular (see red circled regions in Figure 4b) and that adjacent FS branches at potential intersections develop gaps instead. These gaps are more clearly noticeable in CASTEP calculations because of the use of much coarser grid (0.005 Å^-1) compared to that used in QE (0.0015 Å^-1).

3.2. Electronic Band Structures and Fermi Surfaces with Pressure

The folded nearly free interlayer band crossing the FL (predominantly with Ca 4s-orbital character) shifts in position towards higher energy values with increasing pressure, as shown in Figure 5. The equivalent calculation results for rhombohedral unit cells are displayed in supplemental Figure S8. Figures 5 and S8 demonstrate that the highest critical temperature (T_c) for CaC₆ is achieved when the interlayer band is half-filled, after which there is a notable decline as the band increases in energy. Similarly, the charge transfer from calcium to carbon in CaC₆ exhibits a comparable pattern, increasing with pressure, peaking around 7.5 GPa, and then decreasing thereafter (see Table S1).

Figure 6 shows volumetric representations and cross sections of the FS along the $\mathit{\Gamma }-\mathit{Z}-\mathit{L}$ plane as a function of pressure for the rhombohedral primitive unit cell (see also Figures 4, S7 and S8). The colour bar represents the relative value of the Fermi velocity of the electrons. For the rhombohedral unit cell, the apex of the FS of predominant Ca 4s-orbital character (i.e., the top of the slightly distorted sphere) is clearly observed moving towards $\frac{\mathit{\Gamma }\mathit{Z}}{2}$ with increased pressure (see Figure 6e). The apex reaches $\frac{\mathit{\Gamma }\mathit{Z}}{2}$ at about 7.5 GPa and continues below $\frac{\mathit{\Gamma }\mathit{Z}}{2}$ at pressures equal to and above 8 GPa, which matches the evolution of the Fermi level (E = 0) relative to the midpoint energy $E_{\frac{\Gamma Z}{2}}$ of the cosine band in the EBS’s shown in Figures 5 and S8.

3.3. Bonding/Antibonding, Superlattice and Electronic Topological Transitions

Figure 7a–d show the FSs as a function of pressure for the double supercell, constructed by the folding procedure of the rhombohedral reciprocal unit cell as described in Figure 4c. Figure 7e,f (enlarged view of 7e) show results obtained from calculations at 4 GPa using the double hexagonal supercell (with CASTEP software). A stack of two adjacent reciprocal unit cells is shown in Figure 7e.

As shown earlier in Figure 4 for 0 GPa, gaps are also generated in the double hexagonal supercell calculation for 4 GPa, because of the requirement that bands must intersect BZ boundaries perpendicularly in order to guarantee continuity of the bands and their derivatives. These FS gap regions relate to the establishment of superconducting gaps, simultaneously with the establishment of key nesting relationships [28,29,57], which link the EBS, the FS and special phonon vectors from phonon dispersions (PDs). A detailed discussion of such nesting phenomena is provided in a separate, companion article [28].

In Figure S10 (Supplement), at approximately 7.5 GPa, there is a notable change in the experimentally determined superconducting transition temperature (T_c) for CaC₆ (see also Figure 5 of reference [28]). This pressure coincides with a significant reduction in the diameter of the almost spherical Ca 4s-orbital FS to the point where the pockets (labelled 0 in Figure 7a,b) just disappear, and pockets (labelled 1 in Figure 7a–d) become disconnected at the BZ boundary. These alterations in the topology of the FS and the EBS are key to understanding the electronic transport properties of CaC₆ as a function of pressure.

3.3.1. Open and Closed Fermi Surfaces—Topological Transition with Pressure

Open and closed are two types of topologically different FSs [25]. Only closed cross-sections can appear for closed FS, causing closed trajectories for the movement of electrons in the presence of a magnetic field. In contrast, open FSs can generate both open and closed cross sections. Figure 8 displays periodically repeated FSs for CaC₆ calculated for the 2c double superlattice.

As demonstrated in Figure 7a–f, the FS, characterized by a nearly pure 4s orbital nature, extends to the adjacent BZ until pressure reaches 7.5 GPa, at which point a singularity is observed at the connection point. At 16 GPa, the FS with a pure 4s character is confined entirely within the BZ of the double hexagonal supercell. Below 7.5 GPa, the folded FS produces necking of spherical Fermi regions where continuous domains along the c*-direction can be identified (see Figure 7a–f).

As the pressure is increased, these spherical regions become disconnected, breaking up what was previously (at lower pressure) an open, interconnected FS path (Figure 7a,b) and transitioning into a closed FS loop (Figure 7c,d; see also Figure 8) [25,58]. We posit that at about 7.5 to 8 GPa, the coupled superlattice delineated nesting relationship that we have proposed [28] moves from the two internal open FSs to two open FSs, once the inner FS becomes disconnected and closed as shown in Figure 8.

3.3.2. Isoenergetic Fermi Surfaces—Electron Dynamics and Magnetic Fields

The isoenergetic FSs for the 2c double hexagonal supercell of CaC₆ at 0 GPa is displayed in Figure 9. These isoenergetic surfaces are closely linked to the dispersion of the electronic bands [25] and play a significant role in understanding electron dynamics. They can be used to determine the effective mass under a magnetic field [25]. Since experimental determination of superconductivity always requires exposure of the sample to an external electromagnetic field, either electric for a resistive transition determination, or magnetic for a Meissner effect measurement, or electromagnetic via exposure to electromagnetic radiation, at the local level, a range of isoenergetic FSs are always involved.

The topological changes described above are often referred to as Electronic Topological Transitions (ETT) [25,59]. This is a concept which largely evolved from earlier ideas introduced by Lifshitz [25,60]. This concept is typically associated with superconductivity as a function of pressure, as emphasised by Bianconi and colleagues [61,62,63,64,65,66]. Thus, continuity inside and across BZ boundaries of interconnected open FS regions for free-electron-like coupled electron transport appears to be a requirement for superconductivity. When these FS regions are open [67], they appear to correlate with the presence of superconductivity. When they are closed there may be a connection to charge density waves [24].

3.4. Fermi Velocity and Connection to the Superconducting Gap

Such FS topological transitions described above, involving shifts from open to closed FSs, or vice versa, particularly under pressure, display strong resemblance to Lifshitz transitions [68]. These are often associated with changes of T_c and/or even the destruction of superconductivity [61,62,63,64,65,66]. Figure S10 (Supplemental) is a graph of the experimentally determined superconducting transition temperature (T_c) for CaC₆, reproduced and adapted from reference [19]. A sharp drop in T_c takes place at about 8 GPa.

The following fundamental relationship is known from the BCS model for superconductivity, where the order parameter or gap energy $∆\left(0\right)$ is related to the Fermi velocity $v_F$, the reduced Planck’s constant $\mathrm{\hslash }$ and the coherence length $\xi \left(0\right)$[69] as follows:

$$∆\left(0\right)=\frac{{\mathrm{\hslash }v}_F}{\pi \xi \left(0\right)}$$

(5)

Different sections of the FSs display clear differences in Fermi velocities $v_F$; therefore, a k-dependent gap energy $∆\left(k\right)$ anisotropy is expected and the coherence lengths must be correspondingly different. From the colour scale of Figure 6, the Fermi velocity in red areas is about twice the Fermi velocity in green areas. Therefore, coherence length anisotropy is expected, and in agreement with experimental results for CaC₆.

Emery et al. [6] have observed different coherence lengths in the ab and c directions, ${\xi }_{ab}$ and ${\xi }_c$, of 35 nm and 13 nm, respectively. Since the gradient of the FS determines the main direction of movement of the FS under an external field [25,70], we expect that for the Ca 4s dominated spherical FS, a shorter coherence length in the c-direction is accessible. For the predominantly C 2p_z dominated FS, the determined coherence length is predominantly in the ab-direction. Therefore, the green and red FS regions (in terms of Fermi velocities) would have, on average, gaps with a ratio of:

$$\frac{{(v}_F/2)/13}{\raisebox{1ex}{$v_F$}\!\left/ \!\raisebox{-1ex}{$35$}\right.}=\frac{35}{26}=1.35$$

(6)

Therefore, a change in T_c is expected related with the Fano-Feshbach resonance of the gap between different areas of the Fermi surface if, for different reasons, the anisotropy influences a specific measurement, or if the dominant FS that determines properties alternates between green and red Fermi velocity regions. Cross-sectional views of the $\mathsf{\Gamma }-\mathbf{X}-\mathbf{X}1$ plane as a function of pressure is displayed in Figure S11 (Supplemental). Figure S11 shows that the Fermi surface exhibits the most hybridised characteristics at a pressure of 7.5 GPa. This hybridization is likely crucial as it facilitates the intersection and convolution of the FSs for carbon and calcium, forming what is referred to as intercalant FSs [11,21]. According to Sugawara et al. [21], these intercalant FSs and resultant interlayer bands are essential for the stabilization of a superconducting state in CaC₆.

Based on the discussion involving Fermi velocities and coherence lengths above and the discussion on ETT with open or closed FS pockets in section 3.4, the gap energies can be different in different regions of the FS; in particular, at the Lifshitz transion for “appearing” at a new Fermi surface area and at “opening a neck” topological transitions. The critical temperature show minima and maxima at the anti-resonance suppression and resonance amplification,respectively, for pair transfer between different Fermi surface areas [61,62,66,71,72]. Thus, explaining the transition temperature dependence with pressure for CaC₆ near the ETT controlled by the superstructure.

4. Conclusions

A comprehensive analysis of the electronic properties for CaC₆ under varying pressures has been presented. Using DFT calculations, we have investigated the evolution of topologies for the Fermi surface and electronic band structures of CaC₆. The Fermi surface of CaC₆ consists of three bands, with two of these bands intersecting and gradually decoupling with increasing pressure. This evolution of the FS is concurrent with a sharp drop in superconducting properties, indicating a strong correlation between FS and EBS topology and superconductivity in CaC₆.

Our analysis reveals that the highest T_c for superconductivity in CaC₆ occurs when the cosine-shaped interlayer band is half-filled. Folding the FS and EBS at the mid-point of this characteristic cosine-shaped band by introducing a 2c supercell (producing a halving of the rhombohedral c* direction), the FS of the double superlattice defines interconnected pockets, or open loops, which are favourable for nearly free electron coupled movement and superconductivity when the pressure is below 7.5 GPa. At pressures above 8 GPa, the original FS pockets become disconnected and transform into closed loops, manifesting substantial changes in superconducting properties.

By invoking a superlattice perspective of material properties, motivated by the bonding/antibonding character encapsulated in the cosine-modulated electronic band, we not only clarify the bonding/antibonding interaction, but also uncover a Lifshitz or electronic topological transition with significant implications for electronic properties with change in external conditions such as pressure. This analysis brings the pressure behaviour and characteristics of CaC₆ in line with those of many other layered and non-layered superconductors, based on key fundamental principles and mechanisms that influence superconducting properties_.

Our study enhances understanding of the pressure-dependent electronic properties of CaC₆ and offers insight into the mechanisms that control the critical temperature for superconductivity in this compound. The findings not only advance our knowledge of superconductivity in graphite intercalation compounds by identifying important EBS and FS characteristics as requirements for superconductivity, but also opens potential new paths for optimising the anisotropic superconducting gap in graphite intercalated compounds and for exploring novel materials with improved superconducting properties under pressure.