# Activation of C–H and B–H bonds through agostic bonding: an ELF/QTAIM insight

Emilie-Laure Zins, Bernard Silvi, Mohammad Esmail Alikhani

► To cite this version:

Emilie-Laure Zins, Bernard Silvi, Mohammad Esmail Alikhani. Activation of C–H and B–H bonds through agostic bonding: an ELF/QTAIM insight. Physical Chemistry Chemical Physics, 2015, 17 (14), pp.9258-9281. 10.1039/C4CP05728G. hal-01141959v2

HAL Id: hal-01141959

https://hal.sorbonne-universite.fr/hal-01141959v2

Submitted on 1 Mar 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

CrossMark

← click for updates

Cite this: Phys. Chem. Chem. Phys.,
2015, 17, 9258

# Activation of C–H and B–H bonds through agostic bonding: an ELF/QTAIM insight†

Emilie-Laure Zins,<sup>ab</sup> Bernard Silvi<sup>cd</sup> and M. Esmaïl Alikhani*<sup>ab</sup>

Agostic bonding is of paramount importance in C–H bond activation processes. The reactivity of the σ C–H bond thus activated will depend on the nature of the metallic center, the nature of the ligand involved in the interaction and co-ligands, as well as on geometric parameters. Because of their importance in organometallic chemistry, a qualitative classification of agostic bonding could be very much helpful. Herein we propose descriptors of the agostic character of bonding based on the electron localization function (ELF) and Quantum Theory of Atoms in Molecules (QTAIM) topological analysis. A set of 31 metallic complexes taken, or derived, from the literature was chosen to illustrate our methodology. First, some criteria should prove that an interaction between a metallic center and a σ X–H bond can indeed be described as “agostic” bonding. Then, the contribution of the metallic center in the protonated agostic basin, in the ELF topological description, may be used to evaluate the agostic character of bonding. A σ X–H bond is in agostic interaction with a metal center when the protonated X–H basin is a trisynaptic basin with a metal contribution strictly larger than the numerical uncertainty, i.e. 0.01 e. In addition, it was shown that the weakening of the electron density at the X–H<sub>agostic</sub> bond critical point with respect to that of X–H<sub>free</sub> well correlates with the lengthening of the agostic X–H bond distance as well as with the shift of the vibrational frequency associated with the ν<sub>X–H</sub> stretching mode. Furthermore, the use of a normalized parameter that takes into account the total population of the protonated basin, allows the comparison of the agostic character of bonding involved in different complexes.

Received 8th December 2014,
Accepted 23rd February 2015

DOI: 10.1039/c4cp05728g

www.rsc.org/pccp

## I. Introduction

In the context of catalytic processes, understanding the molecular mechanism involving transition metals is of paramount importance. Indeed, since the pioneering studies of Chatt and Davidson¹ on one hand, and Bergman² and Jones³ on the other hand, a considerable number of investigations have demonstrated that transition metal complexes may be involved in C–H bond activation. A thorough study of structures and reactivity of alkyl transition metal complexes has led to the introduction of a new concept: the term “agostic” was proposed to characterize the formation of a 3 center–2 electron (3c–2e) interaction

leading to the activation of a C–H bond around transition metal centers.⁴⁻⁷ The choice of a specific term to underline the importance of the activation of the C–H bond proved to be remarkably appropriate. Indeed, thorough investigation of various catalytic processes further lead to the development of numerous reactions involving agostic interactions, such as the alkane oxidative addition,⁸ C–H bond elimination,⁹ transcyclometallation,¹⁰ cyclometallation of benzoquinoline,¹¹ and Ziegler–Natta polymerization.¹²,¹³ In processes such as cyclometallation using d-block transition metals, it was proposed that agostic interactions may play a decisive role. Indeed, these interactions may either favor or prevent some cyclometallations.¹⁴ A mechanistic study has shown how agostic bonding may lead to deactivation processes of the Grubb catalyst.¹⁵ Thus, a thorough understanding of the agostic bonding may help in designing new reagents that may react in a specific way.

Despite its paramount importance in catalytic processes, the identification of an agostic interaction is far from obvious. Indeed, even if some approaches were proposed to experimentally or theoretically characterize these 3c–2e interactions, a consensual tool that may qualitatively describe the strength of every kind of agostic bonding is still missing.

The formation of an agostic bonding comes from an interaction between a C–H σ bond and an unoccupied orbital of a

<figure><img src="image_3.png" /></figure>

<footnote>ª Sorbonne Universités, UPMC Univ. Paris 06, MONARIS, UMR 8233, Université Pierre et Marie Curie, 4 Place Jussieu, case courrier 49, F-75252 Paris Cedex 05, France. E-mail: esmail.alikhani@upmc.fr</footnote>

<footnote>ᵇ CNRS, MONARIS, UMR 8233, Université Pierre et Marie Curie, 4 Place Jussieu, case courrier 49, F-75252 Paris Cedex 05, France</footnote>

<footnote>ᶜ Sorbonne Universités, UPMC Univ. Paris 06, Laboratoire de Chimie Théorique (LCT), UMR 7616, Université Pierre et Marie Curie, 4 place Jussieu, case courrier 137, F-75252 Paris Cedex 05, France</footnote>

<footnote>ᵈ CNRS, Laboratoire de Chimie Théorique (LCT), UMR7616, Université Pierre et Marie Curie, 4 place Jussieu, case courrier 137, F-75252 Paris Cedex 05, France</footnote>

<footnote>† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4cp05728g</footnote>
hypovalent transition metal center. Criteria were proposed to determine whether an organometallic complex contains an agostic interaction. More specifically, the following geometric parameters were established: a C-H agostic bonding should be characterized by a distance between the metallic center and the hydrogen atom in the 1.8-2.3 Å range, as well as an M-H-C angle in the 90-140° range.⁷ On the other hand, the presence of sterically constrained ligands or pincer ligands, for instance, may lead to profound geometric distortions.¹⁶ Furthermore, close to the threshold values, it may be difficult to conclude whether a complex contains an agostic interaction or not. Thus, these geometric criteria are not unambiguous, without mentioning the difficulty in determining them in some cases, specifically in dynamic systems.

As far as other “weak” interactions are concerned, criteria derived from topological studies were proposed to classify hydrogen bonding into three categories, and to propose a quantitative scale for these types of interactions.¹⁷ Since there are some common points between agostic and hydrogen bonding,¹⁶ the use of a similar approach to characterize agostic bonding comes naturally.

Several reviews aimed at presenting experimental and theoretical tools to identify agostic bonding.¹⁸⁻²⁰ In the present article, we will first briefly recall the different definitions that were/are considered for agostic bonding. Their theoretical studies are then briefly summarized, with a particular emphasis on topological approaches. The following part of this article is dedicated to the presentation of a new methodological approach aiming at estimating the agostic character of bonding. Representative examples of several kinds of agostic bonding taken from the recent literature on the topic are used to illustrate this approach. The topological description of different types of M···H-X, interactions with X = B or C, will be presented. A classification of agostic bonding based on the strength of the interaction is then proposed. To this end, the use of statistical descriptors to qualitatively evaluate the strength of an agostic bonding is validated by comparison with experimental parameters.

## II. Agostic bonding: several definitions and nuances

Despite the fundamental importance of the concept of agostic bonding in terms of reactivity and chemical activation of σ bonds, the definition of these interactions is neither unique nor unambiguous. The main definitions associated with these interactions and the context in which they were introduced are presented herein.

### 1. Historical definition and evolution

The “agostic” term was initially introduced to describe a specific intramolecular 3c-2e interaction in which a σ C-H bond “could act as a ligand to a transition metal center”.⁷ Thus, a specific term was coined for interactions involving the C-H bond to underline their importance in catalytic processes.

<figure><img src="image_4.png" /><figcaption>Fig. 1 Types of agostic bonds that will be investigated in the present study.</figcaption></figure>

First of all, β agostic interactions were characterized and thoroughly investigated both theoretically and experimentally. α C-H agostic bonding appears to be more rare. Even if in the initial definition, the term “agostic” was specifically devoted to the C-H bond, intramolecular 3c-2e interactions between a metallic center and a σ bond involving H bonded to a heteroatom are now almost always called “agostic”. The agostic bonds that will be discussed in the present article are presented in Fig. 1.

From an experimental point of view, the characterization of an agostic bond is mainly based on four criteria:¹⁸

* crystallographic data,

* NMR chemical shifts to high field $δ = -5$ to $-15$ ppm,

* reduced NMR coupling constants $J = 75$ to 100 Hz,

* low vibrational frequencies $ν_{C-H} = 2700-2300$ cm⁻¹.

As far as geometric properties are concerned, Brookhart and Green proposed in their seminal article,⁷ the criteria summarized in Table 1 may be used to distinguish between an agostic and an anagostic interaction:

* Based on QTAIM calculations, five topological criteria were proposed to characterize agostic bonding:¹⁸

* a triplet of concomitant topological objects: a bond critical point, a bond path and an interatomic surface,

* a ring critical point, that is a signature of a structural instability,

* a Laplacian of the electron density of the bond critical point $∇²ρ_{BCP}$ in the 0.15-0.25 a.u. range,

* a negative net charge for an hydrogen atom involved in an agostic bonding,

* a dipolar polarization that is 15-30% larger for an agostic hydrogen compared with a non-agostic one.

However, these criteria are not unambiguous since they cannot be applied to any type of agostic bonding, as it will be further discussed in part III.

<table><caption>Table 1 Geometric criteria used in the literature to distinguish between an agostic and an anagostic interaction⁷</caption><thead><tr><th></th><th>M···H-C agostic</th><th>M···H-C anagostic</th></tr></thead><tbody><tr><td>Distance M···H (Å)</td><td>1.8-2.3</td><td>2.3-2.9</td></tr><tr><td>Angle MHC (°)</td><td>90-140</td><td>110-170</td></tr></tbody></table>

<figure><img src="image_5.png" /></figure>
## 2. Other definitions

In addition with the rigorous initial definition, the term “agostic” is also sometimes used in a broader definition,²¹ to characterize intermolecular interactions, especially in the context of sigma alkane,²²,²³ silane²⁴,²⁵ and sigma boranes.²⁵,²⁶

In an attempt to fully characterize some systems potentially containing weak agostic interactions, experimental and theoretical approaches were applied to numerous and various organometallic complexes characterized by geometric distortions. While in some cases clear agostic bonds were identified, in some other systems the approaches used failed to characterize such interactions. The term “anagostic” was proposed to define systems where an interaction between a σ bond and a metallic center leads to a geometric distortion of the structure whereas some considered criteria are not met to label this interaction as “agostic”.²⁷⁻²⁹

In most of the above-mentioned studies, the agostic interaction occurs in an organometallic complex centered on a metallic atom. In a wider sense however, this concept was employed to qualify situations in which organic molecules interact with metallic surfaces during catalytic processes in the heterogeneous phase. For instance, mechanisms involving agostic interactions were proposed in the chemisorption of hydrogen, hydrocarbons and intermediates on Pt(111).³⁰,³¹ The formation of a particularly rare agostic interaction was even proposed in the context of the water dissociation of a Pt(111) surface, in combination with hydrogen bonding.³²

# III. Theoretical characterization of agostic bonding

## 1. Equilibrium geometry

Since agostic interactions are generally defined on the basis of geometric criteria, geometry optimization may appear as a suitable tool to investigate this type of interactions. Such an approach may complement X-ray structure analyses. In some cases however, the complexes cannot be isolated in the crystalline form (especially in the case of dynamic systems), and a theoretical approach may be the only way to determine whether an agostic interaction is involved in a given reaction pathway. In such a case, particular care should be taken in the choice of the level of theory. In the case of density functional theory calculations, the choice of the functional, the description of the metallic atom as well as the basis set are crucial in the geometrical parameters that will be obtained. Furthermore, in a rigorous theoretical characterization of agostic bonds, such geometry optimizations would only correspond to the first step of a more complete investigation involving either a description in terms of molecular orbitals, or a topological study. As an example, Cho has investigated the agostic structure of titanium methylidene hydride at 15 different levels of theory using Hartree-Fock, Density Functional Theory and Post-Hartree-Fock calculations. He has shown the importance of the choice of the basis set on the agostic distortions, and on the delocalization energies as calculated using the Natural Bond Orbital

(NBO) method (vide infra).³³ The influence of the basis set superposition error (BSSE) was also underlined by several studies. As an example, the α-agostic interaction involved in Cp*Rh(CO)-(η²-alkane) was investigated at various levels of theories,³⁴,³⁵ and the lack of BSSE correction in the conventional ab initio methods leads to an overestimation of the agostic bonding, while some DFT approaches underestimate this interaction.³⁶ Multiconfigurationally quantum chemical methods (complete active space self-consistent field (CASSCF)/second-order perturbation theory (CASPT2)) is a powerful approach to compare the agosticity of several metallic centers. Such an approach was used by Roos et al. to explain experimental results in the case of methylidene metal dihydride complexes.³⁷

## 2. NBO

Agostic interactions may be seen as donor-acceptor systems, with an electron transfer from a σ bond of the ligand to an unoccupied orbital of the metallic atom. In this respect, techniques for the studies of molecular orbitals may shed some light onto agostic bonding. Natural Bond Orbital theory is a powerful approach to characterize donor and acceptor groups in a molecule, and to quantify delocalization energies. In the natural bond orbital approach, a suite of mathematical algorithms are used to diagonalize the density matrix generated during Hartree-Fock, Density Functional Theory or Post Hartree-Fock calculations. These mathematical transformations allow the generation of orbitals in terms of effective atom-like constituents within the molecular environment.³⁸⁻⁴⁰ This approach allows the study of delocalization from atomic orbitals. Toward the natural population analysis, the NBO approach offers a quantitative way to characterize the donor-acceptor relationship. An interesting example of the use of NBO for a quantitative study of agostic bonding is given in the work of Cho et al.³³,⁴¹ Based on a systematic study, a methodology was proposed to analyze the C-H···M agostic bonds in terms of donor-acceptor analysis as defined in the NBO theory.⁴² Thus, the NBO method provides a valuable way to evaluate the agostic character of interactions.⁴³

## 3. QTAIM

In addition to an orbital description of agostic interactions, topological tools were applied in an attempt to qualitatively and quantitatively investigate these bonds. Two topological approaches were particularly applied to the study of agostic bonding: the Bader's quantum theory of atom in molecules (QTAIM) and the electron localization function (ELF). These topological descriptions are based on the dynamical analysis of gradients, from the electronic density, or from functions of electronic density. A vector field is derived from gradients of a potential, and critical points are defined on the basis of the vector fields.

The QTAIM approach is based on a partition of the molecular space into non-overlapping regions within which the local virial theorem is fulfilled. This implies, that the kinetic energy of each of these regions has a definite value which is achieved if and only if the integral of the Laplacian of the electron density, $\nabla^2\rho(r)$, over each region vanishes, a condition which is fulfilled

<figure><img src="image_6.png" /></figure>
for boundaries which are zero-flux surfaces for the gradient of the density.⁴⁴ In this method, gradient fields of the electron density $\nabla\rho(r)$ are studied. Since an atomic center corresponds to a local maximum of electron density, each atomic center acts as an attractor, and field lines define a basin. Basins are separated from each other by zero-flux surfaces. On each zero-flux surface, gradient lines converge toward a critical point called a bond critical point (BCP). The presence of such a BCP is one of the criterion that was proposed by Popelier et al. to characterize an agostic bond in the QTAIM framework.¹⁸,⁴⁵ On the other hand, the identification of BCP may be difficult in the case of weak agostic interactions. In their systematic study of a set of 20 crystal structures potentially characterized by agostic bonding, Thakur and Desiraju noticed that NBO may be more relevant than QTAIM in describing these weak interactions.³³ Numerous studies have now proven that QTAIM can indeed describe $\beta$ C-H agostic interactions,⁴⁶,⁴⁷ but is not suitable for the study of $\alpha$ C-H agostic interactions or weaker interactions such as C-C agostic bonds.⁴⁸ Indeed, the use of new local criteria seem to be the only way to detect weak agostic bonding by means of the QTAIM approach.⁴⁹

## 4. ELF

The ELF approach⁵⁰⁻⁵² intends to provide a partition into basins of attractors which closely match the VSEPR electronic domains, or in other words Lewis's representation. The assumption that groups of electrons can be localized within space filling non-overlapping domains implies that the variances (the squared standard deviation) of the domain populations are minimal. It was pointed out that it is convenient to define a localization function $\eta(r)$ such as:⁵³

$$ \frac{\delta\sigma^2(\bar{N})}{\delta V} = \oint n \cdot \nabla\eta(r)ds = 0 \quad (1) $$

enabling carrying out the partition into basins of the attractors of the gradient dynamical system of $\eta(r)$. It was further proposed⁵⁴ to minimize the Frobenius norm of the covariance matrix:

$$ \|\sigma^2(\Omega_A, \Omega_B)\|_F = \left( \sum_{A,B} (\sigma^2(\Omega_A, \Omega_B))^2 \right)^{1/2} \quad (2) $$

where the matrix element $\sigma^2(\Omega_A, \Omega_B)$ is the covariance of the populations of the basins labelled $\Omega_A$ and $\Omega_B$. The covariance operator is defined as:

$$ \hat{\sigma}^2(\Omega_A, \Omega_B) = \hat{N}(\Omega_A)\hat{N}(\Omega_B) - \bar{N}(\Omega_A)\bar{N}(\Omega_B) \quad (3) $$

with $\hat{N}(\Omega_A) = \sum_{i}^{N} \hat{y}(r_i)(\hat{y}(r_i) = 1, r_i \in \Omega_A$, and $\hat{y}(r_i) = 0, r_i \notin \Omega_A$)

