Next Article in Journal
AFM and FTIR Investigation of the Effect of Water Flow on Horseradish Peroxidase
Next Article in Special Issue
Comparison of Bifurcated Halogen with Hydrogen Bonds
Previous Article in Journal
Investigation of Solubility Behavior of Canagliflozin Hydrate Crystals Combining Crystallographic and Hirshfeld Surface Calculations
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 26
Issue 2
10.3390/molecules26020304
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metal–Metal Bond in the Light of Pauling’s Rules

by
Elena Levi
1,*,
Doron Aurbach
1 and
Carlo Gatti
2,3,*
1
Department of Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
2
CNR-SCITEC Istituto di Scienze e Tecnologie Chimiche “Giulio Natta”, sezione di via Golgi, via Golgi 19, I-20133 Milano, Italy
3
Istituto Lombardo Accademia di Scienze e Lettere, via Brera 28, I-20121 Milano, Italy
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(2), 304; https://doi.org/10.3390/molecules26020304
Submission received: 22 December 2020 / Revised: 5 January 2021 / Accepted: 6 January 2021 / Published: 8 January 2021
(This article belongs to the Special Issue Chemical Bonding: A Commemorative Special Issue Honoring Professor Linus Pauling)
Browse Figures
Review Reports Versions Notes

Abstract

:
About 70 years ago, in the framework of his theory of chemical bonding, Pauling proposed an empirical correlation between the bond valences (or effective bond orders (BOs)) and the bond lengths. Till now, this simple correlation, basic in the bond valence model (BVM), is widely used in crystal chemistry, but it was considered irrelevant for metal–metal bonds. An extensive analysis of the quantum chemistry data computed in the last years confirms very well the validity of Pauling’s correlation for both localized and delocalized interactions. This paper briefly summarizes advances in the application of the BVM for compounds with TM–TM bonds (TM = transition metal) and provides further convincing examples. In particular, the BVM model allows for very simple but precise calculations of the effective BOs of the TM–TM interactions. Based on the comparison between formal and effective BOs, we can easily describe steric and electrostatic effects. A possible influence of these effects on materials stability is discussed.

1. Introduction

In 1929, based on the first structural studies, Pauling formulated five simple rules to rationalize chemical bonding in minerals [1]. One of them is a bond order (BO) conservation principle or a rule of local electroneutrality, which states that the cation charge in the crystal structure should be compensated by the negative charge of adjacent anions and vice versa. Consequently, the BO sum, Σ BOij, or the bond valence sum, BVS = Σ sij, (sum of the bond strengths according to Pauling) around atom i should be equal to its valence, Vi:
Σ BOij = Σ sij = Vi
According to Equation (1), the BO can be defined as a number of valence electrons of an i atom associated with a given bond. In 1947 Pauling related the BOij to respective interatomic distance, Rij, by empirical constants [2]. This correlation can be written as follows:
sij = exp [(R0 ij − Rij)/bij],
The constants R0 ij and bij are known as bond valence (BV) parameters, transferable for a given atom pair in different compounds. They represent an effective repulsion between i and j atoms and the softness of the bond, respectively [3,4].
Equations (1) and (2) are widely used in crystal chemistry in the framework of the bond valence model (BVM) [3,4], mainly in order to verify the crystal structure solution or to determine the oxidation state of atoms. The BV parameters for a large combination of the atom pairs were proposed [5]. For compounds with localized bonds, it was shown that, in good accordance with Equations (1) and (2), the BVSs calculated with these BV parameters are close to the expected number of valence electrons (Although some studies mark the limited predictive power of the Pauling Rules [6]). In contrast, the Pauling rules and the BVM were for a long time deemed as inapplicable to compounds with delocalized bonds, e.g., metal cluster compounds (i.e., materials with metal–metal, TM–TM (TM = transition metal), interactions) [4,7]. Only recently, the efficacy of the model for such compounds was shown [8,9,10,11,12,13,14,15,16,17].
This paper briefly summarizes the advance in the application of the BVM to compounds with metal–metal bonding, also presenting novel examples and focusing on the validity of Pauling’s rules for these materials. It is shown herein how structural and bonding peculiarities of these materials can be described and understood based on Pauling’s approach.

2. Methods: Determination of the BV Parameters for the Metal–Metal Bonds

As is clear from Equation (2), the accuracy of the BVS values depends mostly on the choice of BV parameters. A few attempts to use quantum chemistry considerations for Pauling’s rules and calculation of the BV parameters are known [18,19,20,21,22], but till now, the BVM remains empirical. Different sets of the BV parameters for the cation-anion pairs were proposed [5], but all of them are based on the direct use of Equations (1) and (2). It means that the BV parameters for a given atom pair are chosen to satisfy a local electroneutrality in a large number of compounds with this pair. As discussed below, such an approach is convenient only for solids without steric constraints. Since metal–metal bonds are commonly too short of matching their atomic surrounding in the crystal structure; it is difficult to find the TMn cluster compounds without lattice strains. In this case, the direct method cannot be used, but Equations (1) and (2) can be applied indirectly based on the valence compensation (see Equation (3) below). In this case, the BVS of the TM–TM bonds, Σ s TM–TM, is accepted to be equal to the difference between the total number of the TM electrons responsible for bonding, lTM, (e.g., six valence electrons for Mo atom in Mo6-cluster compounds) and the BVS of the TM–L bonds, Σ sTM-L. The latter can be calculated with conventional BV parameters [23].
For the determination of the BV parameters, we also used the large sets of quantum chemistry data on the effective BOs for the TM–TM bonds, which are available in the literature. In this case, the BV parameters were fitted based on the BO exponential trend vs. respective interatomic distances. The approximate values of R0 ij for the TM pairs were also calculated using empirical constants proposed by O’Keeffe and Brese in 1991 for most of the chemical elements (bij of 0.37 Å is assumed to be the same for all the i and j pairs) [24]. It was shown that all three methods (based on BO conservation, on quantum chemistry data and O’Keeffe’s constants) result in the close values of the BV parameters [16], but the first approach seems the most reliable for the moment. In particular, the evaluation of BV parameters from quantum chemical data could largely benefit from the (hopefully future) availability of a large set of BO data computed using wavefunctions of comparable quality and a single BO recipe. However, it should also be mentioned that steric and electrostatic effects in cluster compounds are commonly so pronounced that the choice of the BV parameters is less crucial as compared to the case of localized bonds.
Details of the electron counting for various cluster compounds can be found in our previous works (see, for example, Ref. [14]).

