Bond-Centered Orbital Basis

• Bond-centered orbital bases are representations that localize electron density on bonds, dimers, or multi-center units instead of individual atoms.
• Various methods such as IAO/IBO construction, Wannier function localization, and wavelet-based approaches achieve chemical clarity by focusing on bond regions.
• These representations improve charge analysis and computational efficiency in both molecular systems and correlated materials like VO₂.

A bond-centered orbital basis is an orbital representation in which chemically relevant degrees of freedom are organized around bonds, bond regions, dimers, trimers, or other pairwise units rather than solely around isolated atoms. In the literature surveyed here, that idea appears in several distinct but related forms: localized occupied orbitals derived from atom-centered minimal bases, bond-centered Wannier functions used as correlated subspaces, pairwise-antisymmetrized many-electron bases assembled from diatomic building blocks, and correlation-based clusterings of localized orbitals that define bonds as strongly correlated orbital sets. The common objective is to replace opaque delocalized expansions with a basis in which bonding, lone pairs, bond polarity, and multi-center structures emerge directly from the wave function or from a low-energy effective Hamiltonian (Knizia, 2013).

1. Bond-centered representations and their scope

In one widely used formulation, the bond-centered basis is obtained in two steps. First, one constructs a chemically interpretable atom-centered basis; second, one localizes the occupied space with respect to that basis so that each occupied orbital is concentrated on as few atoms as possible. In Knizia’s framework, intrinsic atomic orbitals (IAOs) are polarized core and valence orbitals, one set per atom, that exactly span the occupied molecular-orbital space of the self-consistent-field solution. Intrinsic bond orbitals (IBOs) are then obtained by unitary rotation of the occupied space so as to maximize locality in the IAO populations, yielding orbitals that are typically localized on one or two atoms, or on three atoms in 3-center bonds (Knizia, 2013).

Other constructions generalize the same idea. In system-adapted bases built from product plane waves and wavelet localization, the final wavelet-localized orbitals are cell-centered, but linear combinations of neighboring cell functions provide bond-like orbitals and capture both occupied and correlation space (Baker et al., 2017). In first-principles hybrid atomic-orbital bases, hybrid atomic Wannier orbitals remain atom-centered but are locked to the immediate atomic neighborhood, so that each nearest-neighbor bond is described predominantly by one orbital pair and one tight-binding parameter (Hossain et al., 2020). In bond-centered dynamical mean-field theory for VO$_2$, the correlated orbitals are explicitly centered on V–V bond midpoints and are used as the local impurity degrees of freedom (Mlkvik et al., 2024). At the many-electron level, the finite-basis Spectral Theory of Chemical Bonding represents polyatomic Hamiltonians by pairwise-antisymmetrized bond blocks built from orthogonalized products of atomic eigenstates, so the basis is bond-centered in a pair-Hamiltonian rather than a one-electron-orbital sense (Mills, 2022).

This variety indicates that “bond-centered orbital basis” is not a single algorithm. It is a family of representations in which chemically meaningful pair or multi-center units are made explicit at the level of orbitals, Wannier functions, or many-electron pair states.

2. IAOs and IBOs as the canonical bond-centered construction

The most direct formal route from a conventional SCF calculation to a bond-centered basis is the IAO/IBO construction. Starting from occupied molecular orbitals $\lvert i\rangle$ in a large basis $B_1$, Knizia introduces a minimal free-atom basis $B_2$ containing only core and valence shells. After projection into and back out of $B_2$, one obtains depolarized occupied orbitals $\lvert \tilde{i}\rangle$ and the associated projector $\tilde{O}$. Each free-atom orbital $\lvert \tilde{\rho}\rangle$ is then mapped into the full basis through

$$\lvert \rho \rangle = \Bigl( O \tilde O + (1-O)(1-\tilde O) \Bigr) P_{12}\lvert \tilde \rho\rangle,$$

