Peter C. Müllera,
Linda S. Reitz
a,
Stefan Engel
b,
Richard Dronskowski
*a and
Oliver Janka
*b
aInstitute of Inorganic Chemistry, RWTH Aachen University, Jülich-Aachen Research Alliance (JARA-CSD), RWTH Aachen University, 52056 Aachen, Germany
bInorganic Solid State Chemistry, Saarland University, Campus C4.1, 66123 Saarbrücken, Germany. E-mail: oliver.janka@uni-saarland.de
First published on 28th August 2025
Intermetallic compounds belong to an important class of materials, not only due to the sheer number of compounds known but also due to their application in everyday life. These compounds possess their very own peculiarities, especially when it comes to chemical bonding. To address this point, bonding analyses based on Crystal Orbital Bond Index (COBI) values, Löwdin charges, and – for the first time – ab initio oxidation numbers (ONai) were conducted, all extracted from delocalized plane-wave functions. From the integrated COBI values, to be understood as quantum-chemical bond orders, the differences and similarities in the bonding behavior of the elements, binary and ternary compounds in the Ca–Al–Pt system were analyzed. It became apparent that the Al–Pt interactions, almost regardless of the respective compounds, show significant covalency, while Ca–Al and Ca–Pt interactions are of ionic nature in most cases. Homoatomic Al–Al or Pt–Pt interactions, however, tend to be ambivalent, depending on the respective crystal structure of a given compound.
Disobeying the classical valence rules leads to some problems because unless the underlying wave function is known and analyzed the nature of the chemical bond cannot be understood.29–38 What are the contributions to a certain interaction between two atoms? Ionic interactions or rather covalency? Does a concept of, say, electronegativity still apply? When going to the extremes, e.g., Cs2Pt39,40 or CsAu,40–42 transparent compounds with salt-like behavior can be observed; however, when ΔEN is not as striking as in these examples, is there still a polarity in these intermetallic phases? The term ‘polar intermetallics’ is used quite frequently to describe compounds where a certain polarization within the structure is either assumed or proven. The question of how to properly address the bonding in solids dates back to almost a century: Laves43 and Pauling44–46 introduced atomic packing and “resonance” concepts for crystalline solids, while Hume-Rothery47 addressed the electronic (metallic) state and critical valence-electron concentrations.
Science has come a long way, since nowadays chemical bonding48 may be addressed computationally based on density-functional theory, either resting on the wave function49–51 or on the electron density.52–56 The main difference between both approaches lies in the examined quantity:34 a density-based method such as Bader's quantum theory of atoms in molecules (QTAIM) semi-classically partitions the electron density according to its topology. Even though this looks attractive since the density (1) is easily accessible from DFT calculations and (2) also experimentally observable from, e.g., X-ray diffraction, the density lacks essential information, which is crucial for a well-grounded evaluation. As the density relates to the absolute square of the wave function, the phase information (i.e., the sign of the wave function) is completely lost in the density, thus we cannot distinguish between bonding and antibonding, at least not from the wave function.
In a solid-state context, there is yet another hindrance: common DFT programs employ a delocalized plane-wave basis lacking local (chemical) information that would require atom-centered basis functions. In order to extract these bonding data, the program LOBSTER57 conducts a so-called projection from a plane-wave onto an atomic-orbital basis. Therefore, we regain chemistry in terms of, e.g., Löwdin charges58 and crystal orbital bond index (COBI) values, the latter including the phase information.59 Both established tools allow for a bond classification in terms of ionicity and covalency. Additionally, in this recent contribution, we also introduce wave-function derived oxidation numbers that aim to resolve the schism between quantum-mechanical charges and empirical oxidation states.
The case study in this paper will be the elements Ca, Al and Pt, selected binaries of the systems Ca–Al, Ca–Pt and Al–Pt, as well as all reported ternary compounds of the entire Ca–Al–Pt system. Fig. 1 depicts a Gibbs triangle showing all reported binary and ternary phases based on the Pearson database60 and the recent reports. Besides the crystal structures of the three elements Ca, Al and Pt,61–63 CaAl2 (ref. 64) and CaAl4,64 and CaPt2 (ref. 65) as well as Al2Pt,66 and AlPt67 were selected from the respective binary phase diagrams. Furthermore, CaAlPt, Ca2AlPt2 and CaAl2Pt were selected as representatives of ternary compounds. Our essential goal is to quantify the bonding nature in these intermetallic compounds, namely (1) to check if interatomic distances can be used as an identifier for bonding or non-bonding scenarios and (2) to identify similarities amongst the crystallographically significantly different compounds and highlight specific peculiarities that arise based on these quantum-chemical calculations.
All corresponding electronic structures, based on the optimized calculations, were projected from plane waves onto a local orbital basis using the LOBSTER (Local Orbital Basis Suite Towards Electronic-structure Reconstruction) program package.50,57,75–77 By this method, the local density-of-states matrices (i.e., energy-resolved wave-function eigenvectors) become available, which enable the calculation of the crystal orbital bond index (COBI),59 a generalized solid-state molecular bond index introduced by Wiberg78 and Mayer,79 as well as gross populations and atomic charges in a Löwdin-style formalism, in addition to ab initio oxidation numbers.58
In the present work, we do address this issue, however, by introducing a weighting factor that allocates the bonding electrons to a more electronegative atom, in the spirit of heterolytic electron partitioning. This extension of the Löwdin population analysis then yields ab initio oxidation numbers ONai that are defined as the difference between the number of electrons in the neutral atom Ne,A and the population after heterolytic bond splitting:
![]() | (1) |
As the electronegativity is not straightforwardly accessible from quantum-mechanical calculations, we derive it by referring to the Hamilton matrix elements Hμν and constructing the weighting factor wμν:
wμν = 1 + erf(α(Hνν − Hμμ)) | (2) |
The error function and the parameter α† determine the “smearing” of the electronic partitioning between two atoms, and it ensures that minute energy differences between Hμμ and Hνν of, say, 0.1 eV as well as numerical noise do not cause an unreasonably drastic electron transfer. In the present study, we chose a value of α = 10 eV−1 which, for an energy difference of 0.1 eV, leads to the assignment of 92% electron density to the more electronegative atom and the remaining 8% to the less electronegative atom. A transfer of more than 99% is then performed at a difference of 0.165 eV.
The algorithm formulated in eqn (1) and (2) is visualized in Fig. 2 for the simple example of molecular ammonia, NH3. Fig. 2a shows the calculation of formal charges that result from a homolytic splitting of bonds. As such, the nitrogen atom is left with five electrons, and each hydrogen atom keeps one electron, leading to all atoms being formally neutral. When calculating oxidation numbers (cf. Fig. 2b), all N–H bonds are split heterolytically, so the bonding electrons are assigned to the more electronegative bonding partner, nitrogen in this case. This ionic limit then leads to oxidation numbers of N−III and H+I.
Going from classical Lewis formulae to quantum chemistry, we can quantify the number of electrons located on atoms and bonds by the population analyses already introduced above. As can easily be seen in Fig. 2c, setting the weight factor α to zero, splitting the bonding electrons (indicated by solid red and blue lines) is symmetric, so in case of, say, the 1s–2px/y bond, both bonding partners receive half of the 0.59 bonding electrons, 0.29 each. In sum, this separation leads to Löwdin populations and charges of −0.91 for nitrogen and +0.30 for hydrogen.
If the weighting factor α is chosen to be larger, the partitioning of the bonding population becomes heterolytical. As the 1s orbital of hydrogen is significantly higher in energy compared to the valence orbitals of nitrogen, the electrons are shifted strictly towards nitrogen, in line with the empirical expectation from electronegativities. For the 1s–2px/y bond, the total 0.59 bonding electrons are thus shifted completely to nitrogen, mirroring the recipe of empirical oxidation numbers. In the end, the ONai derived from this algorithm (H+0.75 and N−2.25) match the empirical numbers a lot closer than the respective Löwdin charges. At this point, it should be noted that fractional ONai are an inherent feature of the population analyses used in the calculus. While integer oxidation numbers are a consequence of the classical derivation, this does not apply to quantum chemistry. Analogously, Mulliken/Löwdin charges and populations as well as bond orders by means of ICOBI have fractional values in the vast majority of cases, especially in intermetallic phases as presented in the following.
![]() | ||
Fig. 3 (a) Bond length plotted versus the ICOBI per bond. (b) Ab initio oxidation number plotted versus Löwdin charge. |
The ionic bonding analyses in this study are conducted based on Löwdin charges and, additionally, using the newly introduced ab initio oxidation numbers ONai. Their relation is shown in Fig. 3b. In all examined cases, Ca purely serves as an electron donor with charges and oxidation numbers close to the ideal +2. For Al and Pt, simple expectations are not met, however. Based on electronegativities of Al and Pt, one would assume (partially) cationic Al and anionic Pt. Fig. 3b reveals that all of the examples in our examination possess negatively charged aluminum and positively charged platinum atoms, in contrast to electronegativities of Al and Pt, respectively, 1.61 and 2.28 on Pauling's scale.80 Since Löwdin charges do not reflect the energetic orbital arrangement and symmetrically attribute the density of bonding electrons to both bonding partners, charges do not trivially turn into oxidation states. This disagreement is solved by introducing a weighting scheme based on atomic orbital energies (see the explanation above), ultimately leading to the formulation of ab initio oxidation numbers. Although these oxidation numbers are often fractional, not integer like the classical counterpart, the overall fit is significantly improved. The majority of Al atoms reside in the top left part of Fig. 3b and Pt atoms in the bottom right part. This picture can be traced back to an increased bonding electron transfer from Al to Pt that is revealed when going from Löwdin charges to oxidation numbers. As most of the systems provided in Fig. 3 follow the expected course of bond strength vs. bond length, we will focus the following more detailed discussion of our results on the most interesting members of the Ca–Al–Pt family. For reference, data on all compounds are summarized in the SI.
Based on the interatomic distances, one would expect rather weak covalent interactions to the neighbors, which can easily be verified using the ICOBI values of the nearest-neighbor bonds: Al–Al has the largest value of 0.23, Ca–Ca (0.12) and Pt–Pt (0.09) are smaller, indicating even weaker covalent bonding in these metals.81 If we sum up all of the respective bond indices, we arrive at the total bond capacity of an atom that equals its valence, similar to the empirical bond valence sum.82 Including only the 12 nearest-neighbor interactions, the valences are 2.76 for Al, 1.44 for Ca, and 1.08 for Pt.
This relation of rather small metal–metal bond orders and valences can easily be rationalized: if we consider elemental calcium, each Ca atom directly bonds with 12 nearest neighboring atoms, meaning that two valence electrons per Ca atom are distributed over all of these bonding partners. Consequently, the expected number of “shared” electrons between two Ca atoms will be around 2/12 ≈ 0.17, corresponding to a Ca–Ca bond order of ca. 0.08, as mirrored by ICOBI.
Naturally, none of these compounds show any electron transfer as all Löwdin charges are ±0. This agrees with the ab initio oxidation numbers (ONai) that all have a value of ±0, a trivial consequence of an element.
CaAl2 and CaPt2 both adopt the cubic Laves phase (MgCu2 type; Fdm).85 Here, Al4/Pt4 tetrahedra are found, which are connected via all four corners to a network. The Ca atoms reside in cavities of the said framework (Fig. 5a and b). The Al–Al distances in CaAl2 are 281 pm and therefore are well in line with elemental Al.
Despite the similar interatomic distance, the ICOBI (0.43) of this contact is larger by a factor of two, and this is directly related to the coordination number of six, as compared to CN = 12 in the fcc structure. Speaking of ionicity, one would expect the formation of aluminides, which is directly reflected in significant Löwdin charges. The Ca atoms exhibit positive charges with +1.54 for CaAl2, and Al anions are formed with a charge of −0.77. This charge transfer is even more pronounced when looking at the ONai that are +1.86 and −0.93, corresponding to Ca+II and Al−I. A recent review on binary alkaline-earth trielides utilizing the Bader formalism observed a similar charge transfer.86
In isostructural CaPt2, the Pt–Pt distances are 270 pm and therefore are significantly shorter compared to elemental Pt, suggesting covalent bonding interaction. Indeed, the same picture is evident as found in CaAl2: covalently bonded Pt–Pt and rather ionic Ca–Pt interactions. Yet, there are certain differences addressed in the following: while the formation of cationic Ca and anionic Pt is given by both Löwdin charges and ONai, the quantitative charge transfer does not match the empirical expectations. Based on electronegativity, one would assume more negative Pt compared to Al, but both the Löwdin charge and ONai are smaller than the respective values of Al in CaAl2.
From a quantum-mechanical point of view, however, elemental Al (3s2 3p1) possesses a less than half-filled valence shell, while the valence of Pt (6s2 5d8) is more than half-filled. Following the (generalized) octet rule and correlation arguments, adding an electron to Pt is energetically less favorable than adding an electron to Al. This destabilization is also visible in the (I)COBI of the Pt–Pt contacts. Comparing elemental Pt and CaPt2, the ICOBI is larger for the shorter contact in the binary compound, but only by a small amount – especially in comparison with Al/CaAl2. Pushing electrons from Ca onto Pt leads to the population of antibonding levels in the Pt–Pt interactions, thereby weakening the individual bonds.
CaAl4 adopts a monoclinic structure (CaGa4 type, C2/m) that can be derived from the tetragonal BaAl4 type structure by a group-subgroup formalism.64 Once again, the Al atoms form a network with the Ca atoms residing in cavities (Fig. 5c). Here, one would naively consider the shorter Al–Al contacts to be covalent bonding interactions, while the longer ones would be considered non-bonding. Interestingly, the rather long Al1–Al1 contacts still exhibit significant ICOBI values of about 0.23. The Al1–Al2 distances are significantly shorter, leading to an increase in the ICOBI values to slightly above the half bond order. The shortest distance, finally, is the Al2–Al2 interaction between the layered fragments with a stunningly large ICOBI value of 0.81, indicating a very strong covalent bond, approaching the classical single bond order. The Löwdin charges finally reflect the coordination environments: while Al1 exhibits in a wider sense a coordination number of 12 (Al1@Al24Al14Ca4), Al2 has a coordination number of 9 (Al2@Al2Al14Ca4) with overall shorter distances and therefore a higher overall Löwdin charge. This trend has also been observed by charge transfer analyses based on the definition by Bader.64,86 Interestingly, the ONai further differentiate between both Al sites resulting in an anionic Al2 (−1.27) and a cationic Al1 (+0.31), solid-state disproportionation, so to speak. In order to understand this behavior, a fundamental solid-state periodic property needs to be recalled, namely the Madelung field. Pure electrostatics leads to an attractive (=stabilizing) force between cations and anions. Considering the rather short Al2–Ca distance relative to Al1–Ca, it is safe to assume a stronger stabilization of Al2, in turn lowering the orbital energies of Al2, anion-like. Thus, the bonding electrons of the Al1–Al2 bonds move to the more electronegative Al2 atoms, and this ultimately leads to the formation of Al1 cations and Al2 anions, at least formally.
Finally, the compounds from the Al–Pt system should be addressed (Fig. 6). When looking at the bond indices in Al2Pt (Fig. 6a), the Al–Pt contacts are significantly shorter compared to the Al–Al interactions, even though both exhibit similar ICOBI values of 0.27 and 0.26. The ionic nature of this compound, as given by the Löwdin charges, does not reflect the empirical expectation from electronegativities. This tentative disagreement can immediately be solved by the respective oxidation numbers ONai. Using orbital energies as criteria and heterolytic charge allocation, bonding electrons of the Al–Pt bonds are transferred from Al 3p to the lower-lying Pt 5d orbitals. This way, the positively charged Pt (+0.26) adopts an anionic oxidation number (−0.82) and the negatively charged Al (−0.13) has a cationic oxidation number (+0.41). Such a charge transfer of 0.54 electrons per Al atom to Pt is only revealed by the combination of an atom-centered (ionic) and a bond-centered (covalent) analysis, as explained above.
In AlPt (Fig. 6b), a contrasting picture is observed. There are three distinct Al–Pt contacts, suggesting similar bonding based on the distances. The first two interactions indeed have identical ICOBI values of 0.25, and the third interaction, which is almost identical in length, only has an ICOBI value of 0.17. This puzzling difference shall be investigated using the energy-dependent COBI, as shown in Fig. S2. At a quick glance, we identify antibonding levels directly below the Fermi level, and their amount is largest in the weakest Al–Pt bond, so the small ICOBI is not due to decreased orbital interaction but a shift from occupied bonding to antibonding levels. Note that this destabilizing part is counteracted by the remaining Al–Pt bonds that have a multiplicity of three in contrast to the “weakest” bond appearing only once per polyhedron. On a descriptive level, we may infer that the remaining contacts are stabilized at the cost of the weak bond resulting in a net stabilizing effect. Ionicity, on the other hand, is more transparent. As also found for Al2Pt above, the Löwdin charges show counter-intuitive Al anions and Pt cations. This discrepancy in terms of electronegativities is again resolved by the oxidation numbers being significantly smaller than the Löwdin charges. This time, however, the oxidation numbers are almost neutral with ±0.05, and they do not mirror distinct ions.
In general, it can be observed that the charge transfer from Al to Pt increases with increasing Al content. This is in line with recent studies on the Al–Pt system based on QTAIM,87,88 although the absolute values reported in literature are larger than our Löwdin charges and ONai.
We start with equiatomic CaAlPt89,90 which crystallizes in the orthorhombic crystal system with the space group Pnma and adopts the TiNiSi type structure. The crystal structure can be described based on a network formed by the Al and Pt atoms with the Ca atoms residing in cavities (Fig. 7a). The Al–Pt distances are relatively short within the network; additionally, elongated Al–Al distances can be observed; however, no Pt–Pt interactions are present. When compared with the distance discussions above, at least Al–Pt bonds are present. Since the Al to Pt ratio is 1:
1, forming Al–Pt interactions is a necessity, while Al–Al and Pt–Pt interactions are not necessarily required. The Ca–Al and Ca–Pt distances are relatively short, but clearly distinct to at least the Al–Pt contacts and in the range of the sum of covalent radii (vide infra). Based on the different distances, one can interpret this structure as a polyanionic [AlPt]δ− network with Caδ+ cations residing in the cavities. This is a frequently observed picture, especially in the Al-rich compounds.
The Löwdin charges underline the picture of a polyanionic network (Ca: +1.26; Al: −1.08; Pt: −0.17); however, as already observed for the binary compounds in the Al–Pt system, i.e., the Al atoms carry a higher negative charge counterintuitive to the electronegativities. In this case, both the Al and the Pt atoms are negatively charged and also show an anionic oxidation number. The ONai, however, are in better agreement with the trend derived from the electronegativities and arrive at −1.01 for Pt while Al yields −0.65. The interatomic Al–Pt distances are in a similar range as the ones observed in the binaries, but overall higher ICOBI values are observed. This can be attributed to the additional electron transfer from the Ca atoms onto the polyanion, filling bonding levels. Interestingly, the longest Al–Pt contact exhibits the strongest covalency (0.43). Besides the Al–Pt contacts also one relatively long Al–Al distance is present. This interaction would not be considered bonding with respect to the distances found in elemental Al or the sum of the covalent radii, but the ICOBI value of 0.26 clearly indicates bonding interactions.
We now draw our attention to monoclinic Ca2AlPt2 (Fig. 7b, Ca2SiIr2 type). Here, there are chains of Pt atoms with alternating shorter and longer distances. These chains run parallel to the ab plane with an angle of 63.6° between them. The Al atoms connect two of the chains with rather short Al–Pt contacts, always bridging the longer Pt–Pt distances. The Ca atoms reside in the created cavities. This compound was reported by Doverbratt and coworkers who also analyzed the bonding.91 The authors highlight the linear platinum chains that distort pairwise into shorter dumbbells and longer Pt⋯Pt contacts. When looking at the ICOBI values (Fig. 7b), 0.30 can be found for the short Pt–Pt distance, while only 0.03 is calculated for the long distance. This is in line with Doverbratt et al.;91 however, their bond orders were 0.48 and 0.26 according to the non-quantum-chemical bond-length bond-strength approach by Brese and O'Keeffe.82 They state that “the interactions within the anionic substructures are essentially nonbonding (or slightly bonding) and indicate that the electrostatic repulsions are suppressed when replacing the anionic Ge bridging elements by cationic Al atoms”. This can be clearly seen by the ICOBI values of 0.36 for the Al–Pt interactions. The ionic bonds in Ca2AlPt2 suggest a similar decomposition as found in the previously discussed CaAlPt. Ca is cationic, and Al and Pt form an anionic network with Al being slightly more negatively charged than Pt. As found before, the oxidation numbers reveal a polarization of the Al–Pt bonds that results in more negative Pt and less negative Al.
When finally going to CaAl2Pt (orthorhombic, MgAl2Cu type, Cmcm, Fig. 8a),92 a compound with Al being the majority element is discussed, so can distinct Al–Al bonding be observed? When only focusing on the arrangement of the Al atoms, we see corrugated honeycomb layers all in boat conformation with the Pt atoms residing in the center of each hexagon. As described before, the Al–Pt contacts are the shortest ones observed in the structure, while the Al–Al distances are slightly longer. This generates [Al2Pt] layers which are separated by the Ca atoms leading to Ca–Pt distances of 314 pm and Ca–Al distances of 323 and 342 pm. The shortest distance between the layers not involving Ca atoms is Al–Al = 349 pm. From a crystal-chemical point of view, one would discard this distance as not being involved in any bonding. Interestingly, the same distance elongates to 415 pm in SrAl2Pt92 and even further to 471 pm in BaAl2Pt.92 While CaAl2Pt is quite stable towards moisture, SrAl2Pt and BaAl2Pt decompose rapidly, so does the short Al–Al distance stabilize CaAl2Pt? In the rare-earth representatives, this contact decreases down to 281 pm for ScAl2Pt,93 which is well in line with what is considered a bonding interaction based on structural considerations. In the relaxed DFT model (Fig. 8a), the distance in CaAl2Pt is slightly shorter (340 pm); however, the ICOBI value is 0.31, that is about a third of a single bond. This value is especially remarkable since the shorter Al–Al contact with a length of 292 pm has a similar ICOBI of 0.34 despite the significant bond-length difference, which addresses a considerable covalency to the longer interaction. Talking of ionicity, the Löwdin charges suggest an anionic Al–Pt network that contains Ca cations. The positive sign of the Pt charge (+0.22) turns into an anionic ONai (−0.18) that is primarily caused by a shift of bonding electrons from the Ca–Pt bonds to Pt. Al has the same values for the Löwdin charge and oxidation number.
Another interesting Al–Al interaction can be observed in Ca2Al3Pt (Fig. 8b). Here, the Al atoms form a 63 Kagome net with the Pt atoms connecting two layers. Within the Al layer, there are two different interactions, the shorter one being 272 and the longer being 285 pm. Puzzling, the shorter interaction exhibits the smaller ICOBI value of 0.37, while the significantly longer distance shows a significantly higher covalency, with 0.52 being the highest ICOBI value for an Al–Al interaction observed. The lower ICOBI for the shorter bond can be traced back to (occupied) antibonding contribution, below the Fermi energy. Note that a shorter bond length may increase orbital overlap, but this also holds for antibonding interactions that decrease the bond strength as in the present example. Why, however, does this unexpected course of bond strength vs. bond length appear in the first place for a compound with only one crystallographic site per element? To answer, we take a quick look at the crystal structure: the triangles formed by the longer Al–Al bonds are capped by a calcium atom (dashed lines in Fig. 8b), essentially forming a trigonal pyramid; this structural feature does not exist for the shorter Al–Al contacts. It is safe to assume that the cationic presence of Ca2+ stabilizes the bonds in the anionic Al-network, thus leading to an increased ICOBI despite the longer Al–Al distance. Charges and ONai of Ca2Al3Pt match our findings from the previous discussion.
Footnotes |
† We note that the entire DFT procedure leading to the orbital picture is from first principles, without empirical parameters. To arrive at the chemists’ somewhat arbitrary heterolytical bond splitting, a likewise arbitrary parameter α cannot be avoided such as to not cause unreasonable shifts of bonding electrons. |
‡ Note that we chose covalent radii here on purpose. While the metallic radii exactly match the interatomic distances by definition, we want to contextualize covalent bonding properties that are better described by covalent radii as their sum corresponds to the length of a single bond that has ICOBI = 1. |
This journal is © The Royal Society of Chemistry 2025 |