3. Results and Discussion

3.1. Validity of Equation (2) for Localized and Delocalized Bonds

To validate the BVM application to the delocalized bonds, the first question that should be discussed is “What was the reason to reject Pauling’s rules for compounds with metal–metal interactions?”. Cotton, in the book “Multiple bonds between metal atoms,” wrote [7]: “It is a general qualitative rule in chemistry that bond lengths and bond orders are inversely related… However, there is no a priori reason to expect that similar procedures will (or will not!) work in the very different realm of metal-to-metal bonds. Experience is the only test, and experience thus far has shown that M–M bonds cannot usefully be treated in such a way. We condemn as foolish and hopeless any effort to associate a unique, quantitative bond order with each and every metal–metal internuclear distance.” As an argument for this hard statement, Cotton showed that the same formal BO of the metal–metal bond might be associated with effectively different bond lengths. Indeed, e.g., the length of the Cr–Cr quaternary bond (formal BO = 4) in Cr2 complexes ranges from 1.830 to 2.541 Å [25].
The problematic point in Cotton’s argument is his correlation of the bond lengths to the formal BOs, while Pauling’s Equation (2) relates them to the effective BOs. (The formal BOs for the cation–anion bonds are those expected from the formal oxidation state of the ions. For the metal–metal bonds, the formal BOs should be evident from the structure of molecular orbitals. It is defined as half of the difference between the numbers of bonding and antibonding electrons.) In contrast to the localized bonds, the effective BOs of the metal–metal interactions can significantly differ from the formal values. As it was shown [26], the difference may reach 1 valence unit (v.u.). As an example, Figure 1 presents formal and effective BOs of the Cr–Cr bonds as a function of their lengths. The effective BOs (Mayer, Wiberg, delocalization indices, etc.) were calculated by different authors for various compounds (the references of these studies can be found in [14]). As can be seen, in contrast to the formal BOs, the calculated BOs clearly follow exponential decay, in spite of the high dispersion of the data. A similar analysis of the effective BOs for various atomic pairs obviously confirmed the validity of Equation (2) for both localized and delocalized bonds [13]. The BV parameters of the TM–TM bonds obtained by fitting the exponential curves allow for a simple but careful calculation of the effective BOs from the TM–TM distances by Equation (2) [14,15,16].

3.2. Validity of Equation (1) for Compounds with Metal–Metal Bonds

3.2.1. Stretching of Metal–Metal Bonds Evident by Valence Violations

In 1992 Brown proposed important comments to the rationalization of Equation (1), showing that this equation is valid only for the unstressed bonds [27]. Lattice strains caused by the steric mismatch between different atoms in the crystal structure result in the deviation of the BVS Σ sij, from the expected Vi value. Thus, the difference Σ sij − Vi may serve as a measure of the lattice strain and material instability [3,4,27]. A great advantage of this method is that, in combination with structural analysis, it shows the source of material instability. The problem is that, for the compounds with localized bonds, the Σ sij − Vi difference is commonly small. As a result, in practice, this method needs a very careful choice of the BV parameters, as well as experimental interatomic distance, Rij, otherwise the Σ sij − Vi difference may be caused by inaccuracy or better indeterminacy in the BVS calculations. In spite of this difficulty, the method was successfully used for perovskites [3,4,27] and electrode materials [28] to explain or even to predict their instability.
The existence of stresses in the metal–metal bonds embedded in the inorganic atomic framework was firstly marked by Schäfer and Schnering [29]. Corbett defined a steric mismatch between short metal–metal bonds and closed-shell anion surrounding as “matrix effect” and tried to use Equation (2) to describe it [30,31]. It is interesting to compare the BV parameters proposed by Corbett based on the respective distances in metals with those calculated by Levi et al. using the conservation principle, as well as numerous quantum chemistry data of effective BOs available in the literature [14] (Table 1). As can be seen, two sets of the BV parameters are relatively close to each other, but the R0 values proposed by Corbett are systematically bigger, resulting in the underestimation of the stretching in the metal–metal bonds.
Using the BVM analysis and structural approach, we can explain the bonding peculiarities of cluster compounds, e.g., the expansion of the TMn clusters with the size of surrounding anions. At first glance, such expansion is strange because the cluster size should depend only on the number of bonding (and antibonding) electrons, but not on the anions’ size. However, from structural considerations, it is clear that the higher the ligand size, the larger is the void formed by anions, and the higher is the stretching of the TM–TM bonds in the TMn cluster that occupies this void. To illustrate this statement, Figure 2a compares the BVS and the formal number of valence electrons for the Mo–Mo bonds in isostructural compounds, Cu2Mo6L14 (L = Cl, Br and I). A structure of the cluster unit of this compound is presented in the upper inset of Figure 3. An increase in the ionic radius from Cl (1.81 Å) to Br (1.96 Å) and to I (2.20 Å) results in longer Mo–Mo bonds in the Mo6 cluster, and, respectively, in lower BVS of these bonds, while the formal valence remains the same. Thus, a general conclusion can be drawn: the stretching of the metal–metal bonds is responsible for the unusually high difference between their formal and effective BOs (Figure 1 and Figure 2). The higher is the stretching; the lower are the effective metal–metal BOs (and respective BVSs).

