Abstract
Topological semimetals in BaAl4-type structure show many interesting behaviors, such as charge density wave (CDW) in SrAl4 and EuAl4, but not the isostructural and isovalent BaAl4, SrGa4, and BaGa4. Here using Wannier functions based on density functional theory, we calculate the susceptibility functions with millions of k-points to reach the small q-vector and study the origin and driving force behind the CDW. Our comparative study reveals that the origin of the CDW in SrAl4 and EuAl4 is the strong electron-phonon coupling interaction for the transverse acoustic mode at small q-vector along the Î-Z direction besides the maximum of the real part of the susceptibility function from the nested Fermi surfaces of the Dirac-like bands, which explains well the absence of CDW in the other closely related compounds in a good agreement with experiment. We also connect the different CDW behaviors in the Al compounds to the macroscopic elastic properties.
Similar content being viewed by others
Introduction
The interplay between charge density wave (CDW) and topological phases has received much attention due to the perspective of coupling many-body interaction with non-trivial band structure topology. For example, the nesting Weyl points gapped by chiral symmetry breaking have been proposed to give arise to topological superconductivity1,2 or CDW with the latter to realize axion electrodynamics3 for the case of opposite chirality or a topological phase with monopole harmonic order4 for the same chirality. The Fermi surface nesting (FSN) among Weyl points has been recently studied for the CDW in the quasi-1D (TaSe4)2I5,6 and the long-wavelength helical magnetic order in NdAlSi7. Very recently CDW has been observed in the chiral CoSi8,9 with Kramer-Weyl points. CDW has also been used10 to engineer topological phases to realize Dirac points close to the Fermi energy (EF). However, besides FSN, the emergence of CDW also needs the important ingredient of electron-phonon coupling11,12 (EPC). The traditional view of CDW comes from Peierlsâ instability13 with perfect FSN in 1D to induce Kohn anomaly or soft phonon mode, which results in a structural transition at low temperature. This ideal 1D picture of CDW has been pervasive in condensed matter physics, but was recently substantially revised for real materials14,15. Density functional theory16,17 (DFT) calculations18,19 and experiments20,21 on NbSe2 and TaSe2 have shown that the CDW in these compounds is determined by the real part of linear response or susceptibility function, ReÏ(q), with a large EPC interaction to soften the acoustic phonon mode, while not by the FSN from the imaginary part of susceptibility function, ImÏ(q), which gives a nesting vector not necessarily corresponding to the CDW q-vector in these specific compounds.
Here by studying the topological semimetals SrAl4, EuAl4, and the related compounds, we take a closer look at the microscopic mechanism for the CDW with topological non-trivial bands and find the critical role of EPC interaction or matrix element play for the CDW in SrAl4 and EuAl4 beyond just FSN of the nodal line Dirac-like bands. Topological semimetal compounds in the BaAl4-type structure22 have shown a range of interesting behaviors with a martensitic tetragonal-to-monoclinic structural transition in CaAl423, CDW in SrAl4 and EuAl424,25,26,27,28, and multiple magnetic transitions, topological Hall effects and Skyrmions in EuAl4 and related magnetic compounds at different temperature25,29,30,31,32,33,34,35. The CDW in SrAl4 and EuAl4 has been measured to have a small q-vector along the Î-Z direction. But there is still a lack of understanding of the origin and driving force behind this CDW. Even more interesting is among the non-magnetic isostructural and isovalent series of SrAl4, BaAl4, SrGa4 and BaGa4, a recent experiment24 has found that only SrAl4 has CDW at 243âK. Such distinctly different behaviors regarding CDW make these topological semimetals the perfect subjects to test the proposed mechanism14,15 for CDW as being dominated by EPC interaction instead of FSN.
In this paper, we report a comparative computational study to reveal the origin of the CDW observed in SrAl4 and EuAl4, but its absence in BaAl4 and the two Ga compounds. We calculate both ReÏ(q) and ImÏ(q) using Wannier functions with millions of k-points to reach the small q-vector. Although all of them are topological semimetals with nodal lines in the absence of spin-orbit coupling (SOC) and host Dirac points (DPs) with SOC as protected by the 4-fold rotational symmetry, these DPs are ~0.2âeV above the Fermi energy (EF). The dominating features at EF instead are the Fermi surface (FS) shells formed by the tip and dip of the nodal line Dirac-like band dispersion with two nearby crossings of the EF. With the imperfect FSN from Ï(q), the FS shells give the maximum of the real part of the susceptibility function, ReÏmax, along the Î-Z direction with a small finite q-vector for all the three Al compounds, but not the Ga compounds. Then among SrAl4, BaAl4, and EuAl4, the EPC calculations from density functional perturbation theory36 (DFPT) show a larger EPC interaction or matrix element in SrAl4 and EuAl4 than BaAl4 for the transverse acoustic (TA) mode at about the same q-vector of ReÏmax along the Î-Z direction to induce the TA mode to be imaginary. Our study reveals that the origin of the CDW in SrAl4 and EuAl4 is the strong EPC interaction for the TA mode at small q-vector along the Î-Z direction, besides the ReÏmax from the nested FS shells of nodal line Dirac-like bands, which explains well the absence of CDW in the other closely related isostructural and isovalent compounds, in a good agreement with experiment24. The different behaviors of the TA mode CDW can also be better understood in connection to the different macroscopic shear modulus and Poisson ratio. We find that the electron charge density redistribution between the Al network and different cation layers indicates the more ionic interaction in BaAl4 than SrAl4 and EuAl4, thus can explain these different CDW and elastic behaviors. It is interesting to find the link between the TA mode CDW with microscopic EPC interaction and the macroscopic elastic properties.
Results and discussion
Topological band structures
Using SrAl4 as an example, the body-centered tetragonal (Pearson symbol tI10) crystal structure in space group 139 (I4/mmm) is shown in Fig. 1a. The bulk first Brillouin zone (BZ) and (001) surface BZ with the high symmetry k-points are displayed in Fig. 1b. The two Al sites are non-equivalent with 4d at (0.5, 0.0, 0.25) and 4e at (0.0, 0.0, z). These Al sites form a network of buckled square lattices sandwiched between the electron-donating Sr layers at the 1a site. Such structure is the binary version of the Fe-based superconductors and the related 122 compounds, where the two non-equivalent 4d and 4e sites are replaced by 3d transition metals (Fe, Co, and Ni) and main group V elements (As, P), respectively. Noticeably, the shortest Al-Al distance of 2.598âà (DFT 2.616âà ) is between the two 4e sites from the two neighboring buckled square lattices along the 4-fold c-axis, which is slightly shorter than the Al-Al distance of 2.690âà (DFT 2.688âà ) among the same buckled square lattice (see Supplementary Table 1).
The DFT-calculated band structure of SrAl4 is plotted in Fig. 1c. The highest valence and lowest conduction bands according to simple filling are highlighted in blue and red, respectively. The multiple sawtooth Dirac-like dispersion of the valence and conduction bands near the EF point to a topological semimetal with nodal lines. Indeed, as plotted in Fig. 1e for the band gap contours on the (110) plane without SOC, there are a nodal loop around the Z point and also nodal lines away from the Î-Z direction toward the X point. More nodal lines and loops are plotted on (001) and (100) planes in Supplementary Fig. 1. These nodal lines and Dirac-like dispersion are protected by the different sets of mirror symmetries, which are similar to the topological semimetal series of ZrSiS and the related compounds with square lattices of p-orbitals37. With SOC, these crossings in nodal lines are all gapped out except for the one along the Î-Z direction as zoomed in Fig. 1d. The point group symmetry is D4h. We have used Vasp2trace38,39 to analyze the elementary band representations. From the double group representation of D4h for the spinful system, the top valence band along the Î-Z direction with the 4-fold rotation switches between two different 2-dimensioinal irreducible representations of Î7 and Î6, thus there is no mixing guarantied at the band crossing, or equivalently speaking, the pair of DPs is protected by the 4-fold rotational symmetry40,41,42,43,44,45. Our results also agree with the previous study22 on BaAl4 for the DPs along the Î-Z direction.
The locations of the DPs are also indicated in the BZ by the red dots in Fig. 1b at the momentum energy of (0, 0, ±0.1912 à â1; EFâ+â0.2186âeV). When projected on the (001) surface, the pair of DPs overlap, but there are still topological surface states converging to the DP projection as shown in Fig. 1f along the \(\bar{\varGamma }\)-\(\bar{X}\) direction. The spin texture in Fig. 1g shows that these surface states are spin-momentum locked as expected for topological non-trivial surface states. Furthermore, we have also plotted these topological surface states in the other direction of \(\bar{\varGamma }\)-\(\bar{M}\) in Fig. 1h with quite different band dispersion, but still they converge to the same DP projection. Because the DPs are 0.22âeV above the EF, the physical properties such as CDW is not directly determined by the DPs, but rather the 3D FS of the nodal line Dirac-like band dispersion near the EF.
Susceptibility functions and phonon dispersions
The 3D FS of SrAl4 is plotted in Fig. 2a. The valence band FS with hole pockets (yellow outside and blue inside) are centered around the Z point, while the conduction band FS with electron pockets (purple outside and green inside) are around the Î point. Because of the Dirac-like dispersion, when the bands cross the EF twice nearby along the same high symmetry directions with a dip or tip of the Dirac-like bands, they form FS of thin shells instead of space-filled objects. Such features can be clearly seen when the 3D FS are separated into the hole and electron pieces with a near top-view in Fig. 2b, c, respectively. To study FSN and CDW from band structure, the bare susceptibility function14, Ï(q), based on DFT single-particle Kohn-Sham bands needs to be calculated in real and imaginary parts,
where \(f\left({\varepsilon }_{n,{{{{{\boldsymbol{k}}}}}}}\right)\) and \(\delta \left({\varepsilon }_{n,{{{{{\boldsymbol{k}}}}}}}\right)\) are the Fermi-Dirac distribution and Delta functions, respectively, of εn,k the n-th band energy eigenvalue at the k point with the EF set to zero. The calculation of Ï(q) requires very dense double meshes (k and q). Using the maximally localized Wannier functions46,47 (MLWF), we are able to calculate the 3D susceptibility function Ï(q) efficiently with millions of k-points in the BZ.
The real part of the 3D susceptibility function ReÏ(q)/ReÏmax of SrAl4 is shown in Fig. 2d for the whole BZ with its maximum appearing as two small red spheres along the Î-Z direction with a small finite q-vector. As plotted along the Î-Z direction in 1D in Fig. 2g, the calculated q-vector for ReÏmax is 0.23(Ï/c) for SrAl4, agreeing well with the experimental data of 0.22(Ï/c) (dashed line) for the observed incommensurate CDW at 243 K24. Such a CDW q-vector has usually been associated with FSN, i.e., the sharp peaks in the imaginary part of susceptibility function, ImÏ(q), especially for 1D cases. But as plotted in Fig. 2e along the Î-Z direction, the peak in the ImÏ(q) of SrAl4 around 0.23(Ï/c) is not sharp, indicating an imperfect nesting due to the 3D nature of the FS in real materials. The imperfect nesting vector along the kz direction in SrAl4 arises from the pyramid-shape valence FS shell in Fig. 2b and is also indicated by an orange arrow in the 2D FS contour on the (110) plane in Fig. 2f. The calculated phonon band dispersion of SrAl4 is plotted in Fig. 2h showing an imaginary TA phonon mode (black arrow) along the Î-Z direction with a q-vector of 0.24(Ï/c), which matches that of calculated ReÏmax. This match of the q-vector of the ReÏmax and imaginary TA mode in SrAl4 shows that ReÏmax and the imperfect FSN can be an indicator for CDW. However, we will later show that imperfect FSN and the ReÏmax does not necessarily guarantee a CDW in the case of BaAl4, because BaAl4 lacks a strong enough EPC to induce a soft phonon mode.
For the other three non-magnetic isostructural and isovalent compounds of BaAl4, SrGa4, and BaGa4, their band structures and FS are plotted in Fig. 3 for comparison to SrAl4. BaAl4 has a very similar band structure and FS to SrAl4, as shown in Fig. 3a, d, only with small difference at the BZ boundaries. BaAl4 retains the pyramid-shape valence FS shell around the Z point and the interlocked pattern of the hole and electron FS shells. The ImÏ(q) of BaAl4 as plotted in Fig. 2e, also has a maximum along the Î-Z direction, which has a larger q-vector and with a more extended plateau than that of SrAl4. In contrast, the ReÏmax of BaAl4 in Fig. 2g is at a smaller q-vector than that of SrAl4. This shows a more imperfect FS nesting in BaAl4 than SrAl4 and the mismatch of the peaks in ImÏ(q) and ReÏ(q) is not unexpected, because the latter include contributions from the bands away from the EF. As plotted in Supplementary Fig. 2, the phonon band dispersion of BaAl4 does not have an imaginary phonon mode or CDW near the ReÏmax q-vector.
Next for SrGa4 in Fig. 3b, e, besides the extra band inversion at the S point, the main difference is the upshift of the valence bands along the Î-Z direction. The inner shell of the hole pockets at the Z point almost merge into a single large pocket with the next valence band touching the EF along the Z-S1 direction to form a small ring-shape pocket (red) around the Z point. As the results, the ReÏmax q-vector along the Î-Z direction for SrGa4 is reduced toward qâ=â0 as plotted in Fig. 2g. SrGa4 does not have a CDW either as shown by the absence of imaginary phonon dispersion in Supplementary Fig. 2. Moving to BaGa4 in Fig. 3c, the top valence band along the Î-Z direction is pushed to even higher energy with the next valence band now crossing the EF to give a squarish hole pocket (red) centered at the Z point as seen in Fig. 3f. The imperfect FSN along the Î-Z direction now totally disappears, resulting in no maximum of ReÏ(q) near the Î point for BaGa4 in Fig. 2g. Although the ReÏmax for BaGa4 is now at the Z point, there are no imaginary phonon modes as shown in Supplementary Fig. 2, thus no CDW either. These results from the calculated ReÏ(q) and phonon dispersion show only CDW in SrAl4 but not the other three non-magnetic compounds with the same crystal structure and number of valence electrons, which agree well with experimental observations24.
Electron-phonon coupling and CDW
Despite the similar band structure, FS and also the finite small q-vector for the ReÏmax in BaAl4 when compared to SrAl4, the absence of CDW in BaAl4 demands a closer look. To reveal the origin of the observed CDW in SrAl4 but not BaAl4, using density functional perturbative theory (DFPT), we have calculated the phonon lifetime γqv due to EPC, i.e. the imaginary part of the phonon self-energy, which is proportional to both the EPC matrix elements \({\left|g\left(q\right)\right|}^{2}\) and ImÏ(q). As plotted in Fig. 4a with different electronic smearing, the γqv for the TA mode along the Î-Z direction in SrAl4 is two to three times larger than that in BaAl4. To put the EPC of the TA mode into perspective among all the other modes, γqv can be converted into the dimensionless mode-resolved EPC strength as in \({\lambda }_{{qv}}={\gamma }_{{qv}}/\left(\pi N\left({E}_{F}\right){\omega }_{{qv}}^{2}\right)\), where N(EF) is the electronic density of states at EF. With the electronic smearing of 0.04âRy to stabilize the imaginary TA mode in SrAl4, the calculated λqv for both SrAl4 and BaAl4 are plotted on top of their phonon dispersion in Fig. 4c, d for comparison. As seen from the green shade for λqv, the largest difference is that SrAl4 (Fig. 4c) has a much larger λqv than BaAl4 (Fig. 4d) for the TA mode. Unlike SrAl4, the λqv of BaAl4 for the TA mode only has similar magnitude to the other modes near the Î and Z points. The eigenvectors of the TA mode (Eu) and three such zone-center optical modes Eg (2.3âTHz), B1g (6.9âTHz), and A1g (10.8âTHz) of SrAl4 are shown in Fig. 4e. The Eg mode corresponds to the in-plane motion of Al network, while the B1g and A1g modes at higher frequency correspond to out-of-plane motion of Al network. EPC interaction also directly affects the size of the phonon softening via the real part of phonon self-energy. The λqv decorated plots show that SrAl4 is subjected to a larger size of phonon softening than BaAl4. Given the small frequency of the TA mode near Î, the larger phonon softening in SrAl4 can induce this mode to be imaginary and give the CDW instability. As plotted in Fig. 4b for the dispersion of TA mode with different electronic smearing, in contrast to BaAl4, this mode in SrAl4 has a substantial softening with the decreasing smearing size (mimicking the electronic temperature) due to the stronger EPC interaction. At the smearing of 0.02âRy, the TA mode becomes imaginary at a q-vector around 0.24(Ï/c), in a good agreement with the experimentally observed qCDW24.
When Sr is replaced by Ba going from SrAl4 to BaAl4 and with the similar FSN, the conventional wisdom is that with a larger mass of Ba than Sr giving lower phonon frequencies, the phonon softening to imaginary mode should be easier and CDW should also exist in BaAl4, assuming EPC interaction stay almost the same. However, both these assumptions are not true for the TA mode along the Î-Z direction here. Firstly, comparing the phonon dispersion of SrAl4 (Fig. 4c) to BaAl4 (Fig. 4d), the highest phonon energy is reduced from 10.8âTHz in SrAl4 to 9.9âTHz in BaAl4 and the overall phonon energy range is rescaled by the mass factor as expected. Along the Î-Z direction, the longitudinal acoustic (LA) mode is also reduced from 3.0âTHz in SrAl4 to 2.3âTHz in BaAl4 at the Z point. However, the TA mode changes in the opposite direction and increases from 1.4âTHz in SrAl4 to 1.8âTHz in BaAl4 at the Z point as also shown in Fig. 4b. The first optical Eg mode at Î for the in-plane motion of Al network also increases from 2.3âTHz in SrAl4 to 2.7âTHz in BaAl4. Secondly, comparing the TA mode lifetime γqv between BaAl4 and SrAl4 in Fig. 4a, the lifetime of the TA mode in BaAl4 becomes smaller than that in SrAl4 for both large and small q-vectors. Note the lifetime γqv accounting for the EPC interaction has no mass dependence through Ïqν by definition. It evaluates the deformation of potential experienced by the electrons on FS due to the structural distortion from phonon modes as convoluted by the ImÏ(q). The same EPC matrix elements of the deformed potential are also used to calculate the real part of the Greenâs function for the phonon softening via convolution with ReÏ(q). In Fig. 4a, lifetime γqv shows that the EPC interaction is two to three times smaller in BaAl4 than SrAl4 for the TA mode along the Î-Z direction. Thus, when going from SrAl4 to BaAl4, both the increase in the single-particle (bare) phonon frequency and the decrease in EPC interaction of the TA mode work against the phonon softening in BaAl4. As the results, unlike SrAl4, there is no CDW in BaAl4.
The above analysis focusing on EPC interaction besides imperfect FSN can also explain the behavior of EuAl4, which has a CDW at ~140âK above the Neel temperature, albeit with an even larger mass of Eu than Ba. The band structure and FS of the non-magnetic EuAl4 (no Eu 4âfâ) in Fig. 5a, b resemble those of SrAl4 in Figs. 1c, 2a, respectively. The DP has the momentum-energy at (0.0â±â0.1928âà â1; EFâ+â0.1508âeV). The FS shell of the EuAl4 conduction band is similar to that of SrAl4, except for the larger electron pocket at the N point than SrAl4. The FS shell of the EuAl4 valence band is also similar to that of SrAl4, giving a similar imperfect FSN in Ï(q) as plotted in Fig. 5d, e with a slightly smaller nesting q-vector at 0.19(Ï/c) than SrAl4. Noticeably, the ReÏmax of EuAl4 has a relatively narrow peak similar to SrAl4 at the small q-vector, rather than the extended plateau of BaAl4. The 3D ReÏ(q)/ReÏmax of EuAl4 plotted in the full first BZ in Fig. 5c shows the maximum is indeed only along the Î-Z direction. For the TA mode dispersion plotted in Fig. 5e in comparison to SrAl4 and BaAl4, although in term of atomic mass Sr < Ba < Eu, the TA mode of BaAl4 is the highest among the three in the order of SrAl4â~âEuAl4â<âBaAl4 at the Z point. Then although they all have similar ReÏmax q-vector along the Î-Z direction in Fig. 5f, the three EPC interactions of the TA mode are quite different. As plotted in Fig. 5g, the γqv of EuAl4 has an interesting behavior as a function of q when compared to SrAl4 and BaAl4. At small-q (<0.2), the γqv of EuAl4 increases fast like SrAl4, but then plateau like BaAl4 and next decreases at large-q to be even smaller than BaAl4. Although the maximum of γqv for the TA mode in EuAl4 is smaller than that of BaAl4, its value at small-q (<0.2) instead is larger than BaAl4 and closer to that of SrAl4. The small q-vector range of the EPC interaction is the most relevant to CDW because of the small ReÏmax q-vector. Thus, the EPC strength shown by γqv at small-q (<0.2) is in the order of SrAl4â>âEuAl4â>âBaAl4, explaining why both SrAl4 and EuAl4 have CDW, while BaAl4 does not, even though all three have imperfect FSN and ReÏmax q-vector.
The TA mode CDW here involves a local shear distortion perpendicular to the c-axis. It is very interesting to notice that among the calculated bulk elastic properties (see Table 1), the bulk modulus (B) of 50.6âGPa for BaAl4 is slightly smaller than the 52.9âGPa for SrAl4, reflecting a larger cation size of Ba than Sr, giving both larger a and c lattice constants with a slightly smaller c/a ratio. However, the shear modulus (G) of 36.8âGPa for BaAl4 is larger than the 29.1âGPa for SrAl4. This corresponds to a much smaller Poisson ratio of 0.207 for BaAl4 than the 0.268 for SrAl4, which means a compression along the c-axis has a less in-plane expansion in response for BaAl4 than SrAl4, i.e., the in-plane interaction is stiffer for BaAl4 than SrAl4. Then moving to EuAl4 with the smallest lattice constants among the three, both B of 57.1âGPa and G of 33.9âGPa increase comparing to SrAl4. But the G of EuAl4 is still smaller than that of BaAl4, resulting in a Poisson ratio of 0.252 similar to that of SrAl4, not BaAl4. Thus, the in-plane interaction in EuAl4 is still softer than BaAl4.
To better understand these differences from electronic structure, we have plotted and compared the electron charge density difference of \(\rho \left({{XAl}}_{4}\right)-\rho \left(X\right)-\rho \left({{Al}}_{4}\right)\) for Xâ=âSr, Ba and Eu, respectively in Fig. 5hâj. The charge density redistributions between the Al network and the different cation layers show that there is more electron transferred from the Ba layer to Al network and also more charge accumulation (yellow) at the boundary between the Ba and Al network than the cases of Sr and Eu. The more ionic character of the interaction in BaAl4 with more in-plane charge accumulation makes it harder for the in-plane shear distortion between the Al network and Ba layer, which explains a much smaller Poisson ratio for BaAl4 than SrAl4 and EuAl4. This also means the TA mode softening for CDW with the local shear distortion in BaAl4 is much harder than that in SrAl4 and EuAl4. It is interesting to find the connection between the CDW with microscopic EPC interaction and the macroscopic elastic properties, because the EPC matrix element calculation involves all the bands like the elastic property calculation, in contrast to the susceptibility function calculation involving only the bands near EF.
Conclusions
In conclusion, using the approaches based on density functional theory, we have explained the origin of charge density wave (CDW) in SrAl4 and EuAl4 in comparison to the BaAl4, SrGa4, and BaGa4. Although all of them have nodal lines in the absence of spin-orbit coupling (SOC) and become Dirac semimetals with SOC, the Dirac points are ~0.2âeV above the Fermi energy (EF). The Dirac-like dispersion of the valence and conduction bands form Fermi surface (FS) shells providing imperfect FS nesting (FSN) along the Î-Z direction for the three Al compounds, but not the Ga compounds due to the second valence band being pushed up to cross the EF. The susceptibility functions, Ï(q), have been calculated efficiently on a dense mesh of millions of k-points with the maximally localized Wannier functions (MLWF). The imperfect FSN is verified by the maximum of the real part of susceptibility function, ReÏmax, along the Î-Z direction with a small q-vector for the Al compounds, while it is absent in the Ga ones. Then among the three Al compounds, the electron-phonon coupling (EPC) calculations show a large EPC interaction at small q-vector for the transverse acoustic (TA) mode along the Î-Z direction for SrAl4 and EuAl4, which provides the driving force for the CDW to soften the TA mode, while it is absence in BaAl4. Our study reveals that the origin of the CDW in SrAl4 and EuAl4 is the strong EPC interaction for the TA mode along the Î-Z direction at small q-vector, besides the ReÏmax from the nested FS shells of nodal line Dirac-like bands, which explains well the experimental observations. We also connect the different TA mode CDW distortion to the different macroscopic shear modulus and Poisson ratio. We find that the electron charge density redistribution between the Al network and different cation layers indicates the more ionic interaction in BaAl4 than SrAl4 and EuAl4, thus explains these different CDW and elastic behaviors.
Methods
Density functional theory16,17 (DFT) calculations have been performed with PBE48 exchange-correlation functional including SOC using a plane-wave basis set and projector augmented wave method49, as implemented in the Vienna Ab-initio Simulation Package50,51 (VASP). We use a kinetic energy cutoff of 300âeV, Î-centered Monkhorst-Pack52 (11âÃâ11âÃâ11) k-mesh, and a Gaussian smearing of 0.05âeV. The ionic positions and unit cell vectors are fully relaxed with the remaining absolute force on each atom being less than 1âÃâ10â4 eV/à . The DFT-calculated lattice parameters with PBEâ+âSOC are within ~2% of the experimental data, showing good agreement (see Supplementary Table 1). Phonon band dispersions have been calculated with finite difference method using PHONOPY53 on a (3âÃâ3âÃâ4) supercell of the conventional cell with 360 atoms on the k-mesh of (4âÃâ4âÃâ2). To calculate susceptibility functions efficiently, maximally localized Wannier functions (MLWF)46,47 and the tight-binding model have been constructed to reproduce closely the band structure within ±1âeV of the Fermi level by using Group II sd and Group III sp orbitals. The susceptibility functions in Eqs. (1, 2) have been calculated with four bands, two valence and two conduction bands, around the EF on the dense (120âÃâ120âÃâ90) k and q-mesh using the MLWFs, where the Fermi-Dirac distribution is sampled at the temperature of 100âK and the Delta functions are approximated with Gaussian functions with a smearing of 0.02âeV. The surface spectral functions have been calculated with the surface Greenâs function methods54,55 as implemented in WannierTools56. Electron-phonon coupling (EPC) interaction has been calculated with density functional perturbative theory (DFPT) as implemented in Quantum Espresso57 (QE) using ultra-soft pseudopotentials with a kinetic energy cutoff of 50âRy, a (3âÃâ3âÃâ3) q-mesh and (9âÃâ9âÃâ9) k-mesh with different smearing. In Supplementary Fig. 3 we show the convergence of calculated phonon band dispersions of SrAl4 with respect to increased k-mesh to (18âÃâ18âÃâ18) and q-mesh to (4âÃâ4âÃâ4) in QE using DFPT and also a larger supercell of (4âÃâ4âÃâ4) of the conventional cell with 640 atoms and increased k-mesh in VASP using the finite displacement method with PHONOPY. The second-order elastic properties have been calculated using the stress-strain methodology58,59. The crystal structures and electron charge density difference have been visualized with VESTA60.
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
References
Qi, X.-L., Hughes, T. L. & Zhang, S.-C. Topological invariants for the Fermi surface of a time-reversal-invariant superconductor. Phys. Rev. B 81, 134508 (2010).
Li, Y. & Haldane, F. D. M. Topological nodal cooper pairing in doped weyl metals. Phys. Rev. Lett. 120, 067003 (2018).
Wang, Z. & Zhang, S.-C. Chiral anomaly, charge density waves, and axion strings from Weyl semimetals. Phys. Rev. B 87, 161107 (2013).
Bobrow, E., Sun, C. & Li, Y. Monopole charge density wave states in Weyl semimetals. Phys. Rev. Res. 2, 012078 (2020).
Gooth, J. et al. Axionic charge-density wave in the Weyl semimetal (TaSe4)2I. Nature 575, 315â319 (2019).
Shi, W. et al. A charge-density-wave topological semimetal. Nat. Phys. 17, 381â387 (2021).
Gaudet, J. et al. Weyl-mediated helical magnetism in NdAlSi. Nat. Mater. 20, 1650â1656 (2021).
Li, G. et al. Chirality locking charge density waves in a chiral crystal. Nat. Commun. 13, 2914 (2022).
Rao, Z. et al. Charge instability of topological Fermi arcs in chiral crystal CoSi. Sci. Bull. 68. 165â172 (2024).
Lei, S. et al. Band engineering of dirac semimetals using charge density waves. Adv. Mater. 33, 2101591 (2021).
Grimvall, G. The Electron-Phonon Interaction in Metals. (North-Holland, Amsterdam, 1981).
Giustino, F. Electron-phonon interactions from first principles. Rev. Mod. Phys. 89, 015003 (2017).
Peierls, R. E. Quantum Theory of Solids. (Oxford University Press, New York, 1955).
Johannes, M. D. & Mazin, I. I. Fermi surface nesting and the origin of charge density waves in metals. Phys. Rev. B 77, 165135 (2008).
Zhu, X., Cao, Y., Zhang, J., Plummer, E. W. & Guo, J. Classification of charge density waves based on their nature. Proc. Natl Acad. Sci. USA 112, 2367â2371 (2015).
Hohenberg, P. & Kohn, W. Inhomogeneous electron gas. Phys. Rev. 136, B864âB871 (1964).
Kohn, W. & Sham, L. J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 140, A1133âA1138 (1965).
Calandra, M., Mazin, I. I. & Mauri, F. Effect of dimensionality on the charge-density wave in few-layer 2H-NbSe2. Phys. Rev. B 80, 241108 (2009).
Johannes, M. D., Mazin, I. I. & Howells, C. A. Fermi-surface nesting and the origin of the charge-density wave in NbSe2. Phys. Rev. B 73, 205102 (2006).
Weber, F. et al. Extended phonon collapse and the origin of the charge-density wave in 2H-NbSe2. Phys. Rev. Lett. 107, 107403 (2011).
Valla, T. et al. Quasiparticle spectra, charge-density waves, superconductivity, and electron-phonon coupling in 2H-NbSe2. Phys. Rev. Lett. 92, 086401 (2004).
Wang, K. et al. Crystalline symmetry-protected non-trivial topology in prototype compound BaAl4. Npj Quantum Mater. 6, 28 (2021).
Miller, G. J., Li, F. & Franzen, H. F. The structural phase transition in calcium-aluminum compound (CaAl4): a concerted application of Landau theory and energy band theory. J. Am. Chem. Soc. 115, 3739â3745 (1993).
Nakamura, A. et al. Characteristic Fermi surfaces and charge density wave in SrAl4 and related compounds with the BaAl4-type tetragonal structure. J. Alloy Compd. 654, 290â299 (2016).
Araki, S. et al. Charge density wave transition in EuAl4. J. Phys. Soc. Jpn 83, 015001 (2014).
Shimomura, S. et al. Lattice modulation and structural phase transition in the antiferromagnet EuAl4. J. Phys. Soc. Jpn 88, 014602 (2019).
Kaneko, K. et al. Charge-density-wave order and multiple magnetic transitions in divalent europium compound EuAl4. J. Phys. Soc. Jpn 90, 064704 (2021).
Ramakrishnan, S. et al. Orthorhombic charge density wave on the tetragonal lattice of EuAl4. IUCrJ 9, 378â385 (2022).
Nakamura, A. et al. Transport and magnetic properties of EuAl4 and EuGa4. J. Phys. Soc. Jpn 84, 124711 (2015).
Takagi, R. et al. Square and rhombic lattices of magnetic skyrmions in a centrosymmetric binary compound. Nat. Commun. 13, 1472 (2022).
Moya, J. M. et al. Incommensurate magnetic orders and topological Hall effect in the square-net centrosymmetric EuGa2Al2 system. Phys. Rev. Mater. 6, 074201 (2022).
Meier, W. R. et al. Thermodynamic insights into the intricate magnetic phase diagram of EuAl4. Phys. Rev. B 106, 094421 (2022).
Shang, T. et al. Anomalous Hall resistivity and possible topological Hall effect in the EuAl4 antiferromagnet. Phys. Rev. B 103, L020405 (2021).
Zhu, X. Y. et al. Spin order and fluctuations in the EuAl4 and EuGa4 topological antiferromagnets: A muSR study. Phys. Rev. B 105, 014423 (2022).
Gen, M. et al. Rhombic skyrmion lattice coupled with orthorhombic structural distortion in EuAl4. Phys. Rev. B 107, L020410 (2023).
Baroni, S., de Gironcoli, S., Dal Corso, A. & Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 73, 515â562 (2001).
Schoop, L. M. et al. Dirac cone protected by non-symmorphic symmetry and three-dimensional Dirac line node in ZrSiS. Nat. Commun. 7, 11696 (2016).
Vergniory, M. G. et al. A complete catalogue of high-quality topological materials. Nature 566, 480â485 (2019).
Vergniory, M. G. et al. All topological bands of all nonmagnetic stoichiometric materials. Science 376, eabg9094 (2022).
Yang, B. J. & Nagaosa, N. Classification of stable three-dimensional Dirac semimetals with nontrivial topology. Nat. Commun. 5, 4898 (2014).
Po, H. C., Vishwanath, A. & Watanabe, H. Symmetry-based indicators of band topology in the 230 space groups (vol 8, 50, 2017). Nat. Commun. 8, 50 (2017).
Slager, R. J., Mesaros, A., Juricic, V. & Zaanen, J. The space group classification of topological band-insulators. Nat. Phys. 9, 98â102 (2013).
Kruthoff, J., de Boer, J., van Wezel, J., Kane, C. L. & Slager, R.-J. Topological classification of crystalline insulators through band structure combinatorics. Phys. Rev. X 7, 041069 (2017).
Song, Z. D., Zhang, T. T., Fang, Z. & Fang, C. Quantitative mappings between symmetry and topology in solids. Nat. Commun. 9, 3530 (2018).
Bradlyn, B. et al. Topological quantum chemistry. Nature 547, 298â305 (2017).
Marzari, N. & Vanderbilt, D. Maximally localized generalized Wannier functions for composite energy bands. Phys. Rev. B 56, 12847â12865 (1997).
Souza, I., Marzari, N. & Vanderbilt, D. Maximally localized Wannier functions for entangled energy bands. Phys. Rev. B 65, 035109 (2001).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865â3868 (1996).
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953â17979 (1994).
Kresse, G. & Furthmuller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comp. Mater. Sci. 6, 15â50 (1996).
Kresse, G. & Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169â11186 (1996).
Monkhorst, H. J. & Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 13, 5188â5192 (1976).
Togo, A. & Tanaka, I. First principles phonon calculations in materials science. Scr. Mater. 108, 1â5 (2015).
Sancho, M. P. L., Sancho, J. M. L. & Rubio, J. Quick iterative scheme for the calculation of transfer-matrices - application to Mo(100). J. Phys. F. Met. Phys. 14, 1205â1215 (1984).
Sancho, M. P. L., Sancho, J. M. L. & Rubio, J. Highly convergent schemes for the calculation of bulk and surface green-functions. J. Phys. F. Met Phys. 15, 851â858 (1985).
Wu, Q., Zhang, S., Song, H.-F., Troyer, M. & Soluyanov, A. A. WannierTools: an open-source software package for novel topological materials. Comput. Phys. Commun. 224, 405â416 (2018).
Giannozzi, P. et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys.-Condens Mat. 21, 395502 (2009).
de Jong, M. et al. Charting the complete elastic properties of inorganic crystalline compounds. Sci. Data 2, 150009 (2015).
Golesorkhtabar, R., Pavone, P., Spitaler, J., Puschnig, P. & Draxl, C. ElaStic: a tool for calculating second-order elastic constants from first principles. Comput. Phys. Commun. 184, 1861â1873 (2013).
Momma, K. & Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 44, 1272â1276 (2011).
Acknowledgements
The susceptibility function calculations in this work at Ames National Laboratory were supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division. Topological band structure analysis was supported by the Center for the Advancement of Topological Semimetals, an Energy Frontier Research Center funded by the U.S. Department of Energy Office of Science, Office of Basic Energy Sciences through the Ames National Laboratory under its Contract No. DE-AC02-07CH11358. Electron-phonon coupling calculations were supported by the Ames National Laboratory LDRD. The Ames National Laboratory is operated for the U.S. Department of Energy by Iowa State University under Contract No. DE-AC02-07CH11358.
Author information
Authors and Affiliations
Contributions
P.C.C. and L.-L.W. conceived and designed the work. L.-L.W. designed and performed the ab initio calculations on susceptibility functions and topological band structure analysis. L.-L.W. and N.K.N. performed electron-phonon coupling and elastic calculations. All authors discussed the results and contributed to the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Communications Physics thanks Mina Yoon, Rob Moore, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wang, LL., Nepal, N.K. & Canfield, P.C. Origin of charge density wave in topological semimetals SrAl4 and EuAl4. Commun Phys 7, 111 (2024). https://doi.org/10.1038/s42005-024-01600-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s42005-024-01600-1