Intrinsic bond orbitals (IBO) are localized molecular orbitals giving exact and non-empirical representations of wave functions. They are obtained by unitary transformation and form an orthogonal set of orbitals localized on a minimal number of atoms. IBOs present an intuitive and unbiased interpretation of chemical bonding with naturally arising Lewis structures.^[1] For this reason IBOs have been successfully employed for the elucidation of molecular structures^[1] and electron flow along the intrinsic reaction coordinate (IRC).^[2] IBOs have also found application as Wannier functions in the study of solids.^[3]

Theory

The IBO method entails molecular wave-functions calculated using self-consistent field (SCF) methods such as Kohn-Sham density functional theory (DFT) which are expressed as linear combinations of localized molecular orbitals.^[1][4]

In order to arrive at IBOs, intrinsic atomic orbitals (IAOs) are first calculated as representations of a molecular wave function for which each IAO can be assigned to a specific atom. This allows for a chemically intuitive orbital picture as opposed to the commonly used large and diffuse basis sets for the construction of more complex molecular wavefunctions.^[4]

IAOs are constructed from tabulated free-atom AOs of standard basis-sets under consideration of the molecular environment. This yields polarized atomic orbitals that resemble the free-atom AOs as much as possible, before orthonormalization of the polarized AOs results in the set of IAOs. IAOs are thus a minimal basis for a given molecule in which atomic contributions can be distinctly assigned. The sum of all IAOs spans exactly over the molecular orbitals which renders them an exact representation of the wavefunction.^[4] Since IAOs are associated with a specific atom, they can provide atom specific properties such as the partial charge. Compared to other charges, such as the Mulliken charge, the IAO charges are independent of the employed basis set.^[1]