3.2.2. Compensation of the Valence Violations for the Bonds around Transition Metals

A reasonable question arises: “Does the stretching of the metal–metal bonds decrease the material stability?” The answer is rather negative because the valence deficiency caused by stretching is compensated by the valence excess related to the compression of the metal–ligand bonds. As can be seen from Figure 2b, the effective BVSs for the Mo-L bonds in Cu2Mo6L14 is effectively larger than the formal ones, and their rise with the ionic radius of the ligands corresponds to the BVS drop for the Mo–Mo bonds. The total BVS of the Mo atoms for all three compounds are close to 6, i.e., to the number of valence electrons of the Mo atom (Note that for Mo6 clusters in these compounds, all the electrons are bonding). In the general case, in accordance with the BO conservation principle, the total BVS for the metal atom, which includes the metal–metal (TM–TM) and metal–ligand (TM–L) bonds, is close to the expected number of valence electrons, lTM, responsible of bonding:
BVSTM = Σ sTM-TM + Σ sTM-L = lTM.
Moreover, it was shown that the matrix effect is commonly associated with a more symmetric distribution of valence electrons and electron density around transition metals [10,14]. For example, in the compound Cu2Mo6I14 with an octahedral Mo6 cluster, the Mo atom is bonded to four Mo atoms and five I atoms (see the upper inset in Figure 3). The formal valence distribution is four electrons for the Mo–Mo bonds and two electrons for the Mo-I bonds. Taking into account similar lengths of all the bonds of the same type, the formal BOs should be close to 1 v.u. for the Mo–Mo bond and 0.4 v.u. for the Mo-I bond. Due to the bond strains, in real Cu2Mo6I14, the bond valences are 0.62 and 0.68 v.u., respectively. It is logical to suggest that such symmetric valence distribution caused by the matrix effect impacts the stability of the cluster units. Indeed, it is known that, in spite of the strained bonds, cluster core, Mo6Li8, remains stable even after the dissolution of cluster compounds in different solvents [32].

3.2.3. Clusters as Single Cations with Nonuniform BVS Distribution on Their Ligands

The most unusual result that follows from the BVS calculations for compounds with metal–metal bonds is that, in contrast to the TM atoms, Equation (1) is not valid for the BVS of the separate ligands, but only for a given cluster unit as a whole. Moreover, the BVS distribution on the ligands in many cluster compounds is extremely nonuniform, with a high deviation of the BVS from the expected values. As an example, Figure 3 presents cluster contribution to the valence of the ligands, BVSL = Σ sTM-L, as a function of their distance from the cluster center for the Re6-chalcogenides. (Here, BVSL is related only to the bonds between the ligand and a separate cluster.) As can be seen, this contribution to the inner ligands is about three times higher than that for the outer ligands. This result is not surprising, taking into account the structure of the cluster units: close values of the TM–L distances for the inner and outer ligands and different coordination around chalcogen atoms: three Re atoms for the inner ligands and only one for the outer ligands (see the upper inset in Figure 3). It means that three Re3+ cations contribute their valence electrons in the bonding with the inner ligand, but only one Re3+ cation assists in the bonding with the outer ligand.
It is interesting to note that the BVS distribution on the ligands can be easily understood if, in our bonding scheme, we will replace the TMn cluster with a single cation located in the cluster center (see the low inset in Figure 3). It can be shown that the curves in Figure 3 follow Equation (2), where Rij is the bond length between the imaginary TMn-cation and its ligand, and R0 ij and bij are new BV parameters. These parameters calculated for a set of the TM6-cluster compounds are presented in Table 2. A large value of bij (more than 1 Å in our case) is commonly assigned to a high difference in the electronegativity of the i and j atoms [33]. According to the R0 values, the sizes of the imaginary TMn cations are comparable with the size of such cations like Cs+ or Ba2+, but the formal charge (i.e., Σ BOij = Σ sij = Vi, Equation (1)) is much higher, ranging commonly from 11 to 17 for the Nb6, 12–16 for Mo6, 12–22 for W6 and being equal to 18 for Re6. Due to the large ionic charge, the polarizing power of the clusters should be high. Consequently, the inner ligands located closer to the cluster center should have an effectively larger charge than the outer ones.
To avoid any potential misunderstanding, it is worth recalling that the BVM notion of ionic charge does not bear any direct relation with that of atomic charge customarily adopted in theoretical and computational chemistry. The former, derived from the BVS, measures the actual “valence” of an atom as opposed to its formal value, whereas the latter is simply given by the difference between the atom’s nuclear charge and its fractional number of electrons, evaluated with some suitable recipe (Mulliken’s, Bader’s charges, etc.).