followed by symmetric orthogonalization. The resulting IAOs are atom-centered, minimal, polarized by the molecular environment, and exactly span the occupied space together with the implied minimal virtual complement (Knizia, 2013).

The bond-centered step is the localization of occupied orbitals in the IAO basis. For a localized occupied orbital $\lvert i'\rangle$, the population on atom $\lvert i\rangle$0 is

$\lvert i\rangle$1

and the localization functional is

$\lvert i\rangle$2

Maximizing $\lvert i\rangle$3 over unitary rotations of the occupied space yields the IBOs. The fourth power suppresses fractional distributions and leads to discrete localization, especially important in aromatic systems (Knizia, 2013).

This construction is atom-centered at the IAO level and bond-centered at the IBO level. Its chemical content is unusually explicit. Simple bonds can be expressed to $\lvert i\rangle$4 with IAOs on only two atoms; acrylic acid gives 19 occupied MOs mapped to 19 IBOs, of which 16 of 19 are localized on one or two centers with $\lvert i\rangle$5 of the charge; benzene yields completely localized CC and CH $\lvert i\rangle$6 bonds but a delocalized $\lvert i\rangle$7 system requiring at least 4 centers; diborane produces six normal BH $\lvert i\rangle$8 bonds and two BHB 2e-3c bridge bonds; SO$\lvert i\rangle$9 gives one $B_1$0 and one $B_1$1 bond per S–O, with the $B_1$2 bonds having $B_1$3 oxygen character and $B_1$4 sulfur (Knizia, 2013).

The same basis supports charge analysis. In the IAO basis, the partial charge is

$B_1$5

and these charges are reported as basis-set stable and chemically sensible, with a linear correlation to C 1s shifts of $B_1$6, or $B_1$7 without two outliers (Knizia, 2013). This makes the IAO basis simultaneously a population-analysis basis and the reference frame for constructing bond-centered occupied orbitals.

3. Bond-directed localized bases beyond IBOs

Several later developments retain the same strategic pattern—first define a local basis, then identify the bond-carrying subspace—but depart from the occupied-only localization problem.

A system-adapted basis based on product plane waves first approximates natural orbitals cheaply from HF or DFT orbitals multiplied by low-momentum plane waves, then “chops” the delocalized functions into cell components with wavelets, and finally compresses them by principal component analysis into wavelet-localized orbitals (WLOs). These WLOs are cell-centered, but the bonding pattern is encoded in their intercell combinations. In one-dimensional benchmark systems, the method requires only 2–3 basis functions per electron to achieve chemical accuracy, and a localized WLO basis preserves this compactness: for H$B_1$8 near equilibrium, PCA with $B_1$9 reduces the basis to 14 WLOs with energy error $B_2$0 kcal/mol at $B_2$1; for H$B_2$2, $B_2$3–$B_2$4 WLOs give energy error $B_2$5 kcal/mol (Baker et al., 2017). This suggests a bond-centered basis can be viewed not only as a localized occupied space but also as a localized correlated one-particle space.

A closely related but more explicitly chemical development is the first-principles hybrid atomic-orbital basis built from finite first-moment matrices. Hybrid atomic orbitals are obtained as approximate common eigenstates of the $B_2$6, $B_2$7, and $B_2$8 first-moment matrices, and their hybridization and orientation can be a-priori tuned as per their anticipated neighbourhood. Their Wannier counterparts, the hybrid atomic Wannier orbitals, form an orthonormal multi-orbital tight-binding basis resembling hybrid atomic-orbitals locked to their immediate atomic neighborhood, while spanning the sub-space of Kohn–Sham states. In that basis, nearest-neighbour bonds involve no more than two orbitals irrespective of their orientation, and the spatial extent of self-energy correction to tight-binding parameters is localized mostly within the third nearest neighbourhood (Hossain et al., 2020).