IBOs are constructed as a linear combination over IAOs with the condition of minimizing the number of atoms over which the orbital charge is spread. Each IBO can thereby be divided into the contributions of the atoms as the electronic occupation ${\displaystyle n_{A}(i')}$ of orbital ${\displaystyle i'}$ on atom ${\displaystyle A}$. The localization is performed in the spirit of the Pipek-Mezey localization scheme, maximizing a localization functional ${\displaystyle L}$.^[1]

${\displaystyle L=\sum _{i}^{occ}\sum _{A}^{atom}[n_{A}(i')]^{p}}$

with ${\displaystyle p=4}$ or ${\displaystyle 2}$. While the choice of the exponent ${\displaystyle p}$ does not affect the resulting IBOs in most cases, the choice of ${\displaystyle p=4}$ localizes the orbitals in aromatic systems unlike ${\displaystyle p=2}$.^[1]

The process of IBO construction is performed by unitary tranfomation of canonical MOs, which ensures that the IBOs remain an exact and physically accurate representation of the molecular wavefunction due to the invariance of Slater determinant wavefunctions towards unitary rotations.

${\displaystyle |i'\rangle =\sum _{i}^{occ}|i\rangle U_{ii'}}$

The unitary matrix ${\displaystyle U_{ii'}}$, which produces the localized IBOs upon matrix multiplication with set of occupied MOs ${\displaystyle |i\rangle }$, is thereby chosen to effectively minimize spread of IBOs over the atoms of a molecule. The product is a set of localized IBOs, closely resembling the chemically intuitive shapes of molecular orbitals, allowing for distinction of bond types, atomic contributions and polarization.^[1]

Application in structure and bonding

In his original paper introducing IBOs, Knizia showed the versatility of his method for describing not only classical bonding situations, such as the σ and π bond, but also aromatic systems and non-trivial bonds. The differentiation of σ and π bonds in acrylic acid is possible based on IBO geometries, as are the identification of the IBOs corresponding to the oxygen lone pairs. Benzene provided an example of a delocalized aromatic system to test the IBO method. Apart from the C-C and C-H σ-bonds, the six electron π-system is expressed as three delocalized IBOs. Representation of non-Lewis bonding was demonstrated on diborane B₂H₆, with one IBO stretching over B-H-B, corresponding to the 3-center-2-electron bond.^[1]

IBO of the 3-center-2-electron-bond in B₂H₆ visualized using IBOView. Recreated from reference ^[1].

IBO of acrylic acid visualized using IBOview. Recreated from reference ^[1].

IBOs of benzene visualized using IBOview. Recreated from reference ^[1].

Transition metal compounds

IBO analysis was used to explain the stability of electron rich gold-carbene complexes, mimicking reactive intermediates in gold catalysis. While these complexes are sometimes depicted with a Au-C double bond, representing the sigma donation of the carbene and π backbonding of Au, IBO analysis points towards a minimal amount of π-backbonding with the respective orbital mainly localized on Au. The σ-donating carbene orbital is likewise strongly polarized towards C. Stabilization of the compound thus occurs through strong donation of the aromatic carbene substituents into the carbene carbon p-orbitals, outcompeting the Au-π-backdonation. IBO analysis was thus able to negate the double bond character of the gold-carbene complexes and provided deep insight into the electronic structure of Cy₃P-Au-C(4-OMe-C₆H₄)₂) (Cy = cyclohexyl).^[4]

IBOs of Cy₃P-Au-C(4-OMe-C₆H₄)₂), visualized using IBOview. Recreated from reference.^[4]

The π-backbonding character was again evaluated for gold-vinylidene complexes, as another common type of gold catalysis intermediates. IBO analysis revealed significantly stronger π-backbonding for the gold-vinylidenes compared to the gold-carbenes. This was attributed to the geometric inability of aromatic vinylidene substituents to compete with Au for π-interactions since the respective orbitals are perpendicular to each other.^[5]

Knizia and Klein similarly employed IBO for the analysis of [Fe(CO)₃NO]^–. The even polarization of IBOs between Fe and N points towards a covalently bonded NO ligand. The double bond occurs via two d-p π-interactions and results in a formal Fe⁰ center. Confirmed by further calculations, IBO proved as a fast and straightforward method to interpret bonding in this case.^[6] Making use of the low computational cost, a Cloke-Wilson rearrangement catalyzed by [Fe(CO)₃NO]^– was investigated by constructing the IBOs for every stationary point along the IRC. It was found that one of the Fe-NO π bonds takes active part in catalysis by electron transfer to and from the substrate, explaining the unique catalytic activity of [Fe(CO)₃NO]^– compared to the isoelectronic [Fe(CO)₄]^2–.^[7]

Apart from the above mentioned compounds, the IBO method has been employed to investigate various other transition metal complexes, such as gold-diarylallenylidenes^[4] or diplatinum diboranyl complexes,^[8] proving as a valuable tool to gain insight into the extent and nature of bonding.

Main group compounds

Structures of phosphaaluminirene^[9] (1) and distonic radical ion^[10] (2).

IBO analysis has been employed in main group chemistry to elucidate oftentimes non-trivial electronic structure. The bonding of phosphaaluminirenes was, for example, investigated showing a 3-center-2π-electron bond of the AlCP cycle.^[11] Further application was found for confirming the distonic nature of a phosphorus containing radical cation reported by Chen et al. (see figure). While the IAO charge analysis yielded a positive charge on the chelated P, IBOs showed the localization of the unpaired electron on the other P atom, confirming the spatial separation of radical site and charge.^[12]

Another example is the elucidation of the electronic structure of the hexamethylbenzene dication. Three π bonding IBOs were found between the basal C₅Me₅ plane and the apical C, reminiscent of Cp* coordination complexes. The three π bonds are thereby polarized towards the apical C, which in turn coordinates to a CH₃⁺ cation with its lone pair. IBO analysis therefore revealed the Lewis-acidic and Lewis-basic character of the apical C.^[13]

IBOs of the hexamethyl dication visualized using IBOview. Recreated from reference.^[13]

Applications of IBO for cluster compounds have included zirconium doped boron clusters. IBO analysis showed, that the unusual stability of the neutral ZrB₁₂ cluster stems from several multicenter σ bonds. The B-B σ bonding orbitals extend to the central Zr atom, forming the mulicenter bonds. This example displays the method's aptitude to analyze cluster compounds and multicenter bonding.^[14]

Electron flow along the IRC

IBOs along the reaction coordinate of the SN₂-type self-exchange reaction of H₃CCl and Cl^–. Recreated from reference ^[2].

In 2015, Knizia and Klein introduced the analysis of electron flow in reactions with IBO as a non-empirical and straight-forward method of evaluating curly arrow mechanisms. Since IBOs are exact representations of Kohn-Sham wavefunctions, they can provide physical conformation for curly arrows mechanisms based on first-principles. IBOs usually represent chemical bonds and lone pairs, this method allows for elucidation of bond rearrangements in terms of the elemental steps and their sequence. By calculating the root mean square deviations of the partial charge distributions compared to the initial charge distribution, IBOs taking active part in a reaction can be distinguished from those that remain unchanged along the IRC.^[2]

Knizia and Klein demonstrate the versatility of this method in their original report, first presenting a simple SN₂-type self-exchange reaction of H₃CCl and Cl^–, followed by the migration of π bonds in a substitution reaction (SN₂) and π- to σ-bond transformations in a Claisen rearrangement. Electron flow can be easily followed by observing the migration of an IBO and bond types are easily distinguished based on the geometries of the IBOs. The value of IBO analysis along the IRC^[2] especially shows for complex reactions, such as a cyclopropanation reaction with only one transition state and without intermediates, reported by Haven et al.^[16] Calculations by Knizia and Klein yielded a precise curly arrow mechanism for this reaction.^[2]

Closed-shell systems

Examples of IBO analysis along the IRC included the investigation of C-H bond activation by gold-vinylidene complexes. Through this method, it is possible to discern between concerted and stepwise reactions. The previously thought single step C-H activation reaction was in this case revealed to consist of three distinct phases: i) hydride transfer, ii) C-C bond formation and iii) sigma to pi rearrangement of the lone pair coordinated to Au.^[17]

Comparison of four curly arrow textbook mechanism (left) and seven curly arrow IBO derived mechanism (right) of an epoxidation with a peracid.^[18]

C-H activation by a gold-vinylidene complex. IBO enables distinction of three distinct phases: i) hydride transfer, ii) C-C bond formation and iii) sigma to pi rearrangement of the lone pair coordinated to Au.^[17]

Other reports of IBO analysis along the IRC include the elucidation and confirmation of a previously proposed mechanism for a [3,3]-sigmatropic rearrangement of a Au(I)-vinyl species^[19] or the epoxidation of alkene by peracids.^[18] For the latter, the textbook four curly arrow mechanism was found to be physically inaccurate. Instead, seven changing IBOs were found, yielding an ideal mechanism featuring seven curly arrows.^[18] The combination of IBO analysis with other computational methods, such as natural bond orbital (NBO) analysis for a Ti-catalyzed pyrrole synthesis^[20] or natural localized molecular orbital (NLMO) analysis for an intramolecular cycloaddition of a phosphaalkene to an arene has likewise led to insightful results regarding the specifics of the reaction mechanisms.^[21]  

Schematic representation of hydrogen atom transfer (HAT) versus proton coupled electron transfer (PCET).^[22]

See also

• Localized molecular orbitals

References