3.2.4. Charge Transfer from the Cluster to the Ligands and between the Ligands

In the previous section, we dealt with formal charges of the clusters, but, based on the anion BVSs, we can calculate their effective charges and compare them to the formal ones. For example, the oxidation state of Re in the Re6-cluster compounds is +3, and the formal charge of the Re6-cation is equal to +18, while the effective charge is +26.4 for sulfides, +27.8 for selenides and +31.1 for tellurides (average data for few compounds). For the Mo6-cluster compounds with the Mo oxidation state of +2 (the formal charge of +12), the effective charge is equal to +18.2 for chlorides, +18.5 for bromides and +20.0 for iodides. Thus, the electron redistribution around TM metals in the cluster units, associated with the matrix effect, results in the charge transfer from the cluster (and from the outer ligands) to the inner ligands. As expected from the matrix effect, the value of the charge transfer from the cluster depends on the size of the inner ligands (Figure 4a), but we can also relate it to the ligand electronegativity (Pauling scale) (Figure 4b). The higher the latter, the smaller is the charge transfer. The influence of the outer ligands on the charge transfer is much less pronounced (Figure 5). Moreover, the cluster contribution to the BVSs of the outer ligands is almost unaffected by their composition, allowing easy mutual substitution of these ligands in the synthesis.
Again, it is important to stress that analogously to the notion of “ionic charge”, that of “charge transfer” has in the BVM context a completely different meaning relative to the same term assumed in theoretical/computational chemistry. BVM’s charge transfer is related to an increase of the BVS of subsets of bonding interactions (in the present case, the TM–Li) at the expense of the other ones in a system, while in theoretical/computational chemistry, the term “charge transfer” simply expresses the (fractional) number of electrons transferred from one to another system’s moiety, as a result of the chemical interaction among these, originally isolated, moieties. As a consequence, just opposite to the BVM case, the higher the ligand electronegativity, the larger is the charge transfer from the TM cluster to the ligands in the computational/theoretical chemistry language.

3.2.5. Ligand Valence Violations as a Source of Material Instability

The difference in cations’ contribution to the valence of anions is a normal phenomenon for solids. For example, in orthorhombic NaMnO2, two crystallographically different oxygen atoms have different bonding to the Mn3+ cations, as well as different Mn input to the oxygen valence: about 0.7 and 1.1 v.u. However, this difference is compensated by the respective contribution of Na cations in the oxygen BVS. Otherwise, according to Pauling, the material should be unstable. The BV analysis of cluster compounds shows that it is much more difficult to compensate for the difference between the effective anion charge (BVS) and the formal one, especially for large clusters. Various mechanisms of such compensation for compounds with octahedral clusters are described in refs. [8]. They include (i) formation of the mixed chalcogen-halogen compositions, with the preferential occupation of the inner and outer sites in the cluster units by chalcogen and halogen atoms, respectively; (ii) occupation of the outer sites by chalcogen atoms with low oxidation state; (iii) connectivity of the cluster units by common outer ligands. In spite of this, the anion valence violations are typical for most cluster compounds.
For example, in [Bu4N]4 [Re6S8Cl6], the inner and the outer ligands are presented by sulfur and chlorine, respectively. As a result, the valence violations (Figure 6a) are less pronounced than for pure chlorides or sulfides. In Nb6I11, each outer ligand is common to two adjacent clusters, resulting in minimal valence violations for these ligands (Figure 6b). In Cs4Re6S13 different mechanisms of the valence compensation is working. The outer ligand S1 is common for two adjacent clusters, while the formal valence of the outer ligand S2 is not −2, but −1 (Figure 6c). In this compound, contributions of Cs+ cations in the BVSs of the inner and outer ligands are relatively close and do not effectively diminish the valence violations.
The clear tendency to reduce the valence violations, which appears in a large part of cluster compounds, seems to testify to a negative influence of the valence violations on material stability, but the absence of thermodynamic data on the stability of these compounds prevents the conclusion from being definitive.

3.3. Comparison with the Results of Quantum Chemistry Calculations