where $\hat{N}(\Omega_A)$ is the population operator introduced by Diner and Claverie⁵⁵ and $\hat{N}(\Omega_A) = \int_{\Omega_A} \rho(r)dr$ is the basin population. The minimization is carried out with the help of the local covariance measure for $\sigma$ spin electrons:

$$ \zeta_{\sigma}^{\text{h}}[\rho; r] = \rho_{\sigma}(r) \left[ (h_{\sigma, \sigma}(r, r'))^2 \right] dr' \quad (4) $$

where $h_{\sigma, \sigma}(r, r')$ is the exchange hole, and can be approximated by an expression proportional to the kernel of the Becke and Edgecombe function.⁵⁰

The covariance can be written as the sum of four spin components, i.e.:

$$ \sigma^2(\Omega_A, \Omega_B) = \sigma_{\alpha\alpha'}^2(\Omega_A, \Omega_B) + \sigma_{\alpha\beta}^2(\Omega_A, \Omega_B) + \sigma_{\beta\alpha}^2(\Omega_A, \Omega_B) + \sigma_{\beta\beta}^2(\Omega_A, \Omega_B) \quad (5) $$

With

$$ \sigma_{\alpha\alpha'}^2(\Omega_A, \Omega_B) = \bar{\Pi}_{\alpha\alpha'}(\Omega_A, \Omega_B) - \bar{N}_{\alpha}(\Omega_A)\bar{N}_{\alpha'}(\Omega_B) \quad (6) $$

$$ \sigma_{\alpha\beta}^2(\Omega_A, \Omega_B) = \bar{\Pi}_{\alpha\beta}(\Omega_A, \Omega_B) - \bar{N}_{\alpha}(\Omega_A)\bar{N}_{\beta}(\Omega_B) \quad (7) $$

$$ \sigma_{\beta\alpha}^2(\Omega_A, \Omega_B) = \bar{\Pi}_{\beta\alpha}(\Omega_A, \Omega_B) - \bar{N}_{\beta}(\Omega_A)\bar{N}_{\alpha}(\Omega_B) \quad (8) $$

$$ \sigma_{\beta\beta}^2(\Omega_A, \Omega_B) = \bar{\Pi}_{\beta\beta}(\Omega_A, \Omega_B) - \bar{N}_{\beta}(\Omega_A)\bar{N}_{\beta}(\Omega_B) \quad (9) $$

where $\bar{\Pi}_{\sigma\sigma'}(\Omega_A, \Omega_B) = \int_{\Omega_A} \int_{\Omega_B} \Pi_{\sigma\sigma'}(r_1; r_2)dr_1dr_2$ is the number of $\sigma\sigma'$ electron pairs between $\Omega_A$ and $\Omega_B$. The opposite spin contributions are expected to be negligible with respect to the same spin ones, and, in fact they are null in the Hartree-Fock approximation. Except for specific delocalized bonding situations such as charge-shift bonds, the covariance between non adjacent basins is negligible. The zero flux surfaces of the gradient field of the spin pair composition for which ELF is an excellent approximation, provide a partition which nearly minimizes the covariance between basins limited a same separatrix.⁵⁶

The ELF partition yields basins of attractors clearly related to Lewis's model: core and valence basins. A core basin surrounds a nucleus with atomic $Z > 2$, it is a single basin for the elements of the second period or the union of the basins belonging to the inner shells for heavier elements. It is labelled C(A) where A is the element symbol. In the study of systems involving transition metal elements it is often useful to consider independently the basins of the metal external core (subvalence) shell. The valence basins are characterized by the number of atomic valence shells to which they participate, or in other words by the number of core basins with which they share a boundary. This number is called the synaptic order. Thus, there are monosynaptic, disynaptic, trisynaptic basins, and so on. Monosynaptic basins, labelled V(A), correspond to the lone pairs of the Lewis model, and polysynaptic basins to the shared pairs of the Lewis model. In particular, disynaptic basins, labeled V(A, X), correspond to two-centre bonds, trisynaptic basins, labeled V(A, X, Y), to three-centre bonds, and so on. The valence shell of a molecule is the union of its valence basins. As hydrogen nuclei are located within the valence shell they are counted as a formal core in the synaptic order because hydrogen atoms have a valence shell. For example, the valence basin accounting for a C-H bond is labeled V(C,H) and called protonated disynaptic. The valence shell of an atom, say A, in a molecule is the union of the valence basins whose label lists contain the element symbol A.

The ELF population analysis provides not only the basin populations and the associated covariance matrix but also the probability of finding n electrons in a given basin and the

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /><figcaption>Fig. 2 QTAIM and ELF basins for the free NCCH₃ molecule.</figcaption></figure>

contribution of the QTAIM basins to the ELF basins. The contribution of the atomic basin of A to the ELF disynaptic basin V(A,B) which is denoted by $N[V(A,B)|A]$ is evaluated by integrating the electron density over the intersection of the V(A,B) basin and of the atomic basin of A. Raub and Jansen⁵⁷ have introduced a bond polarity index defined as:

$$p_{AB} = \frac{\bar{N}[V(A, B)|A] - \bar{N}[V(A, B)|B]}{\bar{N}[V(A, B)|A] + \bar{N}[V(A, B)|B]} \quad (10)$$

As an example, Fig. 2 displays the QTAIM and ELF basins as well as the ELF/QTAIM projection in the case of a simple molecule that is involved as a ligand in the following study, namely the NCCH₃ molecule. To obtain an optimal clarity, only two QTAIM atomic surfaces (H and N) are shown. It is interesting to note that there are no ELF core basins for a hydrogen atom. We reported the ELF valence basin populations as well as the corresponding atomic contributions.

# IV. Topological descriptions, probabilistic approach and qualitative estimators

## 1. Context

ELF has proven to be a valuable tool for characterizing all kinds of chemical bonds. In the case of weak interactions such as 2c-3e and hydrogen bonding, the ELF analysis enables us to describe and classify them in terms of topological properties very close and similar to the traditional chemical concepts.

The analysis of basin population was helpful to distinguish between weak, medium and strong hydrogen bonding. It was shown that the core-valence bifurcation index is a suitable criterion to quantitatively describe these interactions.¹⁷,⁵⁸

In the case of 2c-3e bonding, no disynaptic basins are found. However, this is not the sole criterion to determine whether a bonding exists or not, and once again, core-valence bifurcation index has proved to be a very powerful tool. In addition, in this specific case, a topological delocalization index was helpful in quantifying the electron fluctuation.⁵⁹,⁶⁰

These examples prove that, beyond the BCP, topological approaches can be used in a quantitative way to thoroughly characterize electron localization/delocalization involved in different types of interactions. Below we will present a similar approach to characterize agostic interactions.

## 2. Presentation of the approach

**a. A statistical approach.** The bonding in many molecular systems could not be fully described by considering only a perfect bonding electron localization following the Lewis structure. It is principally due to the indiscernibility of the electron which implies an always present delocalization. This electronic delocalization is not observable, and therefore there is no experimental property that allows its direct measurement.

The multivariate analysis is a basic statistical method enabling one to reveal the correlations between different groups of data. It relies upon the construction of the covariance matrix elements defined by

$$\langle \text{cov}(i, j) \rangle = \langle \hat{i}\hat{j} \rangle - \langle \hat{i} \rangle \langle \hat{j} \rangle \quad (11)$$

where $\langle \hat{i} \rangle$, $\langle \hat{j} \rangle$, and $\langle \hat{i}\hat{j} \rangle$ are the averages of the data values and of their product. The diagonal elements of the covariance matrix are the variances $\sigma^2(\hat{i}) = \langle \hat{i}^2 \rangle - \langle \hat{i} \rangle^2$, which measure the statistical dispersion of the data among the group. The standard deviations $\sigma(\hat{i})$ are the square root of the variances. The off-diagonal matrix elements indicate which data are involved in delocalization.

The bonding in most molecular systems can be described by a strict localization of electron pairs. A more realistic picture, closely related to the concept of resonance, is provided by the superposition of electron multiplets distributed among the basins and therefore accounting for the electron delocalization. The population of a given basin $\Omega_A$ appears accordingly as the average of such n-tuplets weighted by the probability, $P_n(\Omega_A)$ of finding n electron in $\Omega_A$:

$$\bar{N}(\Omega_A) = \sum_{n=1}^{N} n P_n(\Omega_A) \quad (12)$$

A measure of the delocalization is given by the covariance matrix whose diagonal matrix elements are the square of the standard deviations of the basin populations whereas the off matrix elements indicate which basins are involved in the delocalization.⁶¹

In our search for several parameters to qualitatively estimate the strength of agostic bonding, we studied a large number of systems. From this study, it was found that four parameters from QTAIM and ELF taken together may allow a comparison of agostic bonding present in different systems between the X-H bond (X = C, B...) and a metallic center M:

* First of all, we will consider only hydrogen atoms as potentially involved in agostic interactions. In the Brookhart and Green's definition,⁶,⁷ even if this is not the only criterion, the presence of three atomic centers sharing two electrons is the fundamental aspect of agostic bonding. Thus, in the topologic description, the total population of the protonated potentially agostic basin V(H) is an obvious important parameter.

<figure><img src="image_9.png" /></figure>
• The projection of ELF on AIM basins gives some information on the atomic contribution in the agostic protonated basin. Hereafter in this paper, this quantity will be labeled as M/X/H in the case of the X-H agostic bond (X = C or B). Hereafter, this information will be used as a clear indicator of the trisynaptic character of a protonated basin. It is worthy to note that an atomic contribution only makes sense if its value is larger than the numerical error, i.e. 0.01 electron.

• The covariance calculated from the ratio between the basin's population of the potentially agostic H atom and the population of the metallic center core basin C(M) is an important parameter to characterize the interaction between these two basins $Cov(V(H)/C(M))$. The covariance thus obtained from ELF topological analysis gives some insights into the delocalization of electrons between the two atoms. Obviously this value depends on the theoretical description of the metallic center: the covariance is smaller when a pseudo potential is used for the metallic center. Taken together with the previous parameters, a covariance larger than 0.03 (in absolute value) is a proof of an agostic interaction between the X-H σ bond and M.

• Furthermore, in the context of comparison between several levels of calculations, the σ² and the variance calculated with the ELF topological approach are an indication of deviation from perfect localization coming from inter-population. Similar values of variances for a same system calculated at different levels of theories thus allows us to ensure that all the theoretical descriptions are consistent with each other.

Furthermore, it was shown that the electron density ($ρ$) of the bond critical point (BCP) could be related to the bond order and thus the bond strength.⁶² We propose the use of three $ρ(BCP)$ to gain some insight into the strength of the agostic interaction. In line with Popelier and Logothetis,¹⁸ we suggest the use of the $ρ(BCP)$ of the M-H$_β$ bonding, when it exists. Additionally, the comparison of the $ρ(BCP)$ values associated with the C-H$_agostic$ and C-H$_free$ allows us to estimate the weakening of the C-H$_agostic$ bond caused by the agostic interaction. Two conditions are necessary for the use of the $ρ(BCP(C-H$_agostic$))$ and $ρ(BCP(C-H$_free$))$ values. First, the carbon bearing the hydrogen atom potentially involved in an agostic interaction should also bear an additional non-agostic hydrogen atom. This condition is often fulfilled. Furthermore, neither of the H$_agostic$ and H$_free$ atoms should be involved in another non-covalent interaction. Provided that these conditions are satisfied, we can propose the following reference values. In the alkyl complexes, the $ρ(BCP(C-H$_free$))$ are characterized by

$$0.28 \leq ρ(BCP(C-H_{free})) \leq 0.29$$

whereas the C-H$_agostic$ are characterized by smaller electron densities at the BCP, with

$$0.200 \leq ρ(BCP(C-H_{agostic})) \leq 0.27$$

From the ELF point of view, a protonated basin is considered as a trisynaptic basin when its population originates from three atomic centers. An H-agostic bond is a protonated

trisynaptic basin whose population is around $2e⁻$. Consequently, it corresponds to the traditional $3c-2e$ interaction in chemistry.

These parameters were not only chosen for their ability to describe an agostic interaction, but also for their phenomenological significances. Indeed, the increase in covariance $Cov(V(H)/C(M))$ signifies that the delocalization of electrons between the metallic center and the agostic protonated basin increases, concomitantly with a weakening of the bond between the hydrogen and other atom X. The H-X bond is thus activated and the agostic character of the interaction between H and M is stronger. In a limit case, when the agostic bonding becomes stronger and stronger, the ELF basin evolves from a trisynaptic to a disynaptic basin, and the covariance increases till the formation of a covalent M-H bond. This corresponds to a metallic hydride.

Moreover, the strength of an agostic H-X bond may be estimated from the QTAIM properties calculated at the H-X bond critical point (BCP). Interestingly, the charge density at the agostic H-X BCP compared to that of a non-agostic H-X bond in the complex (and/or compared to the charge density at the H-X BCP of free ligand) provides an indication to the strength of the agostic interaction.

**b. Choice of a representative set of agostic bonding.** When Popelier and Logothetis carried out their seminal work on the topological characterization of agostic bonding, they have chosen TiCl$_2$-alkyl complexes as models, and they investigated agostic bonding by means of the QTAIM approach.¹⁸ Thus, these simple models have already covered four types of potentially agostic bonds. The QTAIM method was helpful in characterizing an agostic bond in the case of $β-C-H$, $γ-C-H$ and $C-C$ interactions with the metallic center, but not in the case of an interaction with an $α-C-H$. These systems therefore appear as references in the context of setting a new theoretical approach for the characterization of agostic bonding.

From this study as well as from further theoretical investigations, it was found that QTAIM suitably describe $β-C-H$ and $γ-C-H$ agostic bonds, but not the $α-C-H$ ones, even in cases for which experimental data tend to prove that $α-C-H$ agostic bonds were indeed formed.

Some compounds were observed under different agostomers. We will choose a few examples of agostomers based on experimental investigations carried out by Baird et al.⁶³,⁶⁴ In the case of $[Cp_2TiCH_2CHMe t-Bu]⁺$, they observed that the $α$-agostic isomer is preferentially formed although a $β$-agostic isomer could have been formed. Such a situation is relatively rare: when both $α$-and $β$-agostic isomers may be formed, generally the $β$-agostic form is more stable. Another alkyl-titanium complex that may exist as $β$- and $γ$-agostomers will be considered.

The complex formed between acetonitrile and zero-valent nickel that may lead to the formation of a C-H agostic bond⁶⁵ will be investigated.

In an attempt to understand the parameters that influence the formation of agostic bonds, the effect of co-ligands, small changes in the structures and the nature of the metallic center, will be investigated. The effect of the presence of a co-ligand

<figure><img src="image_10.png" /></figure>
will be topologically investigated based on the example of a rhodium thiophosphoryl pincer complex studied by Milstein *et al.*⁶⁶

Different titanium complexes will be compared, based on the compounds studied by Popelier,¹⁸ Baird⁶³,⁶⁴ and Mc Grady.⁶⁷ The comparison of complexes studied by Mc Grady⁶⁷ and Forster⁶⁸ will allow the study of influence of the nature of the metallic center. Further examples derived from the model compounds of Popelier¹⁸ as well as from Sabo-Etienne⁶⁹ will also be discussed.

Lastly, interactions involving heteroatoms will be considered. Indeed, such interactions are sometimes considered as “agostic” whereas some authors consider that all the contribution from a σ bond and an unoccupied orbital of a hypovalent transition metal center cannot be classified under an unique appellation. Our aim here is to determine whether the methodology above presented is able to differentiate between these types of interactions, or whether these interactions are due to a similar effect. For this purpose, titanocene and zirconocene amidoborane complexes,⁶⁷,⁶⁸ dimethylaminoborane complexes,⁶⁹ and mesitylborane complexes⁷⁰ will be considered. The last examples correspond to intermolecular interactions, whereas in the other systems, the weak interaction between a σ X-H bond and the metal center should be considered as an intramolecular interaction.

# V. Application to a representative set of examples

## 1. Influence of the level of theory

To begin with, the influence of the level of theory on the parameters calculated with the ELF topology will be checked. To this end, DFT calculations were carried out, and three different types of functionals were selected:

* B3LYP because this is one of the most popular hybrid functional,

* PBE0 because this is a non-empirical hybrid functional⁷¹ that was widely employed in the context of agostic interactions,

* TPSSh because this meta-GGA hybrid functional can be used for reference calculations, when combined with a suitable basis set.⁷²

In combination with these functionals, six different basis sets were selected:

* the 6-311++G(2d,2p) basis set, in combination with a pseudo-potential LANL-2TZ-f including a triple ξ and an additional diffuse f function for the metallic center,

* the 6-311++G(2d,2p) basis set, in combination with a pseudo-potential LANL-2TZ-p including a triple ξ and an additional diffuse p function for the metallic center,

* the 6-311++G(2d,2p) basis set, in combination with a pseudo-potential LANL-2DZ including a double ξ function for the metallic center,

* the 6-311++G(2d,2p) basis set, in combination with a pseudo-potential LANL-2DZ including a double ξ function for the metallic center,

* the 6-311++G(2d,2p) basis set, without any pseudo potential for the metallic center,

* the 6-311++G(3df,3pd) basis set, without any pseudo potential for the metallic center.

For these tests of the influence of the levels of theories, a simple model molecule taken from the study of Popelier and Logothetis was selected (see below, [TiCl₂CH₂CH₃]⁺).

In a first series of tests, the molecule was re-optimized at each level of theory prior to the topological investigation. In a second series of tests, the complex was optimized using the highest level of theory, namely B3LYP, PBE0, or TPSSh/ 6-31++G(3df,3pd), and ELF calculations were then carried out on the wave function of the single point geometries using the LanL2DZ as the basis set.

Optimized geometries are compared in Table 2 and the topological parameters are summed up in Table 3.

The values presented in Table 2 clearly show that the geometry of the complex is correctly described even with the B3LYP functional, when used in combination with a relatively large basis set including 2d and 2p polarization functions. Indeed, if we compare the distances and the angles calculated at the B3LYP/6-311++G(2d,2p) level of theory with the reference values (calculations at the TPSSh/6-311++G(3df,3pd) level of theory), the errors are 0.001 Å, 0.026 Å and 0.3° for the *d*(Ti-C), *d*(Ti-H<sub>agostic</sub>), and θ, respectively.

Table 3 clearly shows that all the selected levels of theory identify an agostic bonding in this simple TiCl₂-alkyl system. Indeed,

* the total population of the valence basin of Hβ is slightly smaller than 2e⁻,

* the variance value (σ²) does not depend on the level of theory which indicates the stability of the ELF topological procedure in partitioning the molecular space,

* in all the cases, the contribution of the titanium atom to this basin is in the range of 0.07–0.08 e⁻ provided to use an explicit triple zeta quality basis set for metallic center,

<table><caption>**Table 2** Distances (in Å) and angles (in °) for the [TiCl₂CH₂CH₃]⁺ complex optimized at different levels of theory</caption><thead><tr><th>Level of theory</th><th>d(Ti-C) (Å)</th><th>d(Ti-H<sub>agostic</sub>) (Å)</th><th>θ (°)</th></tr></thead><tbody><tr><td>B3LYP/6-311++G(2d,2p)</td><td>2.008</td><td>2.034</td><td>113.8</td></tr><tr><td>B3LYP/6-311++G(3df,3pd)</td><td>2.008</td><td>2.030</td><td>113.8</td></tr><tr><td>PBE0/6-311++G(3df,3pd)</td><td>1.986</td><td>2.015</td><td>114.3</td></tr><tr><td>TPSSh/6-311++G(3df,3pd)</td><td>2.009</td><td>2.008</td><td>113.5</td></tr></tbody></table>

<figure><img src="image_11.png" /></figure>
<table><caption>**Table 3** Topological and geometrical parameters obtained at several levels of theories. The total population of the agostic protonated basin V(Hβ) and the contribution from Ti, C and H (Ti/Cβ/Hβ), the covariance cov(V(Hβ)/C(Ti)) and the variance (σ²) are given. Distances (in Å) and angles (in degrees) are given for the optimized structures</caption><thead><tr><th>Level of theory</th><th>V(H<sub>β</sub>)</th><th>Ti/C<sub>β</sub>/H<sub>β</sub></th><th>Cov(V(H<sub>β</sub>)/C(Ti))</th><th>σ<sup>2</sup></th><th>d(Ti-H<sub>agostic</sub>), θ, d(Ti-C)</th></tr></thead><tbody><tr><td rowspan="7">B3LYP</td><td>LanL2DZ</td><td>1.89</td><td>0.04/0.82/1.03</td><td>-0.08</td><td>0.72</td><td>2.116, 115.6, 1.982</td></tr><tr><td>6-31++G(2d,2p)</td><td>1.95</td><td>0.08/0.78/1.09</td><td>-0.09</td><td>0.74</td><td>2.034, 113.8, 2.003</td></tr><tr><td>6-311++G(2d,2p)</td><td>1.93</td><td>0.08/0.76/1.09</td><td>-0.09</td><td>0.73</td><td>2.034, 113.8, 2.008</td></tr><tr><td>6-311++G(2d,2p)/LanL2DZ</td><td>1.92</td><td>0.03/0.78/1.11</td><td>-0.09</td><td>0.72</td><td>2.025, 113.7, 1.994</td></tr><tr><td>6-311++G(2d,2p)/LanL2DZ</td><td>1.92</td><td>0.04/0.75/1.13</td><td>-0.09</td><td>0.73</td><td>2.005, 113.4, 1.996</td></tr><tr><td>6-311++G(2d,2p)/LanL2TZ(P)</td><td>1.90</td><td>0.04/0.73/1.13</td><td>-0.10</td><td>0.72</td><td>2.009, 113.6, 2.004</td></tr><tr><td>6-311++G(2d,2p)/LanL2TZ(F)</td><td>1.90</td><td>0.05/0.72/1.13</td><td>-0.10</td><td>0.72</td><td>2.005, 113.4, 2.002</td></tr><tr><td rowspan="6">PBE0</td><td>6-311++G(2d,2p)</td><td>1.93</td><td>0.07/0.78/1.08</td><td>-0.09</td><td>0.74</td><td>2.019, 114.4, 1.983</td></tr><tr><td>6-311++G(2d,2p)</td><td>1.92</td><td>0.07/0.77/1.08</td><td>-0.10</td><td>0.74</td><td>2.018, 114.3, 1.986</td></tr><tr><td>6-311++G(2d,2p)/LanL2DZ</td><td>1.91</td><td>0.04/0.75/1.12</td><td>-0.10</td><td>0.73</td><td>2.002, 117.0, 1.975</td></tr><tr><td>6-311++G(2d,2p)/LanL2DZ</td><td>1.91</td><td>0.04/0.75/1.12</td><td>-0.10</td><td>0.73</td><td>1.985, 113.9, 1.977</td></tr><tr><td>6-311++G(2d,2p)/LanL2TZ(P)</td><td>1.90</td><td>0.04/0.74/1.12</td><td>-0.10</td><td>0.74</td><td>1.990, 114.1, 1.985</td></tr><tr><td>6-311++G(2d,2p)/LanL2TZ(F)</td><td>1.90</td><td>0.05/0.73/1.12</td><td>-0.10</td><td>0.73</td><td>1.989, 114.1, 1.982</td></tr><tr><td rowspan="6">TPSSH</td><td>6-311++G(2d,2p)</td><td>1.93</td><td>0.08/0.74/1.11</td><td>-0.10</td><td>0.73</td><td>2.016, 113.6, 2.006</td></tr><tr><td>6-311++G(2d,2p)</td><td>1.93</td><td>0.08/0.74/1.11</td><td>-0.10</td><td>0.72</td><td>2.011, 113.5, 2.009</td></tr><tr><td>6-311++G(2d,2p)/LanL2DZ</td><td>1.93</td><td>0.03/0.77/1.13</td><td>-0.10</td><td>0.73</td><td>2.001, 113.4, 1.995</td></tr><tr><td>6-311++G(2d,2p)/LanL2DZ</td><td>1.90</td><td>0.03/0.72/1.15</td><td>-0.09</td><td>0.72</td><td>1.980, 113.1, 1.996</td></tr><tr><td>6-311++G(2d,2p)/LanL2TZ(P)</td><td>1.91</td><td>0.04/0.72/1.15</td><td>-0.10</td><td>0.73</td><td>1.985, 113.3, 2.005</td></tr><tr><td>6-311++G(2d,2p)/LanL2tz(F)</td><td>1.91</td><td>0.04/0.72/1.15</td><td>-0.10</td><td>0.72</td><td>1.985, 113.3, 2.002</td></tr><tr><td rowspan="3">SP</td><td>B3LYP/LanL2DZ//B3LYP/6-311++G(3df,3pd)</td><td>1.83</td><td>0.03/0.73/1.07</td><td>-0.08</td><td>0.72</td><td>Single-point calculations</td></tr><tr><td>PBE0/LanL2DZ//PBE0/6-311++G(3df,3pd)</td><td>1.89</td><td>0.03/0.82/1.04</td><td>-0.08</td><td>0.73</td><td></td></tr><tr><td>TPSSH/LanL2DZ//TPSSH/6-311++G(3df,3pd)</td><td>1.89</td><td>0.02/0.81/1.06</td><td>-0.08</td><td>0.72</td><td></td></tr></tbody></table>

• the calculated covariance is close to -0.1, clearly indicating a delocalization of the electrons between agostic hydrogen and a metallic center.

These four criteria should be considered all together, the presence or the absence of one of the criterion is not sufficient to drive any conclusion.

We would like to emphasize the fact that, whatever the level of theory selected, the topological analysis of the agostic bonding leads to similar results. Thus, the criteria selected are robust toward the level of theory.

This study thus show that the B3LYP/6-311++G(2d,2p) can be used for a topological investigation of agostic interactions. This level of theory was selected for all the further studies presented below.

## 2. Study of different agostomers

Let us first consider alkyl titanocene compounds [Cp₂TiCH₂CH-Me(t-Bu)]⁺ studied by Baird *et al.*⁶³ Two agostic isomers were characterized: α- and β-agostomers (Fig. 3).

The experimental study clearly shows that, between two diastereoisomers α- and β-agostic, the α-agostic one is more stable. The Table 4 shows, without any surprise, that a BCP is indeed obtained in the case of the β-agostic isomer that is not the case for the α-agostic one. Nevertheless, it is worth noting that the electron density at the C-H bond slightly decreases by 0.05 a.u. due to the agostic deformation for both α- and β-agostic isomers. Simultaneously, the Laplacian of charge density (∇²ρ) and the energy density (H) at the BCP increase, leading to the reduction of the covalent character of the C-H bond.

<figure><img src="image_12.png" /><figcaption>**Fig. 3** α- and β-agostomers observed for an alkyl titanocene compound.⁶³</figcaption></figure>

To further characterize these two isomers, the methodological approach above proposed was applied, and Table 5 summarizes the ELF investigation of these two isomers.

The results presented in Table 5 clearly show that, despite the absence of BCP in the alpha agostomer, the agostic bonding is indeed described by the combined ELF/QTAIM studies. Indeed, the total population of the valence basin of Hₓ is below 2e⁻, the contribution of the Ti atom to the protonated valence basin of Hₓ is 0.04, which is not negligible, and the covariance between two basins Vₕₓ/Cₘ is close to -0.1. In comparison, the total population of the agostic H atom is closer to 2e⁻ and the covariance is slightly smaller (in absolute value) in the case of the β-agostic isomer. As a conclusion, the agostic character is slightly more pronounced in the α-agostic isomer compared with the β-agostic one, for this specific case. Furthermore, these two agostic bonds are relatively weak.

Baird *et al.*⁶⁴ also studied an alkyl-titanium complex that exist under two isomeric forms, one presenting a β-agostic

<figure><img src="image_13.png" /></figure>
<figure><img src="image_14.png" /></figure>

<table><caption>**Table 4** Search for BCP's and RCP's in the α- and β- agostic alkyl titanocene compounds studied by Baird *et al.*<sup>63</sup></caption><thead><tr><th>α-agostomer</th><th>β-agostomer</th></tr></thead><tbody><tr><td>BCP(Ti-Hα): does not exist</td><td>BCP(Ti-Hβ): ρ = 0.03, ∇<sup>2</sup>ρ = +0.11, H = 0.00, ε = 0.37</td></tr><tr><td>BCP(Ti-Cα): ρ = 0.10, ∇<sup>2</sup>ρ = +0.07, H = -0.04</td><td>RCP(Ti-Hβ-Cβ): ρ = 0.03, ∇<sup>2</sup>ρ = +0.13</td></tr><tr><td>BCP(Cα-Hα): ρ = 0.24, ∇<sup>2</sup>ρ = -0.70, H = -0.23</td><td>BCP(Ti-Cα): ρ = 0.10, ∇<sup>2</sup>ρ = +0.03, H = -0.04</td></tr><tr><td>BCP(Cα-H): ρ = 0.29, ∇<sup>2</sup>ρ = -1.05, H = -0.31</td><td>BCP(Cα-H): ρ = 0.29, ∇<sup>2</sup>ρ = -1.04, H = -0.31</td></tr><tr><td>BCP(Cα-Cβ): ρ = 0.24, ∇<sup>2</sup>ρ = -0.51, H = -0.19</td><td>BCP(Cα-Cβ): ρ = 0.24, ∇<sup>2</sup>ρ = -0.49, H = -0.19</td></tr><tr><td>BCP(Cβ-H): ρ = 0.29, ∇<sup>2</sup>ρ = -1.03, H = -0.30</td><td>BCP(Cβ-Hβ): ρ = 0.24, ∇<sup>2</sup>ρ = -0.69, H = -0.22</td></tr></tbody></table>

<table><caption>**Table 5** Topological characterization of the α- and β- agostic alkyl titanocene (Cp = cyclopentadienyl) compounds studied by Baird *et al.*<sup>63</sup></caption><thead><tr><th rowspan="2">Isomer</th><th>BuMe<sup>1</sup>HC</th><th>Cp</th></tr><tr><th>Ti-alpha-Baird</th><th>Ti-beta-Baird</th></tr></thead><tbody><tr><td>V(H)</td><td>1.89</td><td>1.9</td></tr><tr><td>M/X/H</td><td>0.04/0.78/1.07</td><td>0.04/0.80/1.11</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.09</td><td>-0.06</td></tr><tr><td>d(Ti-H<sub>agostic</sub>) (Å)</td><td>2.05</td><td>2.01</td></tr><tr><td>d(C-H<sub>agostic</sub>) (Å)</td><td>1.14</td><td>1.17</td></tr></tbody></table>

<table><caption>**Table 6** Topological characterization of the β- and γ- agostic alkyl titanium complexes studied by Baird *et al.*<sup>64</sup></caption><thead><tr><th rowspan="2">Isomer</th><th>SiMe<sub>3</sub></th><th>Ti-gamma-Si-Baird</th></tr><tr><th>Ti-beta-Si-Baird</th><th>Ti-gamma-Si-Baird</th></tr></thead><tbody><tr><td>V(H)</td><td>1.93</td><td>1.95</td></tr><tr><td>M/X/H</td><td>0.06/0.77/1.10</td><td>0.05/0.81/1.09</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.08</td><td>-0.07</td></tr><tr><td>d(Ti-H<sub>agostic</sub>) (Å)</td><td>2.09</td><td>2.00</td></tr><tr><td>d(C-H<sub>agostic</sub>) (Å)</td><td>1.16</td><td>1.15</td></tr></tbody></table>

bonding, and another presenting a γ- agostic bonding. Table 6 summarizes the ELF/QTAIM characteristics of these two isomers. As in the previous case, the statistical parameters of ELF clearly describe both agostic bonds. These bonds are characterized by BCP's and RCP's in the QTAIM description. The quantitative study of these bonds shows that the σ C-H interactions are relatively weak, as it was the case in the previous α- and β- agostic isomers.

As a conclusion, these examples show that the topological tools of the ELF approach indeed allows the characterization of α-, β- and γ- agostic bonds.

### 3. Topological characterization of a representative set of C-H agostic bonds

A set of complexes was used to probe the agostic character of different features of M···H-C bonding by means of the topological criteria previously presented, and the selected topological criteria are summarized in Table 7.

For all the species the total population of the valence basin of the agostic hydrogen atom is in the 2 ± 0.15 e⁻ range.

Both the covariance Cov(V(H)/C(M)) and the atomic contributions in the valence basin of H(M,X,H) fluctuate depending on the complex.

Table 7 presents complexes containing intramolecular σ C-H agostic bonding. However, the Rh-H₂CO-Milstein<sup>66</sup> and Rh-butene-Milstein<sup>66</sup> complexes could not be easily classified as β agostic species. We note in passing with these two examples that a change in the nature of a co-ligand may strongly affect the agostic character of a C-H bond. Indeed a weak agostic character is topologically predicted in the case of the Rh-butene-Milstein<sup>66</sup> complex, whereas a strong one is predicted for the analogous compound containing H₂CO instead of butene as the co-ligand. These two cases will be further discussed hereafter.

A close look of the data reported in Table 7 shows some trends which be summarized as follows:

* In the view of the metal contribution in the protonated basin, we can classify the complexes into four categories: (1) M = 0.01 e corresponding to an undefined case, (2) 0.01 < M < 0.05 for weak-medium agostic bonding, (3) 0.05 < M < 0.20 for medium-strong agostic interaction, and (4) M > 0.20 for almost pre-dissociated C-H<sub>agostic</sub> or pre-hydride H<sub>agostic</sub>-M.

* One can note that the calculated harmonic vibrational frequency of the C-H<sub>agostic</sub> oscillator is always red-shifted with respect to that of C-H<sub>free</sub>. To a certain extent, the amount of this red-shift reflects the strength of the agostic interaction. It is interesting to note that the harmonic vibrational frequency of C-H<sub>free</sub> is sometimes blue-shifted with respect to the C-H frequency in the free ligand. However, on the ground of C-H vibrational frequency red-shift one can easily distinguish three categories of agostic species: red-shift ≈ 2% for the weakest agostic compound (CpTiNiPr₂Cl₂-McGrady), red-shift ≈ 40% for the strongest agostic compound (Rh-H₂CO-Milstein), and 5% < red-shift < 40% for the other compounds going from weak to strong cases.

* The case of the CpTiNiPr₂Cl₂-McGrady compound: this compound has been previously considered as an agostic case by McGrady *et al.*<sup>67</sup> and also by Scherer and coworkers.<sup>73</sup> We would like to emphasize that the very low metal contribution (0.01 e) in the population of V(H) makes actually impossible to decide the presence or absence of an agostic interaction, because of the numerical uncertainty of our ELF analysis which is just equal to 0.01 e.

This is also consistent with the geometrical properties of the complex: if we refer to the criteria summarized in Table 1, the CpTiNiPr₂Cl₂-McGrady compound is anagostic (d(M-H) > 2.3 Å).

In order to check the possible effect of the dispersion contribution in the electronic structure, we also optimized
<table><caption>**Table 7** Topological characterization of a set of M···H–C complexes reported in the literature as agostic compounds. Optimized geometries have been found at the B3LYP/6-311++G(2d,2p) level of theory</caption><thead><tr><th>Compound</th><th>Ti-Popelier<sup>18</sup></th><th>Ti-alpha-Baird<sup>63</sup></th><th>Ti-beta-Baird<sup>63</sup></th><th>EtTiCl<sub>3</sub>(dmpe)-McGrady<sup>67</sup></th></tr></thead><tbody><tr><td>V(H)</td><td>1.93</td><td>1.90</td><td>1.96</td><td>1.97</td></tr><tr><td>M/X/H</td><td>0.08/0.76/1.09</td><td>0.04/0.79/1.07</td><td>0.04/0.81/1.11</td><td>0.03/0.94/1.00</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.09</td><td>-0.09</td><td>-0.06</td><td>-0.04</td></tr><tr><td>d(Ti-H<sub>agostic</sub>) (Å)</td><td>2.034</td><td>2.049</td><td>2.157</td><td>2.183</td></tr><tr><td>d(C-H<sub>agostic</sub>) (Å)</td><td>1.146</td><td>1.140</td><td>1.151</td><td>1.115</td></tr><tr><td>ω(C-H<sub>agostic</sub>) (cm<sup>-1</sup>)</td><td>2534</td><td>2589</td><td>2429</td><td>2815</td></tr><tr><td>d(C<sub>z</sub>-H<sub>free</sub>) (Å)</td><td>1.0874, 1.0874</td><td>1.085</td><td>1.083</td><td>1.083</td></tr><tr><td>ω(C-H<sub>free</sub>) (cm<sup>-1</sup>)</td><td>3087-3172</td><td>3037-3134</td><td>3034-3197</td><td>3067-3177</td></tr><tr><td colspan="5">QTAIM topological parameters for agostic compound: ρ, ∇<sup>2</sup>ρ, H(ρ) in a.u.</td></tr><tr><td>BCP(C-H<sub>β</sub>)</td><td>0.233, -0.67, -0.22</td><td>0.238, -0.70, -0.30</td><td>0.235, -0.69, -0.22</td><td>0.262, -0.86, -0.26</td></tr><tr><td>BCP(C-H<sub>free</sub>)</td><td>0.287, -1.07, -0.30</td><td>0.288, -1.05, -0.31</td><td>0.288, -1.05, -0.31</td><td>0.285, -1.03, -0.31</td></tr><tr><td>BCP(Ti-H<sub>β</sub>)</td><td>0.045, +0.14, 0.00</td><td>Does not exist</td><td>0.029, +0.11, 0.00</td><td>Does not exist</td></tr><tr><td colspan="5">Some relevant parameters for free ligand:</td></tr><tr><td>d(C-H) (Å)</td><td>1.091</td><td>1.096</td><td></td><td>1.0908</td></tr><tr><td>ω(C-H) (cm<sup>-1</sup>)</td><td>3034-03101</td><td>2991-3110</td><td></td><td>3034-3101</td></tr><tr><td>BCP(C-H<sub>β</sub>)</td><td>0.281, -1.01, -0.30</td><td>0.282, -1.01, -0.30</td><td></td><td>0.281, -1.01, -0.30</td></tr><tr><td colspan="5">Electron density at BCP(C-H) in function of d(C-H)</td></tr><tr><td colspan="5">R<sup>2</sup> = 0.996</td></tr><tr><td colspan="5">R<sup>2</sup> = 0.982</td></tr><tr><td colspan="5">R<sup>2</sup> = 0.989</td></tr></tbody></table>

<table><thead><tr><th>Compound</th><th>CpTiNiPr<sub>2</sub>Cl<sub>2</sub>-McGrady<sup>67</sup></th><th>Rh-butene-Milstein<sup>66</sup></th><th>Rh-H<sub>2</sub>CO-Milstein<sup>66</sup></th><th>Ni-Jones<sup>65</sup></th></tr></thead><tbody><tr><td>V(H)</td><td>2.06</td><td>2.01</td><td>2.12</td><td>2.01</td></tr><tr><td>M/X/H</td><td>0.01/1.07/0.98</td><td>0.05/1.02/0.95</td><td>0.23/0.93/0.96</td><td>0.07/0.98/0.96</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.02</td><td>-0.15</td><td>-0.38</td><td>-0.13</td></tr><tr><td>d(Ti-H<sub>agostic</sub>)</td><td>2.363</td><td>1.939</td><td>1.665</td><td>1.811</td></tr><tr><td>d(C-H<sub>agostic</sub>)</td><td>1.095</td><td>1.131</td><td>1.222</td><td>1.126</td></tr><tr><td>ω(C-H<sub>agostic</sub>)</td><td>3006</td><td>2612</td><td>1950</td><td>2665</td></tr><tr><td>d(C<sub>z</sub>-H<sub>free</sub>)</td><td>1.088, 1.090, 1.090</td><td>1.088</td><td>1.088</td><td>1.088</td></tr><tr><td>ω(C-H<sub>free</sub>)</td><td>3040-3130</td><td>3174-3192</td><td>3174-3182</td><td>3070-3119</td></tr><tr><td colspan="5">QTAIM topological parameters for agostic compound: ρ, ∇<sup>2</sup>ρ, H(ρ) in a.u.</td></tr><tr><td>BCP(C-H<sub>agostic</sub>)</td><td>0.287, -1.05, -0.30</td><td>0.250, -0.77, -0.24</td><td>0.202, -0.47, -0.16</td><td>0.256, -0.82, -0.25</td></tr><tr><td>BCP(C-H<sub>free</sub>)</td><td>0.284, -1.03, -0.30</td><td>0.281, -0.97, -0.28</td><td>0.282, -0.98, -0.85</td><td>0.285, -1.04, -0.30</td></tr><tr><td>BCP(M-H<sub>agostic</sub>)</td><td>Does not exist</td><td>0.055, +0.20, -0.01</td><td>0.107, +0.26, -0.05</td><td>0.057, +0.23, -0.01</td></tr><tr><td colspan="5">Some relevant parameters for free ligand:</td></tr><tr><td>d(C-H) (Å)</td><td>1.093</td><td>1.092</td><td></td><td>1.089</td></tr><tr><td>ω(C-H) (cm<sup>-1</sup>)</td><td>3020-3040</td><td>3050</td><td></td><td>3060-3128</td></tr><tr><td>BCP(C-H<sub>β</sub>)</td><td>0.288, -1.06, -0.31</td><td>0.277, -0.94, -0.28</td><td></td><td>0.283, -1.03, -0.30</td></tr></tbody></table>

<figure><img src="image_15.png" /></figure>
<table><caption>Table 7 (continued)</caption><thead><tr><th>Compound</th><th>CpTiNiPr<sub>2</sub>Cl<sub>2</sub>-McGrady<sup>67</sup></th><th>Rh-butene-Milstein<sup>66</sup></th><th>Rh-H<sub>2</sub>CO-Milstein<sup>66</sup></th><th>Ni-Jones<sup>65</sup></th></tr></thead><tbody><tr><td>Electron density at BCP(C-H) in function of d(C-H)</td><td>Borderline case: if there was an agostic case, it should correspond to a very weak interaction!</td><td></td><td></td><td></td></tr></tbody></table>

the studied structure using two hybrid functionals (wB97XD and B2PLYPD3) which are suitable to treat the very weak non-covalent interactions.

As shown by the results reported in Table 8, topological differences between the C<sub>β</sub>-H<sub>β</sub> so-called agostic bond and C<sub>β</sub>-H<sub>free</sub> within the same compound are minor so that we can confidently exclude a dominant agostic interaction within this complex. This description is also supported by a very small vibrational ω(C-H) frequency shift (less than 2%) with respect to free ω(C-H). Furthermore, we found a BCP between the so-called agostic H atom and one of the two chlorine atoms. At this latter BCP, the electron density is equal to 0.015 e belonging to the hydrogen bonded range.<sup>74</sup>

* Agostic bonding and QTAIM bond critical point: we note that the presence of a BCP(H-Ti) in the case of the EtTiCl<sub>3</sub>(dmpe)-McGrady compound actually depend on the level of theory. Indeed, we found a BCP only at the BP86/6-311++G(d,p) level, whereas there is no BCP(H-Ti) when we use B3LYP, PBE0 or BP86 with 6-311++G(2d,2p) as the basis set.

* Concerning the strength of an agostic bond, we note that the electron density at the BCP(C-H agostic) decreases when the agosticity increases. This trend is graphically shown for each compound in Table 7. A global linear regression graph for all the species will be discussed in Section VI.

<table><caption>**Table 8** Effect of the density functional on the characterization of the CpTiNiPr<sub>2</sub>Cl<sub>2</sub>-McGrady complex using the 6311++G(2d,2p) basis set. No BCP(Ti-H<sub>β</sub>) was found for any of the three functional</caption><thead><tr><th></th><th>wB97XD</th><th>B2PLYPD3</th><th>B3LYP</th></tr></thead><tbody><tr><td>d(C<sub>β</sub>-H<sub>β</sub>) (Å)</td><td>1.096</td><td>1.095</td><td>1.095</td></tr><tr><td>ω(C-H) (cm<sup>-1</sup>)</td><td>3035</td><td></td><td>3006</td></tr><tr><td>BCP(C<sub>β</sub>-H<sub>β</sub>): ρ, ∇<sup>2</sup>ρ in a.u.</td><td>0.286, -1.04</td><td>0.287, -1.05</td><td>0.284, -1.03</td></tr><tr><td>d(C<sub>β</sub>-H<sub>free</sub>) (Å)</td><td>1.091</td><td>1.090</td><td>1.090</td></tr><tr><td>ω(C-H) (cm<sup>-1</sup>)</td><td>3077</td><td></td><td>3040</td></tr><tr><td>BCP(C<sub>β</sub>-H<sub>free</sub>): ρ, ∇<sup>2</sup>ρ in a.u.</td><td>0.290, -1.067</td><td>0.291, -1.08</td><td>0.287, -1.05</td></tr></tbody></table>

## 4. Topological characterization of a representative set of B-H agostic bonds

In Table 9 are gathered five complexes containing intra-molecular σ B-H agostic bonds. These complexes have been a subject of experimental and/or theoretical study.<sup>65,67,68</sup>

For all the amidoborane titanocene or zirconocene complexes, where the nitrogen atom is at the α-position and the boron atom at the β-position, the presence of a B-H protonated basin containing a metallic contribution ranging from 3% to 6% of the V(H) population is a clear indicator of the existence of a β-agostic bond. This is supported by the bond lengthening and frequency red-shift of the B-H<sub>agostic</sub> bond. Compared to the C-H agostic bonding, the B-H agosticity should be considered as medium to strong interaction. This consideration is naturally in line with the decrease of the electron density at the B-H<sub>agostic</sub> bond critical point. It is graphically evidenced on the linear regression graph (Table 9).

Globally, the ELF/QTAIM criteria led to a homogeneous and consistent description of the bonds thus supporting the use of “agostic” for both σ C-H···M and σ B-H···M intramolecular bonding.

## 5. Parameters influencing the agostic character of bonding

To further investigate σ C-H β-agostic interactions and the parameters affecting such bonding, an interesting case-study of rhodium thiophosphoryl pincer investigated by Milstein *et al.* will be detailed below. Indeed, the authors reported the formation of identical agostic bonds with Rh that also interacts either with but-2-ene or with formaldehyde. The agostic hydrogen atoms are characterized by largely different geometrical parameters, and the topological investigation proves that the nature of the co-ligand R is of paramount importance on the agostic interaction (see Scheme 1).

To better understand the co-ligand influence we will carefully analyze and compare the topological properties of the Rh-butene-Milstein compound to those of Rh-H<sub>2</sub>CO-Milstein.

<figure><img src="image_16.png" /></figure>
<table><caption>**Table 9** Topological characterization of a set of M···H–B complexes reported in the literature as agostic compounds. Optimized geometries have been found at the B3LYP/Def2-TZVPD level of theory</caption><thead><tr><th>Compound</th><th>ZrClNH<sub>2</sub>BH<sub>3</sub>-a-<br>Forster<sup>68</sup></th><th>ZrClNH<sub>2</sub>BH<sub>3</sub>-b-<br>Forster<sup>68</sup></th><th>TiHNH<sub>2</sub>BH<sub>3</sub>-a-<br>McGrady<sup>67</sup></th><th>TiHNH<sub>2</sub>BH<sub>3</sub>-b-<br>McGrady<sup>67</sup></th><th>TiCp<sub>2</sub>NH<sub>2</sub>BH<sub>3</sub>-<br>McGrady<sup>67</sup></th></tr></thead><tbody><tr><td>V(H)</td><td>1.89</td><td>1.92</td><td>1.89</td><td>1.85</td><td>1.92</td></tr><tr><td>M/X/H</td><td>0.07/0.22/1.59</td><td>0.07/0.27/1.58</td><td>0.07/0.24/1.58</td><td>0.07/0.26/1.52</td><td>0.12/0.23/1.57</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.13</td><td>-0.13</td><td>-0.14</td><td>-0.14</td><td>-0.15</td></tr><tr><td>d(M-H<sub>agostic</sub>)</td><td>2.030</td><td>2.039</td><td>1.892</td><td>1.855</td><td>1.892</td></tr><tr><td>d(B-H<sub>agostic</sub>)</td><td>1.302</td><td>1.276</td><td>1.300</td><td>1.310</td><td>1.299</td></tr><tr><td>ω(B-H<sub>agostic</sub>)</td><td>1932</td><td>2084</td><td>1959</td><td>1814</td><td>1895</td></tr><tr><td>d(B-H<sub>free</sub>)</td><td>1.203</td><td>1.203</td><td>1.205</td><td>1.204</td><td>1.204</td></tr><tr><td>ω(B-H<sub>free</sub>)</td><td>2488-2542</td><td>2487-2539</td><td>2474-2526</td><td>2480-2530</td><td>2477-2524</td></tr></tbody></table><p>QTAIM topological parameters for agostic compound: ρ, ∇<sup>2</sup>ρ, H(ρ) in a.u.</p><table><tbody><tr><td>BCP(B-H<sub>agostic</sub>)</td><td>0.122, +0.02, -0.11</td><td>0.132, 0.00, -0.13</td><td>0.123, +0.02, -0.11</td><td>0.122, -0.01, -0.11</td><td>0.126, -0.01, -0.12</td></tr><tr><td>BCP(B-H<sub>free</sub>)</td><td>0.175, -0.27, -0.20</td><td>0.174, -0.26, -0.19</td><td>0.170, -0.21, -0.18</td><td>0.171, -0.21, -0.18</td><td>0.174, -0.26, -0.19</td></tr><tr><td>BCP(M-H<sub>agostic</sub>)</td><td>0.057, +0.11, -0.01</td><td>0.056, +0.11, -0.01</td><td>0.059, +0.11, -0.01</td><td>0.062, +0.13, -0.01</td><td>0.057, +0.16, -0.01</td></tr></tbody></table>

Some relevant parameters for free ligand:

<table><tr><td>d(B-H) (Å)</td><td>1.208</td></tr><tr><td>BCP(B-H<sub>β</sub>)</td><td>0.171, -0.24, -0.19</td></tr></table>

Electron density at BCP(B-H)
in function of d(B-H)

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /><figcaption>**Scheme 1** Rhodium thiophosphoryl pincer studied by Milstein *et al.*, R = but-2-ene or formaldehyde.<sup>66</sup></figcaption></figure>

Tables 10 and 11 report the most relevant topological properties of both Rh-butene-Milstein and Rh-H<sub>2</sub>CO-Milstein complexes.

A close look of the data reported in Tables 10 and 11 allows us to summarize the similarities and differences between the titled complexes as follows:

• As it concerns the weak interactions, despite overall agreement between the topologies of ELF and QTAIM, few differences however have been underlined.<sup>49,75–78</sup> Nevertheless,

we would like to emphasize that we believe in the complementarity of these two methods, rather than mutual exclusion. But having said that, we remind the readers that there is no bond critical point between H<sub>α</sub> and the metallic center, while we have a trisynaptic protonated basin accounting for the α-agostic interaction. Topological analysis of the C–C bonding – bond between the carbon of methyl and that of aryl – obtained from both QTAIM and ELF methods are likewise complementary and often clarify each other. In both Rh-butene-Milstein and Rh-H<sub>2</sub>CO-Milstein pincer complexes the C–C bond valence basin is indeed a trisynaptic basin with a non-negligible contribution from the metallic center (0.06 and 0.09 e<sup>-</sup>) which is a clear indication of the η<sup>3</sup>-C–C–H agostic compound. This conclusion supports the analysis advanced in the paper of Milstein and coworkers.<sup>66</sup> As for the QTAIM analysis, it gives a BCP between the rhodium and the C atom of the aryl group in the case of the Milstein pincer bonded with cis-2-butene, while there exists no such a BCP for the other species. This means that the agostic interaction in the Rh-butene-Milstein could be referred to as the traditional β-agostic compound, while this is not the case for the Rh-H<sub>2</sub>CO-Milstein pincer complex.

• The ligand effect is another striking feature of these pincer-R (R = cis-2-butene or OCH<sub>2</sub>) agostic compounds. It is worth noting that both QTAIM and ELF topologies provide the same analysis for the ligand effect.

<figure><img src="image_19.png" /></figure>
<table><caption>**Table 10** ELF and QTAIM topological features of the agostic bonds for the Rh-H₂CO-Milstein and Rh-butene-Milstein compounds<sup>66</sup></caption><thead><tr><th>H-agostic compound</th><th>QTAIM properties<sup>a</sup></th><th>ELF properties<sup>b</sup></th></tr></thead><tbody><tr><td>Pincer (S-C-S)-R:<br>R = cis-2-butene</td><td>BCP(Hβ, Rh): Yes<br>BCP(Cα, Rh): Yes<br>BCP(Cβ, Rh): No<br>RCP(Rh, Hβ, Cβ, Cα): Yes<br>BCP(C(Ligand), Rh): Yes<br>BCP(C(Ligand), Rh): Yes<br>RCP(Rh, C, C): Yes</td><td>V(Rh, Cβ, Hβ): Yes<br>V(Rh, Cα): No<br>V(Rh, Cβ): No<br>V(Cα, Cβ, Rh): Yes<br>V(C(Ligand), Rh): Yes<br>V(C(Ligand), Rh): Yes</td></tr><tr><td>Pincer (S-C-S)-R:<br>R = OCH<sub>2</sub></td><td>BCP(Hβ, Rh): Yes<br>BCP(Cα, Rh): No<br>BCP(Cβ, Rh): No<br>RCP(Rh, Hβ, Cβ, Cα): No<br>BCP(O(Ligand), Rh): Yes</td><td>V(Rh, Cβ, Hβ): Yes<br>V(Rh, Cα): No<br>V(Rh, Cβ): No<br>V(Cα, Cβ, Rh): Yes<br>V1(O(Ligand), Rh): Yes<br>V2(O(Ligand), Rh): No</td></tr></tbody></table>

<sup>a</sup> BCP and RCP are the bond and ring critical points corresponding to the (3, −1) and (3, +1) critical points. <sup>b</sup> V(X,Y) and V(X,Y,Z) stand for the disynaptic and trisynaptic basins which share two and three core basins, respectively.

* In the case of the cis-2-butene ligand, rhodium thio-phosphoryl pincer cation involves in the formation of two non-equivalent disynaptic basins labeled as V(C(Ligand), Rh) in Tables 10 and 11. The metal atom contribution amounts to 33% and 42% of the total averaged population of these metal-ligand bonds. These basins clearly are indicative of the formation of two metal-carbon coordinate covalent bonds. In parallel, we found two bond critical points for two Rh-C bonds and a RCP in the center of a C-Rh-C triangle. The non-negligible negative values of the energy density at the BCPs (−0.034 and −0.043) clearly indicate the non-negligible covalency of these bonds. As a consequence, the formation of these coordinate covalent bonds enriches the valence shell of the metal leading to lower acidity. Accordingly, the agostic interaction is relatively weak (as shown by values presented in Table 7, with a small covariance V(H)/C(M), and a small contribution of the metal in the valence basin of the agostic hydrogen atom). This weak agostic interaction is geometrically confirmed by a small H<sub>agost</sub>-C distance and a large H<sub>agost</sub>-Rh distance.

* Contrarily, in the case of the formaldehyde-rhodium thiophosphoryl pincer cation, the interaction between oxygen and rhodium atoms is almost weak – manifested by small contribution of Rh in the one of the two valence basins of oxygen – and thus without noticeable change in the valence shell of metal. Thus the electron-deficiency of the Rh center is not counterbalanced by a notable electron transfer. As a consequence, the agostic bonding between Rh and H will be enhanced. Geometrically, the H<sub>agost</sub>-Rh distance will be smaller than in the previous case, and the H<sub>agost</sub>-C distance will be larger, thus leading to a more pronounced activation of the H-C bond.

Thus, the presence of a co-ligand can be of paramount importance in the activation of a C-H bond by means of the formation of an agostic bond, and a topological description of the systems may help in understanding these differences.

On the other hand, a change in the nature of the co-ligand do not necessarily lead to a fundamental change in the agostic character of the bonding. As an example, Table 12 allows to follow the geometry and the topological characterization of the valence basin of the agostic hydrogen atom in β- and γ-agostic alkyl titanium complexes studied by Baird *et al.*<sup>64</sup> The initial agostomers contain two cyclopentadienyl ligands. We suggest the substitution of a cyclopentadienyl ligand either by formaldehyde or by a chlorine atom to not lead to a huge distortion of the agostic bonding: H<sub>agost</sub>-Ti and H<sub>agost</sub>-C distances are almost unchanged (see Table 11). As far as the topological description of the valence basin of the agostic hydrogen atom is concerned, almost no change is observed, neither in the contribution of the metal center to the protonated valence basin nor in the covariance values. Thus, in this case, the substitution of the cyclopentadienyl ligand by formaldehyde or by a chlorine atom does not affect the agostic character of the bond.

To further investigate the parameters influencing the formation of agostic bonds, we studied the influence of the metallic center on the agosticity. The β model compound of Popelier and Logothetis was chosen. Table 13 shows that the substitution of the titanium atom by either Zr or Hf does not affect the geometry of the agostic bond, and the topological description of the valence basin of the agostic hydrogen atom remains similar. Thus, in some cases, and specifically in the

<table><caption>**Table 11** The quantitative topological characteristics of different types of agostic bonds</caption><thead><tr><th>Ago stic compound</th><th>QTAIM properties<sup>a</sup> ρ, ∇<sup>2</sup>ρ, H, ε</th><th>ELF properties<sup>b</sup> M, C and H contributions</th></tr></thead><tbody><tr><td>Pincer (S-C-S)-R:<br>R = cis-2-butene<br>Rh-butene-Milstein<sup>66</sup></td><td>BCP(Hβ, Rh): 0.055; +0.204, −0.008, 0.29<br>BCP(Cα, Rh): 0.050, +0.141, −0.008, 0.12<br>RCP(Rh, Hβ, Cβ, Cα): 0.045, +0.172, −0.004, No ε<br>BCP(C(Ligand), Rh): 0.102, +0.192, −0.034, 1.13<br>BCP(C(Ligand), Rh): 0.111, +0.158, −0.043, 0.37<br>RCP(Rh, C, C): 0.100, +0.297, −0.027, No ε</td><td>V(Rh, Cβ, Hβ): 0.05/1.02/0.95<br>V(Cα, Cβ, Rh): 1.25/0.91/0.06<br>V(C(Ligand), Rh): 0.54/0.23<br>V(C(Ligand), Rh): 0.63/0.21</td></tr><tr><td>Pincer (S-C-S)-R:<br>R = OCH<sub>2</sub><br>Rh-H<sub>2</sub>CO-Milstein<sup>66</sup></td><td>BCP(Hβ, Rh): 0.107, +0.255, −0.046, 0.26<br>BCP(O(Ligand), Rh): 0.086, +0.488, −0.011, 0.46</td><td>V(Rh, Cβ, Hβ): 0.23/0.93/0.96<br>V(Cα, Cβ, Rh): 1.24/0.99/0.09<br>V1(O(Ligand), Rh): 2.31/0.03<br>V2(O(Ligand), Rh): 2.81/0.00</td></tr></tbody></table>

<sup>a</sup> The four QTAIM characteristics at a critical point are given by ρ (the charge density), ∇<sup>2</sup>ρ (the Laplacian of charge density), H (the energy density) and ε (the ellipticity). Note that we have only the first three characteristics at a RCP. These quantities are given in atomic units. <sup>b</sup> The X/Y/Z contributions are the atomic contributions in the averaged population of the V(X,Y,Z) basin. These numbers are in electrons.

<figure><img src="image_20.png" /></figure>
<table><caption>**Table 12** Topological characterization of the β- and γ-agostic alkyl titanium complexes as studied by Baird et al.⁶⁴ and differently substituted</caption><thead><tr><th>Isomer</th><th>Characterization</th><th>β</th><th>γ</th></tr></thead><tbody><tr><td rowspan="4">R<sup>1</sup> = R<sup>2</sup> = Cp</td><td>V(H)</td><td>1.93</td><td>1.95</td></tr><tr><td>Ti/C/H</td><td>0.06/0.77/1.10</td><td>0.05/0.81/1.09</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.08</td><td>-0.07</td></tr><tr><td>d<sup>agostic</sup><sub>Ti-H</sub> (Å)</td><td>2.09</td><td>2.00</td></tr><tr><td>d<sup>agostic</sup><sub>H-C</sub> (Å)</td><td>1.16</td><td>1.15</td></tr><tr><td rowspan="4">R<sup>1</sup> = Cp<br>R<sup>2</sup> = OCH<sub>2</sub><br>(TRANS)</td><td>V(H)</td><td>1.98</td><td>1.99</td></tr><tr><td>Ti/C/H</td><td>0.04/0.86/1.08</td><td>0.03/0.88/1.07</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.07 (AIM: -0.06)</td><td>-0.06 (AIM: -0.05)</td></tr><tr><td>d<sup>agostic</sup><sub>Ti-H</sub> (Å)</td><td>2.12</td><td>2.10</td></tr><tr><td>d<sup>agostic</sup><sub>H-C</sub> (Å)</td><td>1.13</td><td>1.12</td></tr><tr><td rowspan="4">R<sup>1</sup> = R<sup>2</sup> = OCH<sub>2</sub></td><td>V(H)</td><td>Starting from β isomer, it converges to γ agostomer.</td><td>1.99</td></tr><tr><td>Ti/C/H</td><td></td><td>0.06/0.85/1.07</td></tr><tr><td>Cov(V(H)/C(M))</td><td></td><td>-0.09 (AIM: -0.07)</td></tr><tr><td>d<sup>agostic</sup><sub>Ti-H</sub> (Å)</td><td></td><td>2.00</td></tr><tr><td>d<sup>agostic</sup><sub>H-C</sub> (Å)</td><td></td><td>1.13</td></tr><tr><td rowspan="4">R<sup>1</sup> = Cp<br>R<sup>2</sup> = Cl<br>(TRANS)</td><td>V(H)</td><td>1.98</td><td>1.97</td></tr><tr><td>Ti/C/H</td><td>0.08/0.80/1.09</td><td>0.06/0.84/1.08</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.09 (AIM: -0.07)</td><td>-0.08 (AIM: -0.06)</td></tr><tr><td>d<sup>agostic</sup><sub>Ti-H</sub> (Å)</td><td>2.06</td><td>2.01</td></tr><tr><td>d<sup>agostic</sup><sub>H-C</sub> (Å)</td><td>1.14</td><td>1.14</td></tr><tr><td rowspan="4">R<sup>1</sup> = R<sup>2</sup> = Cl</td><td>V(H)</td><td>Starting from β isomer, it converges to γ agostomer.</td><td>1.98</td></tr><tr><td>Ti/C/H</td><td></td><td>0.08/0.84/1.06</td></tr><tr><td>Cov(V(H)/C(M))</td><td></td><td>-0.09 (AIM: -0.07)</td></tr><tr><td>d<sup>agostic</sup><sub>Ti-H</sub> (Å)</td><td></td><td>1.99</td></tr><tr><td>d<sup>agostic</sup><sub>H-C</sub> (Å)</td><td></td><td>1.15</td></tr></tbody></table>

<table><caption>**Table 13** Influence of the substitution of the metallic center in the case of the β model compound of Popelier and Logothetis¹⁸</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Geometry</th><th colspan="2">β model compound (Popelier & Logothetis)</th></tr><tr><th>ELF</th><th>AIM</th></tr></thead><tbody><tr><td>M = Ti</td><td>d(Hβ-Cβ) = 1.152<br>d(M-Cα) = 2.010<br>d(M-Cβ) = 2.393<br>d(M-Hβ) = 2.009<br>a(M-Cα-Cβ) = 84.3<br>a(Cα-Cβ-Hβ) = 113.5</td><td>V(Cβ-Hβ) = 1.93<br>Ti(0.08)/C(0.75)/H(1.10)<br>Cov(V(Hβ)/C(M)) = -0.10</td><td>Q(Hβ) = -0.10<br>Q(Ti) = +1.9<br>Cov(Ti/Hβ) = -0.08</td></tr><tr><td>M = Zr</td><td>d(Hβ-Cβ) = 1.156<br>d(M-Cα) = 2.148<br>d(M-Cβ) = 2.553<br>d(M-Hβ) = 2.149<br>a(M-Cα-Cβ) = 86.3<br>a(Cα-Cβ-Hβ) = 113.9</td><td>V(Cβ-Hβ) = 1.94<br>Zr(0.05)/C(0.74)/H(1.15)<br>Cov(V(Hβ)/C(M)) = -0.10</td><td>Q(Hβ) = -0.15<br>Q(Zr) = +2.6<br>Cov(Zr/Hβ) = -0.08</td></tr><tr><td>M = Hf</td><td>d(Hβ-Cβ) = 1.157<br>d(M-Cα) = 2.140<br>d(M-Cβ) = 2.545<br>d(M-Hβ) = 2.155<br>a(M-Cα-Cβ) = 85.9<br>a(Cα-Cβ-Hβ) = 114.5</td><td>V(Cβ-Hβ) = 1.93<br>Hf(0.04)/C(0.71)/H(1.18)<br>Cov(V(Hβ)/C(M)) = -0.09</td><td>Q(Hβ) = -0.18<br>Q(Hf) = +2.8<br>Cov(Hf/Hβ) = -0.07</td></tr></tbody></table>

<figure><img src="image_21.png" /></figure>
<figure><img src="image_22.png" /></figure>

case of titanocene compounds, the substitution of the metallic center by an atom belonging to the same chemical family, does not affect the agostic character of the bond.

## 6. Characterization of double σ-BH₂ and σ-CH₂ interactions with a metallic center

The term “agostic bonding” was also used in the literature to describe situations in which a small molecule is in interaction with a metallic complex by means of two simultaneous weak interactions.²⁶,⁶⁹,⁷⁰ Sabo-Etienne, Alcaraz *et al.* are particularly active in the study of such intermolecular interactions. To complete our topological study on agostic bonds, we applied our methodology to intermolecular interactions involving dimethylaminoborane⁶⁹ and mesitylborane.⁷⁰

The complexes involved intermolecular interactions involving dimethylaminoborane and studied by Sabo-Etienne *et al.* are presented in Table 14.

In the case of the osmium-containing complexes, an interaction is observed between the σ B-H bond and the metallic center when X = H. The substitution of this hydride by X = Cl dramatically affect the intermolecular interaction. The H atom

trans to the chlorine atom leads to the formation of an hydride and the bond between H₁ and B is broken. On the other hand the σ B-H agostic interaction with the metal is retained for the H₂ atom.

In the case of the ruthenium-containing complexes, an interaction is similarly observed between the σ B-H bond and the metallic center when X = H. As in the case of the osmium-containing complex, the substitution of the X = H atom by a Cl atom causes a distinction between the H₁ and H₂ atoms: the H atom trans to the Cl atom leads to the formation of a stronger interaction with the metal, whereas the interaction between H₂ and B becomes weaker.

For the present study, we selected the osmium-containing complexes for a topological investigation. Indeed, the two osmium-containing complexes allowed us to compare our quantitative approach with the strength of the interaction experimentally observed. Furthermore, these examples gave the opportunity that the topological descriptors herein chosen correctly discriminate a strong σ bond interaction and the formation of a hydride. For the sake of comparison, we also investigate the same system with M = Fe. Table 15 summarizes the topological data obtained for the intermolecular interactions of the four complexes.

In the case of the osmium-containing complexes, the topological criteria selected for the present study are indeed consistent with the formation of an agostic bond, even if the total population of the valence basin of the hydrogen atom is a little bit high for a 3c-2e interaction.

As expected, H₁ and H₂ atoms are equivalent, and the agostic interaction is relatively strong compared with what was expected in the case of intramolecular interactions.

When the X = H atom is substituted by a chlorine atom, the topological description is fully consistent with what is expected from the data available in the literature. In the case of the H₁ atom, the total population of the basin is 1.67 e⁻ and the boron atom is not involved in this basin, which signifies that an

<table><caption>**Table 14** Characterization of the intermolecular interactions between ruthenium- and osmium-containing complexes and dimethylaminoborane investigated by Etienne-Sabo *et al.*⁶⁹</caption><thead><tr><th>Nature of M and X</th><th>Ru, H</th><th>Ru, Cl</th><th>Os, H</th><th>Os, Cl</th></tr></thead><tbody><tr><td>d<sub>H1-M</sub> (Å)</td><td>1.85</td><td>1.62</td><td>1.81</td><td>1.62</td></tr><tr><td>d<sub>H1-B</sub> (Å)</td><td>1.30</td><td>1.55</td><td>1.37</td><td>1.95</td></tr><tr><td>d<sub>H2-M</sub> (Å)</td><td>1.85</td><td>2.00</td><td>1.81</td><td>1.97</td></tr><tr><td>d<sub>H2-B</sub> (Å)</td><td>1.30</td><td>1.27</td><td>1.37</td><td>1.29</td></tr><tr><td>Agostic character</td><td>Agostic</td><td>H₁: strongly agostic<br>H₂: weakly agostic</td><td>Agostic</td><td>H₁: hydride<br>H₂: agostic</td></tr></tbody></table>

<table><caption>**Table 15** Intermolecular interactions involving dimethylaminoborane⁶⁹ and mesitylborane⁷⁰</caption><thead><tr><th>Compound</th><th>Os,H-MeB-Alcaraz<sup>69</sup></th><th>Os,Cl-MeB-Alcaraz<sup>69</sup></th><th>Fe,H-MeB</th><th>Fe,Cl-MeB</th><th>Ru, dimethyl<br>amino-borane-Alcaraz<sup>70</sup></th></tr></thead><tbody><tr><td>V(H)</td><td>2.25</td><td>H₁: 1.67<br>H₂: 2.00</td><td>2.35</td><td>H₁: 2.32<br>H₂: 2.28</td><td>2.16</td></tr><tr><td>M/B/H</td><td>0.18/0.59/1.48</td><td>H₁: 0.39/0.01/1.27<br>H₂: 0.04/0.42/1.54</td><td>0.30/0.54/1.51</td><td>H₁: 0.35/0.46/1.51<br>H₂: 0.22/0.52/1.54</td><td>0.18/0.59/1.37</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.31</td><td>H₁: -0.40<br>H₂: -0.17</td><td>-0.31</td><td>H₁: -0.37<br>H₂: -0.29</td><td>-0.33</td></tr><tr><td>d<sub>H1-M</sub> (Å)</td><td>1.81</td><td>1.62</td><td>1.69</td><td>1.62</td><td>1.87</td></tr><tr><td>d<sub>H1-B</sub> (Å)</td><td>1.37</td><td>1.95</td><td>1.27</td><td>1.29</td><td>1.28</td></tr><tr><td>d<sub>H2-M</sub> (Å)</td><td>1.81</td><td>1.97</td><td>1.69</td><td>1.78</td><td>1.87</td></tr><tr><td>d<sub>H2-B</sub> (Å)</td><td>1.37</td><td>1.29</td><td>1.27</td><td>1.24</td><td>1.28</td></tr></tbody></table>
<figure><img src="image_23.png" /></figure>

hydride is formed, as already reported in the literature⁶⁹ for this case. On the other hand, the topological description of the H₂ atom is consistent with the formation of an agostic interaction. The total population of the valence basin of the hydrogen atom is 2.01 and both Os and B are involved in this basin. The covariance Cov(V(H)/C(M)) as well as the small contribution of the metal in the valence basin of H₂ is clearly consistent with a weaker agostic interaction compared with what is observed for the complex with X = H.

Thus, the topological description herein proposed is qualitatively and quantitatively consistent with what is already reported for these osmium-containing complexes.

When the osmium atom is replaced by Fe, the total population of the valence basin of the hydrogen atoms interacting with the metallic center increases: a total population of 2.35 is calculated when X = H, whereas slightly smaller populations are obtained when X = Cl (2.33 and 2.29 e⁻).

For the “Fe,H-MeB” compound, both H₁ and H₂ atoms are identical, and the contribution of the metallic center to the valence basin of the hydrogen atom is quite large (0.30 e⁻), concomitantly with a large contribution of the boron atom (0.54). Furthermore, the covariance Cov(V(H)/C(M)) is large, thus suggesting that the agostic character of this compound is quite large.

When X = H is substituted by X = Cl, the two H₁ and H₂ atoms are characterized by quantitatively different interactions with the metallic center. Contrary to what was observed in the case of the osmium, the agostic interaction of the two H₁ and H₂ atoms is conserved. The agostic character of the hydrogen atom *trans* to the chloride atom is slightly reinforced, whereas the agostic character of the hydrogen atom *cis* to the chloride atom is slightly weakened.

In the case of a ruthenium-containing dimethylaminoborane compound,⁷⁰ a similar topological description is calculated. Once again, the total population of the valence basin of the hydrogen atom is relatively high, but the contribution of the metallic atom in this basin and the covariance Cov(V(H)/C(M)) are fully consistent with the description of an agostic bond.

During the discussions it is clear that the chlorine substituted compounds lead to the formation of a hydridic bond only when the metallic atom is an osmium atom, whereas a strong agostic interaction is formed with iron and ruthenium atoms.

Thus, from a topological point of view, these intermolecular interactions are similar to the intramolecular, agostic interactions,

and there is no topological reason for not using the same name, “agostic”.

In an attempt to better understand the formation of such intermolecular agostic bonds, model systems were studied:

(1) BH₂Cl forming simultaneously two interactions with Ru(PH₃)Cp,

(2) an isomer of the previous model system in which the BH₂Cl molecule only forms one σ B-H interaction with the metallic center,

(3) an analog of the first model system in which BH₂Cl is substituted by CH₂Cl₂. The formation of such complexes was proposed in the literature⁷⁰,⁷⁹ and was reported with Li as a metallic center,⁸⁰ but not, to our knowledge, in the case of transition metal complexes. On the other hand, cases in which the three hydrogen atoms of a CH₃-R group are simultaneously involved in agostic bonding with a same metallic center, were reported or proposed in the literature.⁸¹,⁸² Thus, a topological description of multiple σ C-H intermolecular agostic bonding is needed.

The topological criteria obtained for these three systems are presented in Table 16.

As already noted in the case of compounds reported in Table 15, relatively large values of total population of the hydrogen basins are observed in the case of the model systems 1 and 3. This corresponds to cases for which the molecule forms simultaneously two σ interactions with the metallic center. In the case of the Model system – 1, the topological criteria suggest the formation of two identical and relatively strong agostic bonds. A similar strong interaction is calculated in the isomeric system forming only one σ B-H interaction with the metallic center (Model system 2). On the other hand, in this case the total population of the valence basin is not particularly high, thus suggesting that the relatively large values of V(H) reported in Tables 15 and 16 are a specific signature of double intermolecular interactions.

In the case of the Model system – 3, both the hydrogen atoms involved in a σ interaction with the metallic center are identical. They are characterized by a lower total population of the valence basin of the hydrogen atoms, in comparison with the Model system – 1. Furthermore, the contribution of the metallic center in these valence basins is particularly weak. This may explain why the experimental formation, the isolation and the characterization of such systems may be difficult.

<table><caption>**Table 16** Model compounds used for the study of intermolecular interactions</caption><thead><tr><th>Compound</th><th>Model system – 1</th><th>Model system – 2</th><th>Model system – 3</th></tr></thead><tbody><tr><td>V(H)</td><td>2.26</td><td>2.07</td><td>2.16</td></tr><tr><td>M/B/H</td><td>0.22/0.56/1.45</td><td>0.19/0.43/1.44</td><td>0.04/01.19/0.93</td></tr><tr><td>Cov(V(H)/C(M))</td><td>-0.41</td><td>-0.35</td><td>-0.14</td></tr><tr><td>d<sup>agostic</sup><sub>H-M</sub> (Å)</td><td>1.73</td><td>1.73</td><td>2.03</td></tr><tr><td>d<sup>agostic</sup><sub>H-B or C</sub> (Å)</td><td>1.34</td><td>1.35</td><td>1.12</td></tr></tbody></table>

<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /></figure>

<figure><img src="image_26.png" /></figure>
# VI. Discussion

## 1. Identification of the existence or non-existence of agostic bonding

To begin with, Scheme 2 summarizes the conditions that should be fulfilled to conclude that an agostic bond exists between the C-H bond and the metallic center.

## 2. A qualitative comparison of different agostic bonds

To further evaluate the capability of the above presented statistical approach to qualitatively characterize the agostic bonding, our theoretical approach will be compared with experimental data. The first necessary step for such a comparison was obviously to find experimental criteria that correctly describe the agostic character of the bonding:

* around different metallic centers,

* involving totally different compounds, and not focused in a well-defined chemical family of compounds,

* for α-, β-, γ-agostomers,

* involving σ C-H and σ B-H agostic bonds, with the possible presence of heteroatoms,

* with constrained geometries in the case of pincer ligands and bi-agostic compounds.

**a. Which kinds of experimental data are available to classify the agostic bonding depending on the strength of the interaction?** Prior to the introduction of the parameters that we chose as a “universal” measurement of the agostic character of the bonding, we would like to briefly overview the experimental approaches currently available for the characterization of such interactions.

The experimental approaches generally used to characterize agostic bonds include NMR shifts, vibrational $\nu(C-H)$ shifts (IR spectroscopy) and X-ray structures.<sup>7,18–20,83</sup> Obviously, the changes in reactivity due to the formation of an agostic bond, and more precisely the activation of the σ bond involved in the interaction with the metallic center, is in itself an evidence that such an interaction exists. More rarely EPR spectroscopy and

visible/UV spectroscopy also proved, in some specific cases, to be able to describe the formation of an agostic bond.<sup>84,85</sup>

These methods will be shortly described below and in Table 17, in the context of the study of agostic bonding.

Experimentally, the formation of an agostic bond is distinguished first and foremost by an activation of the σ bond interacting with the metallic center.<sup>86</sup> In the same chemical family of compounds, it may be possible to qualitatively estimate the agostic character of the bonding, by comparing their reactivity toward a same reagent. On the other hand, such a qualitative approach will be limited to a specific chemical family of compounds and may depend on the reagents. As a conclusion, such an approach, if fully relevant in a screening approach for the most suitable complex in a given reaction process, will not allow the determination of “absolute” agostic character of the bonding.

Spectroscopic methods that were applied to the study of agostic bonding almost cover the whole electromagnetic spectrum, from NMR to X ray spectroscopy.

The comparison of NMR shifts will not allow the comparison of agostic bonds of different chemical natures. Furthermore, only the strongest agostic bonding can be characterized by “classical” one-dimensional NMR experiments. Indeed, all the agostic bonds except the strongest ones are flexible in solution in the timescale of the NMR experiment, thus obscuring the effect.<sup>87,88</sup>

Complementarily to the NMR spectroscopy, EPR spectroscopy was also successfully applied for the characterization of agostic bonds.<sup>85</sup> However, since this method is restricted to the study paramagnetic compounds, EPR spectroscopy is far from an “universal” experimental method that may be used to compare the agostic bonding of a wide set of complexes.

Another spectroscopic approach that was used to experimentally characterize agostic bonding is infrared (IR) spectroscopy.<sup>88</sup> Once again, the characterization of weak interactions in mixtures that may contain agostic and anagostic isomers may be difficult. An elegant procedure to characterize agostic interactions by means of IR spectroscopy, particularly applied by Andrews *et al*.,<sup>89</sup> consists of a coupling with the matrix isolation technique. Species to be analyzed are diluted in an inert medium in the gas phase, and further condensed on surface maintained at cryogenic temperatures. However, such a characterization is far from an easy tool that may be routinely used, without speaking about the experimental difficulty to generate a cryogenic matrix in which an agostic isomer is isolated. On the other hand, vibrational frequencies may be calculated by *ab initio* or DFT procedures once the geometry of the agostic compound is accurately determined.

Complexes of the transition metals are often colored. When the formation of an agostic bonding affects the energies of the d orbitals of the metallic center (especially the HOMO and the LUMO levels), visible/UV spectroscopy may be used to further characterize these interactions. Molybdenum β- and γ-agostomers were recently characterized by this technique.<sup>84</sup> However the characterization of agostic bonds by visible/UV spectroscopy

<figure><img src="image_27.png" /><figcaption>**Scheme 2** Determination of the existence or non-existence of an agostic bonding between a C-H bond and a metallic center.</figcaption></figure>

<figure><img src="image_28.png" /></figure>
<table><caption>**Table 17** Comparison of different experimental approaches that may be used to characterize agostic bonds</caption><thead><tr><th>Method</th><th>Effect leading to a possible characterization of an agostic bond</th><th>General values expected in the case of an agostic bond (and “normal” values expected for all other cases)</th><th>Limitations of the method</th></tr></thead><tbody><tr><td>Reactivity</td><td>Activation of the C-H (or X-H, X = heteroatom such as B) bond.</td><td>Change in the reactivity.<sup>14,86</sup></td><td>The reactivity will obviously depend on co-ligands, geometry and the exact chemical nature of agostic bonding.</td></tr><tr><td>NMR</td><td>Redistribution of bonding-electron density in the formation of an agostic bonding.</td><td>δ = -5 to -15 ppm, <sup>1</sup>J<sub>C-H,agostic</sub> = 75 to 100 Hz for C-H agostic bonding (<sup>1</sup>J<sub>Csp<sup>3</sup>-H<sub>free</sub></sub> = 128 MHz)<sup>7,18,83,91</sup><br>Downfield paramagnetic shifts in the 700-1100 ppm range for axial ligands.<sup>92</sup></td><td>NMR spectra cannot be obtained for dynamical systems.</td></tr><tr><td>EPR spectroscopy</td><td>Influence on the Zeeman electronic effect of the distortion of the geometry.</td><td>Variations in the values of g's and determination of ΔH<sup>o</sup> and ΔS<sup>o</sup>.<sup>85,93</sup></td><td>Limited to paramagnetic compounds for with the single electron is involved in the agostic bonding.</td></tr><tr><td>IR spectroscopy</td><td>Weakening of the C-H bond leading to a reduction in frequency for ν(C-H) vibrational mode.</td><td>ν(C-H) = 2300-2700 cm<sup>-1</sup> (to be compared with the ≈ 2700-3000 cm<sup>-1</sup> range for free ligands).<sup>18-20</sup></td><td>Difficulty to identify a small weak signal that may overlap with other signals.</td></tr><tr><td>Visible/UV spectroscopy</td><td>Valence electron transition that may be affected by the formation of the M···H agostic bonding.</td><td>Depends on the crystal field perturbation induced by the agostic bonding.<sup>84</sup></td><td>Limited to cases for which the agostic bondings affect the color of the compound</td></tr><tr><td>X-ray diffraction</td><td>Geometrical distortion.</td><td>Distance H<sub>agostic</sub>-M = 1.8-2.3 Å (2.3-2.9 Å).<sup>7,94</sup><br>Distance H<sub>agostic</sub>-X larger than the corresponding value for the free ligand<br>Change in the angles between the atoms.<sup>7</sup></td><td>Difficulty to localize hydrogen atoms.<br>Limited to agostomers that may be isolated as crystals.</td></tr></tbody></table>

can only be applied in very specific cases for which the interactions between the metallic center and the σ bond affect the color of the species.

Since an agostic bonding should lead to a distortion of the geometry of a compound, the crystalline X-ray structure is an obvious characterization tool.<sup>90</sup> The low scattering factor for hydrogen atoms and the difficulty to localize hydrogen in the vicinity of a metallic center are some of the limits of the method, in addition to the fact that agostic bonds may exist as intermediate species in liquid solution during reaction processes. In these latter cases, agostomers may not be obtained in the crystalline form. Even if limited to agostomers that may be isolated as crystals, and despite to its limitations X-ray diffraction structures were often reported in the literature for compounds involving agostic bonds. Furthermore, other experimental techniques provide data that may be further used for theoretical studies of the species. Thus geometries of agostomers can be obtained from DFT calculations and post-Hartree-Fock approaches. Theoretical calculations also allow the determination of spectroscopic data that may be compared with experimental values. It is then possible to determine structures for these compounds. From a geometrical point of view, the formation of an agostic bond will affect:

* the angles between the atoms,

* the H<sub>agostic</sub>-M distance,

* the H<sub>agostic</sub>-X distance.

**b. Choice of representative experimental criteria that may describe the strength of agostic bonding.** Among above-mentioned characterizations, only vibrational frequencies and geometrical parameters can easily be obtained for all the types

of agostic bonds, either directly experimentally, or by means of a combined experimental and theoretical investigation.

The (harmonic) vibrational shift of the H-X stretching mode, Δω(H-X), caused by the agostic bonding, is an obvious tool to investigate the strength of the interaction. This shift is defined as follows:

$$ \Delta\omega(\text{H-X}) = \omega(\text{H-X})^{\text{free ligand}} - \omega(\text{H-X})^{\text{agostic compound}} $$

The shift may depend on the nature of the X atom. In order to compare agostic bonds involving chemically different X atoms, we suggest the use of a normalized parameter:

$$ \frac{\Delta\omega(\text{H-X})}{\omega(\text{H-X})^{\text{free ligand}}} $$

Additionally, geometrical parameters such as angles and distances were already proposed to characterize agostic bonds.<sup>7</sup> The angles between the atoms will obviously depend on the system, and will not be directly linked to the agostic character of the bonding.

As already mentioned, the values of the angles cannot be used to compare the agostic character of bondings, in the set of compounds chosen for the present study, because their initial geometries are too much different from each other (Fig. 1). The H<sub>agostic</sub>-M distance will decrease, concomitantly with the increase in the H<sub>agostic</sub>-M distance, during the formation of an agostic bond. On the other hand, these two distances will be strongly affected by the nature of the X atom. In a first
approximation, we can suggest that the ratio between $d_{H-X}^{agostic}$ and $d_{H-M}^{agostic}$, $\frac{d_{H-X}}{d_{H-M}}$ will not strongly depend on the nature of X.

Thus, two normalized experimental parameters may be used to compare the strength of agostic bonding, $\frac{d_{H-X}}{d_{H-M}}$ and $\frac{\Delta\omega(H-X)}{\omega(H-X)^{free ligand}}$. A simpler way to compare agostic bonding in a same family of complexes involving only one chemical type of H-X as agostic bond is to directly use the $d_{H-X}^{agostic}$ and $\Delta\omega(H-X)$ experimental parameters.

**c. May statistical parameters obtained from the ELF/QTAIM study be used to qualitatively estimate the agostic bonding strength?** Once it was checked that the criteria summarized in Scheme 2 are consistent with the existence of an agostic bond, several statistical parameters appear to be particularly relevant to qualitatively compare the strength of agostic bonding.

First of all, the contribution of the metallic center in the protonated agostic V(H) basin, $M_{V(H)}$, is an obvious important parameter.

On the other hand, this parameter depends on the whole structure of the metallic complex. As a consequence, a direct comparison of the $M_{V(H)}$ values will not allow the rigorous quantification of the agostic character of bonds. To overcome this limit, we propose to introduce a normalized parameter, the normalized metal contribution in the V(H) basin ($N_{M_{V(H)}}$). Its value is calculated by the ratio between the metal contribution in the V(H) basin and the total population of the V(H) basin:

$$N_{M_{V(H)}} = \frac{M_{V(H)}}{V(H)} \times 100 \quad (13)$$

The main advantage of this normalized parameter toward the $M_{V(H)}$ value is to take into account the variation of the number of electrons in the V(H) basin that may be induced by other atoms in the vicinity of H.

Fig. 4 shows the variations of $M_{V(H)}$ and $N_{M_{V(H)}}$ as a function of the $(\frac{d_{H-X}}{d_{H-M}})$ ratio. The last point in the plots of the metal contribution corresponds to the hydride atom. This point was not included in the curve of the covariance. For all the other studied systems, all the values were taken into account in the plots.

As already pointed out, the complexes selected for the present study cover a wide range of chemical families, owing to the nature of the metallic center, the nature of the ligand involved in the agostic bonding, the nature of co-ligands and the geometry of the complexes. Because of these chemical differences, a relative spread of the of the $M_{V(H)}$ and $N_{M_{V(H)}}$ was expected. Given this point, a satisfactory linear correlation is observed for $M_{V(H)}$ and $N_{M_{V(H)}}$ with the ratio of the distances. The coefficient of the linear regression clearly confirm that $N_{M_{V(H)}}$ is a suitable criterion to compare agostic bonds in compounds belonging to different chemical families. Furthermore, the linear dependence observed demonstrates that $N_{M_{V(H)}}$ can indeed be used to qualitatively characterize the agostic character of the bonding.

<figure><img src="image_29.png" /><figcaption>Fig. 4 Variation of the $M_{V(H)}$ and $N_{M_{V(H)}}$ parameters as a function of the $\frac{d_{H-X}}{d_{H-M}}$ ratio, with X = B or C, and $N_{M_{V(H)}} = \frac{M_{V(H)}}{V(H)} \times 100$.</figcaption></figure>

As mentioned earlier, the electron density at the bond critical point of the X-H bond is also an indication of the σ-donation to the metallic center and thus of the strength of the agostic bonding. Once again, if directly used, this parameter will depend on the system under investigation. As a consequence, the direct comparison of the $\rho(BCP_{X-H})$ may be a suitable way to classify the strength of agostic bonds for a series of complexes belonging to a same chemical family of compounds, whereas the definition of a normalized value may be useful to compare agostic bonding in a set of complexes belonging to different chemical families. To this extent, we suggest to use the ratio between the $\rho(BCP)$ for the X-H bond in the agostic complex and the $\rho(BCP)$ for the X-H bond in the free ligand:

$$\frac{\rho(BCP(X-H_{agostic}))}{\rho(BCP(X-H_{free}))}$$

Fig. 5 illustrates the possibility of using:

* direct H-X distances to characterize the strength of agostic bonds inside a same chemical family of compounds;

* reduced frequencies of the H-X stretching mode to characterize the strength of agostic bonds bonding in a set of complexes belonging to different chemical families;

* reduced $\rho(BCP)$ for the H-X bond to classify the interaction with a metallic center.

Finally, the X-H bonds can be divided into four categories, depending on the strength of the M-H agostic bonding probably formed with a metallic center. Once the existence of the agostic bonding is ascertained based on the criteria presented in Scheme 2, four main estimators can be used to evaluate the strength of the interaction:

* two experimental parameters, $\frac{d_{H-X}}{d_{H-M}}$ and $\frac{\Delta\omega(H-X)}{\omega(H-X)^{free ligand}}$; alternatively these parameters may be obtained by geometry optimization and theoretical calculations;

* the $N_{M_{V(H)}}$ normalized parameter, based on ELF/QTAIM calculations;

<figure><img src="image_30.png" /></figure>
<figure><img src="image_31.png" /><figcaption>Fig. 5 Variation of the reduced $\frac{\rho(\text{BCP}(X-H_{\text{agostic}}))}{\rho(\text{BCP}(X-H_{\text{free}}))}$ QTAIM parameter as a function of the: (A) $d_{H-X}$ distance, and (B) reduced frequencies of the H-X stretching mode. The strength of the agostic bonding increases from left to right and from top to bottom in both the graphs.</figcaption></figure>

<table><caption>Table 18 A qualitative classification of the strength of agostic bonding based on three criteria</caption><thead><tr><th rowspan="2">Criterion</th><th colspan="4">Characterization of the M-H interaction</th></tr><tr><th>Anagostic</th><th>Weak-medium agostic</th><th>Medium-strong agostic</th><th>Pre-hydride</th></tr></thead><tbody><tr><td>d<sub>H-X</sub></td><td>&lt;0.5</td><td>0.5-0.7</td><td>0.65-0.8</td><td>&gt;1</td></tr><tr><td>d<sub>H-M</sub></td><td>&ge;0.90</td><td>0.90-0.70</td><td></td><td>&lt;0.65</td></tr><tr><td>&Delta;&omega;(H-X)/&omega;(HX)<sup>free ligand</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>N<sub>M(V(H))</sub></td><td>&lt;1</td><td>1-5</td><td>5-15</td><td>&gt;20</td></tr><tr><td>&rho;(BCP(X-H<sub>agostic</sub>))/&rho;(BCP(X-H<sub>free</sub>))</td><td>&ge;0.90</td><td>0.90-0.70</td><td></td><td>&lt;0.65</td></tr></tbody></table>

• the $\frac{\rho(\text{BCP}(X-H_{\text{agostic}}))}{\rho(\text{BCP}(X-H_{\text{free}}))}$ normalized QTAIM parameter.

In Table 18, a global classification is proposed to estimate the strength of an agostic bonding, based on these four parameters.

## VII. Conclusion

Herein we presented a topological quantitative investigation of σ C-H interactions. It was shown that the agostic character of the bonding can be compared even in the case of

• complexes around a different metallic center,

• complexes containing different ligands and co-ligands,

• α-, β-, γ-, and η<sup>3</sup>-CCH agostomers,

• mono- and bis-agostic species,

• inter- and intramolecular agostic bonding.

Thus the present method makes it possible to compare the agostic character of interactions that take place in systems that do not belong to a same chemical family, by means of a simple topological approach. The large set of compounds presented in this article shows that the methodology proposed here can lead to an unambiguous determination of the agostic character, both qualitatively and quantitatively of a wide range of chemical compounds.

Furthermore, the topological description of intramolecular σ C-H and σ B-H interactions, as well as intermolecular σ B-H interactions, is similar, thus suggesting that the use of the same “agostic” term for all these interactions, is appropriate.

It has been evidenced that any X-H···M agostic bond could be fully analyzed using three topological descriptors obtained easily with the ELF method within the free TopMod code.

We have to pay a particular attention on the stability and quality of the wave function derived from an optimization procedure. Indeed, the ELF topological analysis is actually an *a posteriori* approach whose reliability depends on that of the wave function in hand.

The most relevant descriptors are:

1. The existence of a protonated trisynaptic basin, labeled as V(H<sub>ag</sub>),

2. Three contributors participate in the averaged population of this basin: M/X/H,

3. The first and second descriptors are actually similar. We suggest them explicitly in order to underline the trisynaptic character of this basin which accounts for the 3c-2e interactions. Nevertheless, we also suggest the use of dimensionless descriptors to qualitatively estimate the strength of agostic bonding:

• the ratio between the metal contribution and the total population of the protonated basin, in order to emphasize the importance of the relative metal contribution. This $N_{M(V(H))} =$

$$ \frac{M_{V(H)}}{V(H)} \times 100 $$ ratio allows the quantitative classification of the agostic character of bonding;

• the ratio of the electron densities at the BCP, $\frac{\rho(\text{BCP}(X-H_{\text{agostic}}))}{\rho(\text{BCP}(X-H_{\text{free}}))}$ in order to take into account the weakening of the X-H bond due to the interaction with the metallic center.

<figure><img src="image_32.png" /></figure>
4. We suggest the use of the above-mentioned descriptors in combination with two other criteria, ideally derived from experimental data. Indeed, experimental or, if not available, theoretical data, may be used to determine two additional reduced parameters that suitably describe the strength of agostic bonds:

* the weakening of the H-X bond due to the agostic interaction is characterized by a shift in the frequencies of the vibrational stretching mode of the H-X bond. The $\frac{\Delta\omega(H-X)}{\omega(H-X)^{\text{free ligand}}}$ reduced frequency parameter may be obtained from experimental IR spectra, or from (harmonic) calculations;

* the $\frac{d_{H-X}}{d_{H-M}}$ reduced distance may be obtained or derived from X-ray structures. This reduced distance is increasing within the strength of the agostic bonding.

5. The covariance between two basins: the agostic protonated basin and the metallic core basin C(M). This quantity provide a measure of association between two quantities: C(M) and V(H<sub>ag</sub>). Agostic bonding for which the ratio N<sub>M<sub>v(H</sub>)</sub> is similar, the agostic character of the bonds increase with the |C(M)|.

Finally, we would like to emphasize the mutual complementary aspect of the topological results obtained from QAIM and ELF analysis.

## VIII. Computational details

The ELF calculations have been done using TopMod software.⁹⁵ Furthermore, the AIMAll software was used for the quantitative study of the topological QAIM data.⁹⁶

All the necessary wfn files for the topological investigations have been obtained using the Gaussian 09 Rev D.01 quantum chemical package. Calculations have been performed using different density functionals (B3LYP, PBE0 and TPSSh) as well as several basis sets (6-31++G(2d,2p), 6-311++G(2d,2p), TZVP, and LanL2DZ) as implemented in the Gaussian 09 package. In addition, the Def2-TZVP, Def2-TZVPD, LanL2TZ(P), and LanL2TZ(F) basis sets of Ahlrichs et al.,⁹⁷ obtained from "EMSL Basis Set Exchange Library".⁹⁸ Furthermore, in order to take into account the dispersion contribution in some model compounds (mono- and bi-agostic model structures), geometry optimization have been done using the wB97XD range separated hybrid functional⁹⁹ (Table 16) and also using the B2PLYPD3 double hybrid functional including the D3 version of Grimme's dispersion with the original D3 damping function.¹⁰⁰ In order to ensure that an optimized structure corresponds to a true minimum, a frequency analysis was performed. Concerning the geometrical structures of the studied complexes, additional data on the optimized geometries are available in the ESI.†

## Note added in proof

Since this work was submitted one paper dealing with the agostic-bonding in a particular set of pincer complexes has appeared.¹⁰¹ In this work several topological approaches have been used. Conclusions fully agree our description.

## References

1. J. Chatt and J. M. J. Davidson, The tautomerism of arene and ditertiary phosphine complexes of ruthenium(0), and the preparation of new types of hydrido-complexes of ruthenium(II), *J. Chem. Soc.*, 1965, 843–855.

2. A. H. Janowicz and R. G. Bergman, Carbon–hydrogen activation in saturated hydrocarbons: direct observation of M + R–H → M(R)(H), *J. Chem. Soc.*, 1982, **104**, 352–354.

3. W. D. Jones, A. J. Vetter, D. D. Wick and T. O. Northcutt, Alkane Complexes as Intermediates in C–H Bond Activation Reactions, in *Activation and Functionalization of C–H Bonds*, ed. K. I. Goldberg and A. S. Goldman, ACS Symposium Series 885, Washington, DC, 2004, ch. 3, pp. 56–69.

4. M. Brookhart, M. L. H. Green and R. B. A. Pardy, Two-electron, three-centre carbon–hydrogen–cobalt bonds in the compounds [Co(η–C₅Me₄R)(η–C₂H₄)(η–C₂H₄·μ–H)]BF₄, R = Me and Et, *J. Chem. Soc., Chem. Commun.*, 1983, 691–693.

5. M. Brookhart and M. L. H. Green, Carbon-hydrogen-transition metal bonds, *J. Organomet. Chem.*, 1983, **250**, 395–408.

6. M. Brookhart, M. L. H. Green and L. L. Wong, Carbon-Hydrogen-Transition Metal Bonds, *Prog. Inorg. Chem.*, 1988, **36**, 1–124.

7. M. Brookhart, M. L. H. Green and G. Parkin, Agostic interactions in transition metal compounds, *Proc. Natl. Acad. Sci. U. S. A.*, 2007, **104**, 6908–6914.

8. R. H. Crabtree, Organometallic alkane CH activation, *J. Organomet. Chem.*, 2004, **689**, 4083–4091.

9. B. Rybtchinsk, L. Konstantinovsky, L. J. W. Shimon, A. Vigalok and D. Milstein, Solvent-Stabilized Alkylrhodium(III) Hydride Complexes: A Special Mode of Reversible C–H Bond Elimination Involving an Agostic Intermediate, *Chem. – Eur. J.*, 2000, **6**, 3287–3291.

10. I. Omae, Agostic bonds in cyclometallation, *J. Org. Chem.*, 2011, **696**, 1128–1145.

11. E. Clot, O. Eisenstein, T. Dubé, J. W. Faller and R. H. Crabtree, Interplay of Weak Interactions: An Iridium(III) System with an Agostic tert-Butyl but a Nonagostic Isopropyl Group, *Organometallics*, 2002, **21**, 575–580.

12. T. Ziegler, The 1994 Alcan Award Lecture Density functional theory as a practical tool in studies of organometallic energetics and kinetics. Beating the heavy metal blues with DFT, *Can. J. Chem.*, 1995, **73**, 743–761.

13. H. H. Brintzinger, D. Fischer, R. MRLhaupt, B. Rieger and R. M. Waymouth, Stereospecific Olefin Polymerization with Chiral Metallocene Catalysts, *Angew. Chem., Int. Ed. Engl.*, 1995, **34**, 1143–1170.

14. M. Albrecht, Cyclometalation Using d-Block Transition Metals: Fundamental Aspects and Recent Trends, *Chem. Rev.*, 2010, **110**, 576–623.

15. J. Mathew, N. Koga and C. H. Suresh, C–H Bond Activation through σ-Bond Metathesis and Agostic Interactions: Deactivation Pathway of a Grubbs Second-Generation Catalyst, *Organometallics*, 2008, **27**, 4666–4670.

16. M. Lein, J. A. Harison and A. J. Nielson, Identification of non-classical C–H···M interactions in early and late

<figure><img src="image_33.png" /></figure>
transition metal complexes containing the CH(ArO)₃ ligand, *Dalton Trans.*, 2013, **42**, 10939-10951.

17 F. Fuster and B. Silvi, Does the topological approach characterize the hydrogen bond?, *Theor. Chem. Acc.*, 2000, **104**, 13-21.

18 P. L. A. Popelier and G. Logothetis, Characterization of an agostic bond on the basis of the electron density, *J. Organomet. Chem.*, 1998, **555**, 101-111.

19 W. Scherer and G. S. McGrady, Agostic Interactions in d⁰ metal alkyl complexes, *Angew. Chem., Int. Ed.*, 2004, **43**, 1782-1806.

20 M. Lein, Characterization of agostic interactions in theory and computation, *Coord. Chem. Rev.*, 2009, **253**, 625-634.

21 M. Etienne and A. S. Weller, Intramolecular C-C agostic complexes: C-C sigma interactions by another name, *Chem. Soc. Rev.*, 2014, **43**, 242-259.

22 W. H. Bernskoetter, C. K. Schauer, K. I. Goldberg and M. Brookhart, Characterization of a rhodium(I) σ-methane complex in solution, *Science*, 2009, **326**, 553-556.

23 J. A. Calladine, S. B. Duckett, M. W. George, S. L. Matthews, R. N. Perutz, O. Torres and K. Q. Vuong, Manganese alkane complexes: an IR and NMR spectroscopic investigation, *J. Am. Chem. Soc.*, 2011, **133**, 2303-2310.

24 S. Lachaize and S. Sabo-Etienne, σ-Silane Ruthenium Complexes: The Crucial Role of Secondary Interactions, *Eur. J. Inorg. Chem.*, 2006, 2115-2127.

25 G. Alcaraz and S. Sabo-Etienne, NMR: A good tool to ascertain σ-silane or σ-borane formulations?, *Coord. Chem. Rev.*, 2008, **252**, 2395-2409.

26 J. F. Hartwig, C. N. Muhoro and X. He, Catecholborane bound to titanocene. Unusual coordination of ligand σ bonds, *J. Am. Chem. Soc.*, 1996, **118**, 10936-10937.

27 A. T. Çolak, O. Z. Yeşilel and O. Büyükgüngör, A novel anagostic C single bond H···Cu interaction and clustered water molecules in copper(II)-pyridine-2,5-dicarboxylate complex, *J. Mol. Struct.*, 2011, **991**, 68-72.

28 H. V. Huynh, L. R. Wong and P. S. Ng, Anagostic Interactions and Catalytic Activities of Sterically Bulky Benzannulated N-Heterocyclic Carbene Complexes of Nickel(II), *Organometallics*, 2008, 27, 2231-2237.

29 M. G. D. Holaday, G. Tarafdar, A. Kumar, M. L. P. Reddy and A. Srinivasan, Exploring anagostic interactions in 5,15-porphodimethene metal complexes, *Dalton Trans.*, 2014, **43**, 7699-7703.

30 M. L. Yang, Y. A. Zhu, C. Fan, Z. J. Sui, D. Chen and X. G. Zhou, Density functional study of the chemisorption of C1, C2 and C3 intermediates in propane dissocuation on Pt(111), *J. Mol. Catal. A: Chem.*, 2010, **321**, 42-49.

31 G. Papoian, J. K. Nørskov and R. Hoffmann, A comparative theoretical study of the hydrogen, methyl, and ethyl chemisorption on the Pt(111) surface, *J. Am. Chem. Soc.*, 2000, **122**, 4129-4144.

32 T. Jacob and W. A. Goddard III, Agostic Interactions and Dissociation in the First Layer of Water on Pt(111), *J. Am. Chem. Soc.*, 2004, **126**, 9360-9368.

33 H. G. Cho, Agostic Structure of Titanium Methylidene Hydride (CH₂=TiH₂): A Theoretical Investigation for the

Elusive Intramolecular Interaction, *Bull. Korean Chem. Soc.*, 2009, **30**, 1631-1633.

34 S. Niu and M. B. Hall, Theoretical Studies on Reactions of Transition-Metal Complexes, *Chem. Rev.*, 2000, **100**, 353-406.

35 M. D. Su and S. Y. Chu, Theoretical Model for Insertion of the 16-Electron Species (η⁵-C₅H₅)M(L) into Saturated Hydrocarbons. A (η⁵-C₅H₅)M(CO) + CH₄ (M = Ru-, Os-, Rh, Ir, Pd+, Pt+), *Organometallics*, 1997, **16**, 1621-1627.

36 S. Zaric and M. B. Hall, Ab Initio Calculations of the Geometries and Bonding Energies of Alkane and Fluoroalkane Complexes with Tungsten Pentacarbonyl, *J. Phys. Chem. A*, 1997, **101**, 4646-4652.

37 B. O. Roos, R. Lindh, H. G. Cho and L. Andrews, Agostic Interaction in the Methylidene Metal Dihydride Complexes H₂MCH₂ (M = Y, Zr, Nb, Mo, Ru, Th, or U), *J. Phys. Chem. A*, 2007, **111**, 6420-6424.

38 E. D. Glendening, C. R. Landis and F. Weinhold, Natural bond orbital methods, *Wiley Interdiscip. Rev.: Comput. Mol. Sci.*, 2012, 2, 1-42.

39 F. Weinhold and C. R. Landis, Natural bond orbitals and extensions of localized bonding concepts, *Chem. Educ.: Res. Pract. Eur.*, 2001, 2, 91-104.

40 F. Weinhold and C. R. Landis, *Valency and Bonding – a natural bond donor-acceptor perspective*, Cambridge University Press, 2005.

41 H. G. Cho and B. S. Cheong, Agostic Interaction of the Smallest Zirconium Methylidene Hydride: Reproduction of the Distorted Structure Experimentally Observed in Matrix Infrared Spectra, *Bull. Korean Chem. Soc.*, 2009, **30**, 470-481.

42 T. S. Thakur and G. R. Desiraju, Theoretical investigation of C–H···M interactions in organometallic complexes: A natural bond orbital (NBO) study, *THEOCHEM*, 2007, **810**, 143-154.

43 M. Kaupp, The Unusual Fluxional Structure of Tetramethyloxotungsten: Quantum Chemical Structure Predictions for the d0 and d1 Complexes [MOR₄] (M = W, Re; R = H, CH₃), *Chem. – Eur. J.*, 1998, **4**, 2059-2071.

44 R. F. Bader, *Atoms in Molecules: A Quantum Theory*, Oxford University Press, Oxford, 1994.

45 S. Pillet, G. Wu, V. Kulsomphob, B. G. Harvey, R. D. Ernst and P. Coppens, Investigation of Zr-C, Zr-N, and Potential Agostic Interactions in an Organozirconium Complex by Experimental Electron Density Analysis, *J. Am. Chem. Soc.*, 2003, **125**, 1937-1949.

46 J. Saßmannshausen, Agostic or not? Detailed Density Functional Theory studies of the compounds, *Dalton Trans.*, 2011, **40**, 136-141.

47 W. Scherer, V. Herz, A. Brück, C. Hau, F. Reiner, S. Altmannshofer, D. Leusser and D. Stalke, The Nature of β-Agostic Bonding in Late-Transition-Metal Alkyl Complexes, *Angew. Chem., Int. Ed.*, 2011, **50**, 2845-2849.

48 R. F. W. Bader and C. F. Matta, Bonding to titanium, *Inorg. Chem.*, 2001, **40**, 5603-5611.

49 V. Tognetti, L. Joubert, R. Raucoules, T. de Bruin and C. Adamo, Characterizing agosticity using the Quantum

<figure><img src="image_34.png" /></figure>
Theory of Atoms in Molecules: Bond Critical Poins and their local properties, *J. Phys. Chem. A*, 2012, **116**, 5472–5479.

50 A. D. Becke and K. E. Edgecombe, A simple measure of electron localization in atomic and molecular system, *J. Chem. Phys.*, 1990, **92**, 5397–5403.

51 B. Silvi and A. Savin, Classification of chemical bonds based on topological analysis of electron localization functions, *Nature*, 1994, **371**, 683–686.

52 A. Savin, B. Silvi and F. Colonna, Topological analysis of the electron localization function applied to delocalized bonds, *Can. J. Chem.*, 1996, **74**, 1088–1096.

53 B. Silvi, I. Fourré and M. E. Alikhani, The topological analysis of the electron localization function. A key for a position space representation of chemical bonds, *Monatsh. Chem.*, 2005, **136**, 855–858.

54 P. W. Ayers, Electron localization functions and local measures of the covariance, *J. Chem. Sci.*, 2005, **117**, 441–454.

55 S. Diner and P. Claverie, in *Localization and Delocalization in Quantum Chemistry*, ed. O. Chalvet, R. Daudel, S. Diner and R. Malrieu, Dordrecht, 1976, vol. II, pp. 395–448.

56 B. Silvi, The Spin-Pair Compositions as Local Indicators of the Nature of the Bonding, *J. Phys. Chem. A*, 2003, **107**, 3081–3085.

57 S. Raub and G. Jansen, A quantitative measure of bond polarity from the electron localization function and the theory of atoms in molecules, *Theor. Chem. Acc.*, 2001, **106**, 223–232.

58 F. Fuster and S. Grabowski, Intramolecular hydrogen bonds: the QTAIM and ELF characteristics, *J. Phys. Chem. A*, 2011, **115**, 10078–10086.

59 I. Fourré and J. Bergès, Structural and topological characterization of the three-electron bond: the SO radical, *J. Phys. Chem. A*, 2004, **108**, 898–906.

60 I. Fourré and B. Silvi, What can we learn from two-center three-electron bonding with the topological analysis of ELF?, *Heteroat. Chem.*, 2007, **18**, 135–160.

61 J. Poater, M. Duran, M. Solà and B. Silvi, Theoretical evaluation of electron delocalization in aromatic molecules by means of Atoms in Molecules (AIM) and Electron Localization Function (ELF) topological approaches, *Chem. Rev.*, 2005, **105**, 3911–3947.

62 K. B. Wiberg, R. F. W. Bader and C. D. H. Lau, Theoretical analysis of hydrocarbon properties. 2. Additivity of group properties and the origin of strain energy, *J. Am. Chem. Soc.*, 1987, **109**, 1001–1012.

63 A. F. Dunlop-Brière, P. H. M. Budzelaar and M. C. Baird, α- and β-agostic alkyl-titanocene complexes, *Organometallics*, 2012, **31**, 1591–1594.

64 A. F. Dunlop-Brière, M. C. Baird and P. H. M. Budzelaar, [Cp2TiCH2CHMe(SiMe3)]+, an Alkyl–Titanium Complex Which (a) Exists in Equilibrium between a β-Agostic and a Lower Energy γ-Agostic Isomer and (b) Undergoes Hydrogen Atom Exchange between α-, β-, and γ-Sites via a Combination of Conventional β-Hydrogen Elimination–Reinsertion and a Nonconventional CH Bond Activation Process Which Involves Proton Tunnelling, *J. Am. Chem. Soc.*, 2013, **135**, 17514–17527.

65 T. A. Atesin, T. Li, S. Lachaize, W. W. Brennessel, J. J. Garcia and W. D. Jones, Experimental and theoretical examination of C-CN and C-H bond activation of acetonitrile using zerovalent nickel, *J. Am. Chem. Soc.*, 2007, **129**, 7562–7569.

66 M. Montag, I. Efremenko, Y. Diskin-Posner, Y. Ben-David, J. M. L. Martin and D. Milstein, Exclusive C-C oxidative addition in a rhodium thiophosphoryl pincer complex and computational evidence for an η³-C-C-H agostic intermediate, *Organometallics*, 2012, **31**, 505–512.

67 D. J. Wolstenholme, K. T. Traboulsee, A. Decken and G. S. McGrady, Structure and Bonding of Titanocene Amido-borane Complexes: A Common Bonding Motif with Their β-Agostic Organometallic Counterparts, *Organometallics*, 2010, **29**, 5769–5772.

68 T. D. Forster, H. M. Tuononen, M. Parvez and R. Roesler, Characterization of β-B-Agostic Isomers in Zirconocene Amido-boraneComplexes, *J. Am. Chem. Soc.*, 2009, **131**, 6689–6691.

69 G. Bénac-Lestrille, U. Helmstedt, L. Vendier, G. Alcaraz and E. Clot, S. Sabo-Etienne. Dimethylaminoborane coordination to late transition metal centers: Snapshots of the B-H oxidative addition process, *Inorg. Chem.*, 2011, **50**, 11039–11045.

70 G. Alcaraz, E. Clot, S. Sabo-Etienne, L. Vendier and U. Helmstedt, Mesitylborane as a bis (σ-B-H) ligand: an unprecedented bonding mode to a metal center, *J. Am. Chem. Soc.*, 2007, **129**, 8704–8705.

71 C. Adamo and V. Barone, Toward reliable density functional methods without adjustable parameters: The PBE0 model, *J. Chem. Phys.*, 1999, **110**, 6158–6169.

72 J. Tao, J. P. Perdew, V. N. Staroverov and G. E. Scuseria, Climbing the density functional ladder: nonempirical meta–generalized gradient approximation designed for molecules and solids, *Phys. Rev. Lett.*, 2003, **146401**, 1–4.

73 W. Scherer, D. J. Wolstenholme, V. Herz, G. Eickerling, A. Brück, P. Benndorf and P. W. Roesky, On the nature of agostic interaction in transition-metal amido complexes, *Angew. Chem., Int. Ed.*, 2010, **49**, 2242–2246.

74 U. Koch and P. L. A. Popelier, Characterization of C-H-O Hydrogen Bonds on the Basis of the Charge Density, *J. Phys. Chem.*, 1995, **99**, 9747–9754.

75 I. Vidal, S. Melchor, I. Alkorta, J. E Iguero, M. R. Sundberg and J. A. Dobado, On the Existence of r-Agostic Bonds: Bonding Analyses of Titanium Alkyl Complexes, *Organometallics*, 2006, **25**, 5638–5647.

76 B. Silvi, The synaptic order: a key concept to understand multicenter bonding, *J. Mol. Struct.*, 2002, **614**, 3–10.

77 N. Berkaïne, P. Reinhardt and M. E. Alikhani, Metal (Ti, Zr, Hf) insertion in the C-H bond of methane: manifestation of an agostic interaction, *Chem. Phys.*, 2008, **343**, 241–249.

78 P. L. A. Popelier, On the full topology of the Laplacian of the electron density, *Coord. Chem. Rev.*, 2000, **197**, 169–189.

79 J. C. S. da Silva, W. B. de Almeida and W. R. Rocha, Theoretical investigation of the structure and nature of the interaction in metal-alkane σ-complexes of the type [M(CO)₅(C₂H₆)] (M = Cr, Mo, and W), *Chem. Phys.*, 2009, **365**, 85–93.

<figure><img src="image_35.png" /></figure>
<figure><img src="image_36.png" /></figure>

80 D. Braga, F. Grepioni, K. Biradha and G. R. Desiraju, Agostic interactions in organic metal compounds. A Cambridge Structural database study, *J. Chem. Soc., Dalton Trans.*, 1996, 3925–3930.

81 J. Pinkas, V. Vargab, I. Císarovac, J. Kubistaa, M. Horaceka and K. Mach, The first thermally stable half-sandwich titanium zwitterionic complex, *J. Organomet. Chem.*, 2007, **692**, 2064–2070.

82 W. E. Billups, S. C. Chang, R. H. Hauge and J. L. Margrave, Detection of a a-Complex in the Reaction of Cobalt Atoms with Methane, *J. Am. Chem. Soc.*, 1993, **115**, 2041–2042.

83 R. H. Crabtree, *The organometallic chemistry of the transition metals*, John Wiley and Sons Inc, 5th edn, 2009, p. 307.

84 E. F. van der Eide, P. Yang and R. M. Bullock, Isolation of two agostic isomers of an organometallic cation: different structures and colors, *Angew. Chem.*, 2013, **125**, 10380–10384.

85 W. W. Lukens, M. R. Smith III and R. A. Andersen, A π-Donor Spectrochemical Series for X in (Me5C5)2TiX, and β-Agostic Interactions in X Et and N(Me)Ph, *J. Am. Chem. Soc.*, 1996, **118**, 1719–1728.

86 See for instance: (a) R. S. Anju, D. K. Roy, B. Mondal, K. Yuvaraj, C. Arivazhagan, K. Saha, B. Varghese and S. Ghosh, Reactivity of Diruthenium and Dirhodium Analogues of Pentaborane – Agostic versus Boratrane Complexes, *Angew. Chem., Int. Ed.*, 2014, **53**, 2873–2877; (b) M. A. Alvarez, D. Garcia-Vivo, M. E. Garcia, M. E. Martinez, A. Ramos and M. A. Ruiz, Reactivity of the r-Agostic Methyl Bridge in the Unsaturated Complex [*Mo₂*(η⁵-C₅H₅)₂(μ-η¹:η²-CH₃)(μ-PCy₂)(CO)₂]: Migratory Behavior and Methylidyne Derivatives*, *Organometallics*, 2008, **27**, 1973–1975.

87 G. Alcaraz and S. Sabo-Etienne, NMR: a good tool to ascertain σ-silane or σ-borane formulations?, *Coord. Chem. Rev.*, 2008, **252**, 2395–2409.

88 See for instance (a) J. A. Calladine, S. B. Duckett, M. W. George, S. L. Matthews, R. N. Perutz, O. Torres and K. Q. Vuong, Manganese Alkane Complexes: An IR and NMR Spectroscopic Investigation, *J. Am. Chem. Soc.*, 2011, **133**, 2303–2310; (b) O. Torres, J. A. Calladine, S. B. Duckett, M. W. George and R. N. Perutz, Detection of s-alkane complexes of manganese by NMR and IR spectroscopy in solution: (η⁵-C₅H₅)Mn(CO)₂(ethane) and (η⁵-C₅H₅)-Mn(CO)₂(isopentane), *Chem. Sci.*, 2015, **6**, 418–424.

89 See for instance H. G. Cho, X. Wang and L. Andrews, The C-H Activation of Methane by Laser-Ablated Zirconium Atoms: CH2ZrH2, the Simplest Carbene Hydride Complex, Agostic Bonding, and (CH₃)₂ZrH₂, *J. Am. Chem. Soc.*, 2005, **127**, 465–473.

90 See for instance F. M. Conroy-Lawis, L. Mole, A. D. Redhouse, S. A. Litster and J. L. Spencer, Synthesis of

Coordinatively Unsaturated Diphosphine Nickel(II) and Palladium(II) β-Agostic Ethyl Cations: X-Ray Crystal Structure of [NiBu<sub>t</sub><sub>2</sub>P(CH<sub>2</sub>)<sub>2</sub>PBu<sub>t</sub><sub>2</sub>(C<sub>2</sub>H<sub>5</sub>)][BF<sub>4</sub>], *J. Am. Chem. Soc.*, 1991, 1601–1603.

91 C. Elschenbroich, *Organometallics: a concise introduction*, Weinheim, New York, Basel, 2nd edn, 1992, p. 168.

92 K. Rachlewicz, S. L. Wang and C. H. Peng, C.H. Hung and L. Latos-Grazynski. Remarkable paramagnetically shifted 1H and 2H NMR spectra of Iron(II) complexes of 2-aza-21-carbaporphyrin: an evidence for agostic interaction, *Inorg. Chem.*, 2003, **42**, 7348–7350.

93 G. A. Luinstra, L. C. Ten Cate, H. J. Heeres, J. W. Pattiasina, A. Meetsmaand and J. H. Teuben, Synthesis and Reactivity of Tervalent Paramagnetic Titanium Compounds (η⁵-C₅Me₅)₂TiR: Molecular Structure of (η⁵-C₅Me₅)₂TiCH₂CMe₃, *Organometallics*, 1991, **10**, 3227–3237.

94 T. S. Thakur and G. R. Desiraju, Misassigned C-HCu agostic interaction in a copper(II) ephedrine derivative is actually a weak, multicentred hydrogen bond, *Chem. Commun.*, 2006, 552–554.

95 S. Noury, X. Krokidis, F. Fuster and B. Silvi, Computational tools for the electron localization function topological analysis, *Comput. Chem.*, 1999, **23**, 597–604.

96 T. A. Keith, *AIMAll (Version 14.10.27)*, TK Gristmill Software, Overland Park KS, USA, 2014 (aim.tkgristmill.com).

97 F. Weigend and R. Ahlrichs, Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy, *Phys. Chem. Chem. Phys.*, 2005, 7, 3297–3307.

98 (a) D. Feller, The role of databases in support of computational chemistry calculations, *J. Comput. Chem.*, 1996, **17**, 1571–1586; (b) K. L. Schuchardt, B. T. Didier, T. Elsethagen, L. Sun, V. Gurumoorthi, J. Chase, J. Li and T. L. Windus, Basis Set Exchange: A Community Database for Computational Sciences, *J. Chem. Inf. Model.*, 2007, **47**, 1045–1052.

99 J.-D. Chai and M. Head-Gordon, Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections, *Phys. Chem. Chem. Phys.*, 2008, **10**, 6615–6620.

100 L. Goerigk and S. Grimme, Efficient and Accurate Double-Hybrid-Meta-GGA Density Functionals—Evaluation with the Extended GMTKN30 Database for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions, *J. Chem. Theory Comput.*, 2011, 7, 291–309.

101 C. Lepetit, J. Poater, M. E. Alikhani, B. Silvi, Y. Canac, J. Contreras-García, M. Solà and R. Chauvin, The Missing Entry in the Agostic-Anagostic Series: Rh(I)-η¹-C Interactions in P(CH)P Pincer Complexes, *Inorg. Chem.*, DOI: 10.1021/acs.inorgchem.5b00069.