The “maximally valent orbitals” program makes the bond-centered criterion fully explicit. In a basis of custom hybridized atomic Wannier orbitals, the authors evaluate Mayer’s bond order and optimize hybrid orientations so that the variance of orbital-pair bond-order contributions along each nearest-neighbour coordination segment is maximized. The resulting maximally valent hybrid atomic orbitals identify the exact orientation of the major overlapping orbitals such that they maximally represent the covalent interactions throughout the system. In systems with non-ideal bond angles, those optimal orbitals can deviate from the geometrical coordination segments, giving a first-principles realization of bent bonds; the paper also emphasizes that different physical aspects of covalent interactions are not necessarily represented by a single unique set of atomic or bonding orbitals (De et al., 2022).

4. Correlation-based definitions of bond-centered orbital clusters

A different line of work defines a bond-centered basis through correlation structure rather than through localization functionals. In the multiorbital correlation theory of the chemical bond, a proposed bond partition $B_2$9 is judged by the entropy-based correlation

$B_2$0

defined for an arbitrary partition $B_2$1 of a localized orbital set $B_2$2. A bond is then a cluster of orbitals with strong internal multiorbital correlation and weak correlation to other clusters. Successive bipartitioning of the orbital set reveals whether there is a well-defined partition $B_2$3 at which the increase in $B_2$4 changes from small to large; if such a partition exists, it furnishes a rigorous bond clustering. If it does not exist, then the bonding picture is intrinsically ambiguous (Szalay et al., 2016).

This framework supplies a formal interpretation of aromatic and antiaromatic systems. In benzene, pyrrole, furan, and thiophene, the ring $B_2$5 orbitals form a single multiorbital cluster, whereas in borole and cyclobutadiene the corresponding ring orbitals split into smaller units. For C$B_2$6, the non-existence of well-defined multiorbital correlation clustering is explicitly identified as a reason for the debated bonding picture (Szalay et al., 2016).

Orbital-entanglement analysis offers a related diagnostic. One-orbital entropy $B_2$7 and mutual information $B_2$8 quantify how strongly orbitals participate in bond formation and bond breaking. For diatomics, strongly correlated bonding/antibonding pairs appear as pairs with large $B_2$9 and rising $\lvert \tilde{i}\rangle$0 along dissociation; the authors also note that for larger reactive systems the possible necessity for localizing molecular orbitals is currently under investigation (Boguslawski et al., 2013). This suggests that a bond-centered orbital basis can be validated, or rejected, by the way it concentrates entanglement into chemically meaningful units.

An orthogonal-valence-bond reformulation of CASSCF wave functions reaches a similar conclusion from a different angle. Canonical CASSCF MOs are orthogonally localized into orthogonal atomic orbitals, producing OVB configuration state functions whose squared coefficients have a direct probabilistic interpretation. In ethane, ethene, and N$\lvert \tilde{i}\rangle$1, this yields clear single-, double-, and triple-bond patterns. In C$\lvert \tilde{i}\rangle$2, however, no single neutral or ionic configuration dominates near equilibrium, and the electron structure does not allow one to say that C$\lvert \tilde{i}\rangle$3 in its ground state has a double, or triple, or even a quadruple bond (Sax, 2023).

5. Bond-centered correlated subspaces in materials

In correlated materials, a bond-centered orbital basis can be used as the correlated subspace itself. In VO$\lvert \tilde{i}\rangle$4, a bond-centered Wannier basis is constructed by unitary transformation of V-centered $\lvert \tilde{i}\rangle$5 Wannier functions so that the new orbitals are centered on the midpoints of nearest-neighbor V–V bonds along the rutile $\lvert \tilde{i}\rangle$6 axis. For each bond, there are three such orbitals: one predominantly $\lvert \tilde{i}\rangle$7-like and two predominantly $\lvert \tilde{i}\rangle$8-like. Single-site DMFT is then performed with the bond center as the impurity “site”, so that a local Hubbard interaction on each bond-centered orbital effectively includes both on-site and inter-site interactions in the original atomic basis (Mlkvik et al., 2024).