In this section, we try to analyze how the results of BVM application to compounds with metal–metal bonds agree with the data obtained by other methods. In Section 2, it was mentioned that we used a plethora of recent quantum chemistry data to confirm the validity of Equation (2) and to determine the BV parameters of the metal–metal bonds. Naturally, the BO values obtained by the BVM method for the TM–TM bonds agree well with quantum chemistry BO estimates. Moreover, Equation (2) is in good agreement with fundamentals of quantum chemistry, namely with exponential decay of the atomic and molecular wavefunctions (or orbitals), and thus may be qualitatively justified (see Ref. [13] for more details). Nevertheless, it is worth emphasizing that due to the high dispersion of quantum chemistry BOs (see Figure 1) obtained by using different recipes and wavefunction qualities, their exponential decay can be clearly established only in a wide range of interatomic distances and only for a sufficiently large variety of the BO data.
At first glance, it seems very simple to describe steric and electrostatic effects in cluster compounds by comparison of formal and effective BOs, regardless of the BO calculation method. However, with the exception of Corbett’s works mentioned above, these effects are known only due to the recent BVM application. The quantum chemistry studies rather compare the effective BOs obtained by different calculation methods, sometimes very complicated and time-consuming. In addition, they are focused solely on the metal–metal bonds, without analysis of respective valence violations for the metal–ligand bonds. Even in the rare cases that all the bonds were analyzed, the effects were not discussed.
As an example of such a rare case, Table 3 presents the results of the BV analysis based on the data by Baranovski and Korolkov for [Mo6S8(CN)6]6−. The bond lengths, R, and respective BOs calculated by quantum chemistry methods are taken from the original work [34]. Based on these data and geometry of the cluster unit (see the upper inset of Figure 3), we calculated the BVS separately for four Mo–Mo bonds in the Mo6 cluster, the BVS for four Mo–S bonds and a single Mo–C bond, as well as a total BVS for the Mo atom. A comparison between these BVSs and the formal values shows the same features that were discussed above: stretching of the Mo–Mo bonds and compression of the Mo–S bonds, while the total BVS of the Mo atom is close to the number of its valence electrons. The BVSs of the ligands decrease with the distance from the cluster center: they are high for the inner S ligands and low for the outer C ligands. It is worth noticing that the quantum chemistry work of Baranovski and Korolkov is not devoted to steric and electrostatic effects in cluster compounds. In spite of this, the results of this work are in excellent agreement with those obtained by the BVM application.

3.4. The BVM Application to the Electrode Materials with Metal–Metal Bonding

Compounds with metal–metal bonds often have interesting physical properties (optical, catalytic, thermoelectric, etc.). An excellent example is Chevrel phases, MxMo6L8 (L = S, Se and Te, M is various cations). In addition to superconductive and thermoelectric applications [35,36,37,38], the sulfides and selenides were used as unique cathodes in Mg batteries [39]. The latter was suggested as promising analogs of well-known Li batteries, but a slow solid-state diffusion of divalent Mg2+ cations in most of the common hosts hampered their practical use [40]. In contrast, Chevrel phases allow for fast ionic transport of various multivalent cations [41]. One of the reasons for the high ionic conductivity is the presence of the Mo6 cluster, which can easily adopt up to four electrons upon cation insertion into Chevrel phase. This is in contrast to the individual TM atoms of the common hosts, which can simultaneously adopt only one electron. However, the BV analysis revealed additional bonding peculiarity in Chevrel phases, which may affect the ionic transport.
To illustrate this peculiarity, we chose Cu-containing Chevrel phases, CuxMo6S8, with accurately determined crystal structure [35], which is very close to that of MgxMo6S8 [42]. As can be seen in Figure 7, the BVSs of two crystallographically different sulfur atoms in Mo6S8 are completely different, in spite of the fact that both S1 and S2 atoms are inner ligands located at relatively close distances from the cluster center. The reason for the BVS distinction is their different bonding to the clusters (Figure 7a). S2 is connected only to one Mo6 cluster by three Mo–S bonds, while S1 has an additional bond to the adjacent cluster, for which it serves as the outer ligand. Thus, S1 has four Mo–S bonds with similar lengths. Figure 7b shows that Cu (or Mg) insertion into Mo6S8 is associated with the BVS growth for the S2 atom. As a result, the BVS distribution on the ligands becomes more uniform.
Such change of the ionic charges (BVSs) in the anion framework upon cation insertion in the electrode material is very unusual. It is commonly accepted that the charge of the inserted cations is compensated solely by the respective change in the oxidation state of transition metal. In contrast, in the Chevrel phases, the compensation also includes anions, thus facilitating the electron transfer in the crystal structure of the host. This mechanism of the charge distribution in Chevrel phases, first discovered by the BV analysis [8], was later confirmed by theoretical and experimental studies [43,44,45].

4. Conclusions

Recent application of the valence-length correlation proposed by Pauling about 70 years ago permits a simple, fast and zero-cost calculation of the effective BOs for the metal–metal bonds. Their comparison with the formal BOs, in combination with crystal structure rationalization, gives a unique possibility to describe steric and electrostatic effects, typical for the compounds with TM–TM interactions. These effects explain the bonding and structural peculiarities of cluster compounds. They seem to be crucial for the stability of these materials.

Author Contributions

Conceptualization, E.L. and C.G.; BO calculations, E.L.; project administration, D.A.; writing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pauling, L. The principles determining the structure of complex ionic crystals. J. Am. Chem. Soc. 1929, 51, 1010–1026. [Google Scholar] [CrossRef]
  2. Pauling, L. Atomic Radii and Interatomic Distances in Metals. J. Am. Chem. Soc. 1947, 69, 542–553. [Google Scholar] [CrossRef]
  3. Brown, I.D. The Chemical Bond. In Inorganic Chemistry: The Bond. Valence Model; Oxford University Press: Oxford, UK, 2006. [Google Scholar]
  4. Brown, I.D. Recent developments in the methods and applications of the bond valence model. Chem. Rev. 2009, 109, 6858–6919. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Brown, I.D. Available online: https://www.iucr.org/resources/data/datasets/bond-valence-parameters (accessed on 3 November 2016).
  6. George, J.; Waroquiers, D.; Di Stefano, D.; Petretto, G.; Rignanese, G.-M.; Hautier, G. The Limited Predictive Power of the Pauling Rules. Angew. Chem. Int. Ed. 2020, 59, 7569–7575. [Google Scholar] [CrossRef] [Green Version]
  7. Cotton, F.A.; Murillo, C.A.; Walton, R.A. Multiple Bonds between Metal Atoms; Springer: New York, NY, USA, 2005. [Google Scholar]
  8. Levi, E.; Aurbach, D. Lattice strains in the ligand framework in the octahedral metal cluster compounds as the origin of their instability. Chem. Mater. 2011, 23, 1901–1914. [Google Scholar] [CrossRef]
  9. Levi, E.; Aurbach, D.; Isnard, O. Bond-valence model for metal cluster compounds. I. Common lattice strains. Acta Cryst. B 2013, 69, 419–425. [Google Scholar] [CrossRef]
  10. Levi, E.; Aurbach, D.; Isnard, O. Bond-valence model for metal cluster compounds. II. Matrix effect. Acta Cryst. B 2013, 69, 426–438. [Google Scholar] [CrossRef]
  11. Levi, E.; Aurbach, D.; Isnard, O. Electronic effect related to the non-uniform distribution of ionic charges in metal cluster chalcogenide-halides. Eur. J. Inorg. Chem. 2014, 2014, 3736–3746. [Google Scholar] [CrossRef]
  12. Singh, V.; Dixit, M.; Kosa, M.; Major, D.T.; Levi, E.; Aurbach, D. Is it True that the Normal Valence-Length Correlation is Irrelevant for Metal-Metal Bonds? Chem. Eur. J. 2016, 22, 5269–5276. [Google Scholar] [CrossRef]
  13. Levi, E.; Aurbach, D.; Gatti, C. Do the Basic Crystal Chemistry Principles Agree with a Plethora of Recent Quantum Chemistry Data? IUCrJ 2018, 5, 542–547. [Google Scholar] [CrossRef] [Green Version]
  14. Levi, E.; Aurbach, D.; Gatti, C. Bond Order Conservation Principle and Peculiarities of the Metal-Metal Bonding. Inorg. Chem. 2018, 57, 15550–15557. [Google Scholar] [CrossRef] [PubMed]
  15. Levi, E.; Aurbach, D.; Gatti, C. Assessing the Strength of Metal–Metal Interactions. Inorg. Chem. 2019, 58, 7466–7471. [Google Scholar] [CrossRef] [PubMed]
  16. Levi, E.; Aurbach, D.; Gatti, C. Steric and Electrostatic Effects in Compounds with Centered Clusters Quantified by Bond Order Analysis. Cryst. Growth Des. 2020, 20, 2115–2122. [Google Scholar] [CrossRef]
  17. Levi, E.; Aurbach, D.; Gatti, C. A revisit of the bond valence model makes it universal. PCCP 2020, 22, 13839–13849. [Google Scholar] [CrossRef] [PubMed]
  18. Burdett, J.K.; McLarnan, T.J. An orbital interpretation of Pauling’s rules. Am. Mineral. 1984, 69, 601–621. [Google Scholar]
  19. Mohri, F. Molecular orbital study of the bond-valence sum rule using Lewis-electron pair theory. Acta Cryst. B 2003, 59, 190–208. [Google Scholar] [CrossRef]
  20. Mohri, F. Molecular orbital study of bond–valence sum rule for hydrogen bond systems using Lewis-electron pair theory. J. Mol. Struct. Theochem. 2005, 756, 25–33. [Google Scholar] [CrossRef]
  21. Hardcastle, F.; Laffoon, S. Theoretical Justification for Bond Valence--Bond Length Empirical Correlations. JAAS 2012, 66, 87–93. [Google Scholar]
  22. Adams, S. Practical considerations in determining bond valence parameters. In Bond Valences; Springer: Berlin, Germany, 2013; pp. 91–128. [Google Scholar]
  23. Brese, N.E.; O’Keeffe, M. Bond-valence parameters for solids. Acta Cryst. B 1991, 47, 192–197. [Google Scholar] [CrossRef]
  24. O’Keefe, M.; Brese, N. Atom sizes and bond lengths in molecules and crystals. J. Am. Chem. Soc. 1991, 113, 3226–3229. [Google Scholar] [CrossRef]
  25. Wu, L.-C.; Hsu, C.-W.; Chuang, Y.-C.; Lee, G.-H.; Tsai, Y.-C.; Wang, Y. Bond characterization on a Cr–Cr quintuple bond: A combined experimental and theoretical study. J. Phys. Chem. A 2011, 115, 12602–12615. [Google Scholar] [CrossRef] [PubMed]
  26. Roos, B.O.; Borin, A.C.; Gagliardi, L. Reaching the maximum multiplicity of the covalent chemical bond. Angew. Chem. 2007, 119, 1491–1494. [Google Scholar] [CrossRef] [Green Version]
  27. Brown, I.D. Chemical and steric constraints in inorganic solids. Acta Cryst. B 1992, 48, 553–572. [Google Scholar] [CrossRef]
  28. Levi, E.; Aurbach, D. Lattice strains in the layered Mn, Ni and Co oxides as cathode materials in Li and Na batteries. Solid State Ion. 2014, 264, 54–68. [Google Scholar] [CrossRef]
  29. Schäfer, H.; Schnering, H.V. Metall-metall-bindungen bei niederen halogeniden, oxyden und oxydhalogeniden schwerer übergangsmetalle thermochemische und strukturelle prinzipien. Angew. Chem. 1964, 76, 833–849. [Google Scholar] [CrossRef]
  30. Corbett, J.D. Correlation of metal-metal bonding in halides and chalcides of the early transition elements with that in the metals. J. Solid State Chem. 1981, 37, 335–351. [Google Scholar] [CrossRef]
  31. Corbett, J.D. Chevrel phases: An analysis of their metal-metal bonding and crystal chemistry. J. Solid State Chem. 1981, 39, 56–74. [Google Scholar] [CrossRef]
  32. Prokopuk, N.; Shriver, D. The Octahedral M6Y8 and M6Y12 Clusters of Group 4 and 5 Transition Metals. In Advances in Inorganic Chemistry; Sykes, A.G., Ed.; Academic Press: San Diego, CA, USA, 1998; Volume 46, pp. 1–49. [Google Scholar]
  33. Adams, S. Relationship between bond valence and bond softness of alkali halides and chalcogenides. Acta Cryst. B 2001, 57, 278–287. [Google Scholar] [CrossRef]
  34. Baranovski, V.I.; Korolkov, D.V. Electron density distribution, polarizability and bonding in the octahedral clusters [Mo6S8(CN)6]6−, [Re6S8(CN)6]4− and Rh6(CO)16. Polyhedron 2004, 23, 1519–1526. [Google Scholar] [CrossRef]
  35. Yvon, K. Bonding and relationships between structure and physical properties in Chevrel-phase compounds MxMo6X8 (M = metal, x = S, Se, Te). In Current Topics in Materials Science; Kaldis, E., Ed.; Elsevier: Amsterdam, The Netherlands, 1979; Volume 3, pp. 53–129. [Google Scholar]
  36. Fischer, O.; Maple, M.B. (Eds.) Superconductivity in Ternary Compounds I; Springer: Berlin, Germany, 1982; Volume 32. [Google Scholar]
  37. Caillat, T.; Fleurial, J.P.; Snyder, G.J. Potential of Chevrel phases for thermoelectric applications. Solid State Sci. 1999, 1, 535–544. [Google Scholar] [CrossRef]
  38. Peña, O. Chevrel phases: Past, present and future. Physica C 2015, 514, 95–112. [Google Scholar] [CrossRef]
  39. Aurbach, D.; Lu, Z.; Schechter, A.; Gofer, Y.; Gizbar, H.; Turgeman, R.; Cohen, Y.; Moshkovich, M.; Levi, E. Prototype systems for rechargeable magnesium batteries. Nature 2000, 407, 724–727. [Google Scholar] [CrossRef] [PubMed]
  40. Levi, E.; Gofer, Y.; Aurbach, D. On the way to rechargeable Mg batteries: The challenge of new cathode materials. Chem. Mater. 2009, 22, 860–868. [Google Scholar] [CrossRef]
  41. Levi, E.; Gershinsky, G.; Aurbach, D.; Isnard, O.; Ceder, G. New insight on the unusually high ionic mobility in chevrel phases. Chem. Mater. 2009, 21, 1390–1399. [Google Scholar] [CrossRef]
  42. Levi, E.; Lancry, E.; Mitelman, A.; Aurbach, D.; Ceder, G.; Morgan, D.; Isnard, O. Phase Diagram of Mg Insertion into Chevrel Phases, MgxMo6T8 (T = S, Se). 1. Crystal Structure of the Sulfides. Chem. Mater. 2006, 18, 5492–5503. [Google Scholar] [CrossRef]
  43. Thöle, F.; Wan, L.F.; Prendergast, D. Re-examining the Chevrel phase Mo6S8 cathode for Mg intercalation from an electronic structure perspective. PCCP 2015, 17, 22548–22551. [Google Scholar] [CrossRef]
  44. Wan, L.F.; Wright, J.; Perdue, B.R.; Fister, T.T.; Kim, S.; Apblett, C.A.; Prendergast, D. Revealing electronic structure changes in Chevrel phase cathodes upon Mg insertion using X-ray absorption spectroscopy. PCCP 2016, 18, 17326–17329. [Google Scholar] [CrossRef]
  45. Yu, P.; Long, X.; Zhang, N.; Feng, X.; Fu, J.; Zheng, S.; Ren, G.; Liu, Z.; Wang, C.; Liu, X. Charge Distribution on S and Intercluster Bond Evolution in Mo6S8 during the Electrochemical Insertion of Small Cations Studied by X-ray Absorption Spectroscopy. J. Phys. Chem. Lett. 2019, 10, 1159–1166. [Google Scholar] [CrossRef] [Green Version]