This representation allows rutile and M1 VO$\lvert \tilde{i}\rangle$9 to be treated on exactly the same footing, and later work extends the same framework to M2 and T phases and to the full structural phase space of rutile-based VO$\tilde{O}$0 with moderate computational cost and without pre-patterning the structure into dimerized and undimerized V–V pairs. In the M2 phase, the dimerized chains host a singlet-insulator and the zigzag-distorted chains host a Mott-insulator, but the two site types are strongly coupled and the metal-insulator transition occurs concomitantly for both (Mlkvik et al., 27 Mar 2026).

Related materials cases show that bond-centered descriptions can also be cluster-centered. In metallic Li$\tilde{O}$1VS$\tilde{O}$2, the low-temperature ribbon-chain phase is interpreted in terms of multiple three-centered two-electron $\tilde{O}$3 bonds on V trimers, so the natural low-energy basis is built from trimer bonding, nonbonding, and antibonding combinations rather than from individual V sites (Katayama et al., 2018). In GaTa$\tilde{O}$4Se$\tilde{O}$5, the relevant low-energy objects are $\tilde{O}$6 molecular states on Ta$\tilde{O}$7 clusters, but the ordered phase is interpreted as a staggered intercluster dimerization pattern and a spin-orbital valence bond ground state, which is again most naturally expressed in a bond- or dimer-centered language (Yang et al., 2022).

A more formal many-electron analogue appears in the finite-basis pair formulation of the Spectral Theory of Chemical Bonding. There the polyatomic Hamiltonian is assembled from archived diatomic matrices in an orthonormal antisymmetrized product basis obtained by symmetric orthogonalization of atomic-product valence-bond states. The resulting pair blocks function as transferable bond-centered building elements and reproduce the valence and geometry of simple hydrocarbons without introducing three- or four-center two-electron integrals (Mills, 2022).

6. Ambiguities, limitations, and interpretive consequences

The surveyed literature does not support a unique universal bond-centered orbital basis. Instead, it shows that different constructions privilege different physical quantities. IBOs privilege locality in occupied-space atomic populations (Knizia, 2013). Wavelet-localized system-adapted bases privilege compactness of the correlated one-particle space (Baker et al., 2017). Maximally valent orbitals privilege concentration of Mayer/Wiberg bond order into a minimal number of orbital pairs (De et al., 2022). Bond-centered DMFT bases privilege inclusion of intra-dimer or intra-bond correlations in a computationally cheap local impurity problem (Mlkvik et al., 2024). Spatially partitioned QMC bases privilege numerical efficiency by using atomic-centered functions near nuclei and B-splines in the interstitial, i.e. bond, regions (Luo et al., 2018).

This plurality has direct chemical consequences. In systems with ordinary localized $\tilde{O}$8 and $\tilde{O}$9 bonds, these viewpoints often coincide closely. In benzene, the $\lvert \tilde{\rho}\rangle$0 framework localizes cleanly whereas the $\lvert \tilde{\rho}\rangle$1 system remains delocalized (Knizia, 2013). In cyclopropane, perfectly localized two-center bonds coexist with unusual p-rich bond character (Knizia, 2013). In water and ammonia, maximally localized and maximally valent bond orbitals need not have the same orientation (De et al., 2022). In C$\lvert \tilde{\rho}\rangle$2, the absence of a stable multiorbital correlation clustering and the spread of weight over many orthogonal-valence-bond configurations mean that no clean integer bond-centered decomposition is available (Szalay et al., 2016, Sax, 2023).

A plausible implication is that a bond-centered orbital basis is best understood as an analysis-dependent compression of the many-electron problem rather than as an invariant observable. What is stable across the literature is not the detailed orbital shape but the organizing principle: chemically meaningful bonding units are identified by concentrating occupancy, bond order, or correlation into a small set of localized degrees of freedom, and those degrees of freedom may be atom-centered, bond-centered, or cluster-centered depending on the phenomenon under study.