Sample Availability: Samples of the compounds are not available from the authors.
Molecules 26 00304 g001 550
Figure 1. Formal and effective bond orders (BOs) of the Cr–Cr bonds vs. interatomic Cr–Cr distances.
Figure 1. Formal and effective bond orders (BOs) of the Cr–Cr bonds vs. interatomic Cr–Cr distances.
Molecules 26 00304 g001
Molecules 26 00304 g002 550
Figure 2. Effective anion charge (BVSs) and the formal numbers of valence electrons for the Mo–Mo (a) and Mo-L (b) bonds in Cu2[Mo6Li8]La6 (L = Cl, Br and I) with octahedral Mo6 clusters coordinated by Li8-cube of the inner ligands and La6-octahedron of the outer ligands.
Figure 2. Effective anion charge (BVSs) and the formal numbers of valence electrons for the Mo–Mo (a) and Mo-L (b) bonds in Cu2[Mo6Li8]La6 (L = Cl, Br and I) with octahedral Mo6 clusters coordinated by Li8-cube of the inner ligands and La6-octahedron of the outer ligands.
Molecules 26 00304 g002
Molecules 26 00304 g003 550
Figure 3. Re6 cluster contribution to the ligand valence as a function of their distance from the cluster center in cluster compounds with Re6L8-core (L = S, Se, Te). The upper inset shows the cluster unit Re6Li8La6. The rhenium and chalcogen atoms are in red and green, respectively. Re-Re bonds in octahedral clusters are marked in blue. The Re-L bonds around Re, inner and outer ligands are in red, green and pink, respectively. In the low inset, the cluster is replaced by a single cation.
Figure 3. Re6 cluster contribution to the ligand valence as a function of their distance from the cluster center in cluster compounds with Re6L8-core (L = S, Se, Te). The upper inset shows the cluster unit Re6Li8La6. The rhenium and chalcogen atoms are in red and green, respectively. Re-Re bonds in octahedral clusters are marked in blue. The Re-L bonds around Re, inner and outer ligands are in red, green and pink, respectively. In the low inset, the cluster is replaced by a single cation.
Molecules 26 00304 g003
Molecules 26 00304 g004 550
Figure 4. Charge transfer from the cluster to the ligands as a function of the size of the inner ligands (a) and the ligand electronegativity (b).
Figure 4. Charge transfer from the cluster to the ligands as a function of the size of the inner ligands (a) and the ligand electronegativity (b).
Molecules 26 00304 g004
Molecules 26 00304 g005 550
Figure 5. Charge transfer from the cluster to the ligands as a function of the electronegativity of the outer ligands.
Figure 5. Charge transfer from the cluster to the ligands as a function of the electronegativity of the outer ligands.
Molecules 26 00304 g005
Molecules 26 00304 g006 550
Figure 6. Valence violations for anions in [Bu4N]4 [Re6S8Cl6] (a), Nb6I11 (b) and Cs4Re6S13 (c). Here absolute values of the oxidation state are presented.
Figure 6. Valence violations for anions in [Bu4N]4 [Re6S8Cl6] (a), Nb6I11 (b) and Cs4Re6S13 (c). Here absolute values of the oxidation state are presented.
Molecules 26 00304 g006
Molecules 26 00304 g007 550
Figure 7. Different coordination of sulfur atoms in Chevrel phase (a) and the change of their BVSs in CuxMo6S8 (b).
Figure 7. Different coordination of sulfur atoms in Chevrel phase (a) and the change of their BVSs in CuxMo6S8 (b).
Molecules 26 00304 g007
Table
Table 1. Comparison of the bond valence (BV) parameters for the metal–metal pairs used by Corbett [30] and Levi et al. [15].
Table 1. Comparison of the bond valence (BV) parameters for the metal–metal pairs used by Corbett [30] and Levi et al. [15].
BondCorbett, 1981Levi et al, 2019
R0, Åb, ÅR0, Åb, Å
Sc-Sc2.9210.262.6950.28
Ti-Ti2.6380.262.5050.47
V-V2.4640.262.4350.43
Fe-Fe2.3670.262.260.44
Co-Co2.3230.262.110.44
Zr-Zr2.9180.262.890.43
Nb-Nb2.7080.262.640.395
Mo-Mo2.6190.262.510.34
W-W2.6350.262.5350.29
Table
Table 2. Bond valence parameters for the TM6-anion pairs (TM6 is a virtual cation located in the cluster center).
Table 2. Bond valence parameters for the TM6-anion pairs (TM6 is a virtual cation located in the cluster center).
BondR0TM6-L, ÅbTM6-L, ÅBondR0TM6-L, ÅbTM6-L, Å
Nb6-F2.961.36
Nb6-Cl3.661.07W6-Cl3.791.38
Nb6-Br3.9251.01W6-Br3.961.18
Nb6-I4.21.12W6-I4.221.08
Mo6-Cl3.751.25Re6-S4.031.09
Mo6-Br3.9051.14Re6-Se4.21.06
Mo6-I4.171.07Re6-Te4.591.14
Table
Table 3. BV analysis based on the data of Ref. [34].
Table 3. BV analysis based on the data of Ref. [34].
AtomBondsR, ÅBO, vuBVS, vuFormal Valence, vuValence Violations, vu
MoMo-Mo2.7350.6372.5483.3330.785
Mo-L 3.3812.667−0.714
All 5.92960.071
SMo-S2.5390.7312.193−2−0.193
CMo-C2.3050.4570.457−10.543
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Levi, E.; Aurbach, D.; Gatti, C. Metal–Metal Bond in the Light of Pauling’s Rules. Molecules 2021, 26, 304. https://doi.org/10.3390/molecules26020304

AMA Style

Levi E, Aurbach D, Gatti C. Metal–Metal Bond in the Light of Pauling’s Rules. Molecules. 2021; 26(2):304. https://doi.org/10.3390/molecules26020304

Chicago/Turabian Style

Levi, Elena, Doron Aurbach, and Carlo Gatti. 2021. "Metal–Metal Bond in the Light of Pauling’s Rules" Molecules 26, no. 2: 304. https://doi.org/10.3390/molecules26020304

APA Style

Levi, E., Aurbach, D., & Gatti, C. (2021). Metal–Metal Bond in the Light of Pauling’s Rules. Molecules, 26(2), 304. https://doi.org/10.3390/molecules26020304

Article Metrics

Back to TopTop