Charge Transfer in the Electron Donor−Acceptor Complex BH3NH3

Synergistic Charge Transfer Effect in Ferrous Heme-CO Bonding within Cytochrome P450

We conducted ab initio valence bond (VB) calculations employing the valence bond self-consistent field (VBSCF) and breathing orbital valence bond (BOVB) methods to investigate the nature of the coordination bonding between ferrous heme and carbon monoxide (CO) within cytochrome P450. These calculations revealed the significant influence exerted by both proximal and equatorial ligands on the π-backdonation effect from the heme to the CO. Moreover, our VB calculations unveiled a phenomenon of synergistic charge transfer (sCT). In the case of ferrous heme-CO bonding, the significant stabilization in this sCT arises from cooperative resonance between the VB structures associated with σ donation and π backdonation. Unlike many other ligands, CO possesses the unique ability to establish two mutually perpendicular π-backdonation orbital interaction pairs, leading to an intensified stabilization attributed to σ-π resonance. Furthermore, while of a smaller energy magnitude, sCT due to one π-π pair is also present, contributing to the differential stabilization of ferrous heme-CO bonding.

A Theoretical Perspective on the Photochemistry of Boron-Nitrogen Lewis Adducts

Boron-Nitrogen (B-N) Lewis adducts form a versatile family of compounds with numerous applications in functional molecules. Despite the growing interest in this family of compounds for optoelectronic applications, little is currently known about their photophysics and photochemistry. Even the electronic absorption spectrum of ammonia borane, the textbook example of a B-N Lewis adduct, is unavailable. Given the versatility of the light-induced processes exhibited by these molecules, we propose in this work a detailed theoretical study of the photochemistry and photophysics of simple B-N Lewis adducts. We used advanced techniques in computational photochemistry to identify and characterize the possible photochemical pathways followed by ammonia borane and extended this knowledge to the substituted B-N Lewis adducts pyridine-borane and pyridine-boric acid. The photochemistry observed for this series of molecules allows us to extract qualitative rules to rationalize the light-induced behavior of more complex B-N-containing molecules.

Revisiting the bonding evolution theory: a fresh perspective on the ammonia pyramidal inversion and bond dissociations in ethane and borazane

This work offers a comprehensive and fresh perspective on the bonding evolution theory (BET) framework, originally proposed by Silvi and collaborators [X. Krokidis, S. Noury and B. Silvi, Characterization of elementary chemical processes by catastrophe theory, J. Phys. Chem. A, 1997, 101, 7277-7282]. By underscoring Thom's foundational work, we identify the parametric function characterizing bonding events along a reaction pathway through a three-step sequence to establish such association rigorously, namely: (a) computing the determinant of the Hessian matrix at all potentially degenerate critical points, (b) computing the relative distance between these points, and (c) assigning the unfolding based on these computations and considering the maximum number of critical points for each unfolding. In-depth examination of the ammonia inversion and the dissociation of ethane and ammonia borane molecules yields a striking discovery: no elliptic umbilic flag is detected along the reactive coordinate for any of the systems, contradicting previous reports. Our findings indicate that the core mechanisms of these chemical reactions can be understood using only two folds, the simplest polynomial of Thom's theory, leading to considerable simplification. In contrast to previous reports, no signatures of the elliptic umbilic unfolding were detected in any of the systems examined. This finding dramatically simplifies the topological rationalization of electron rearrangements within the BET framework, opening new approaches for investigating complex reactions.

Bidirectional Charge Transfer at the Heme Iron in Reversible and Quasi-irreversible Cytochrome P450 Inhibition

The coordination bonding between inhibitor ligands and heme iron plays a critical role in disrupting the essential catalytic functions of cytochrome P450 enzymes (P450s). Despite its intrinsic importance and consequential implications for human health, our current understanding of coordination bonding in P450 inhibition remains limited. To address this knowledge gap, we conducted a systematic theoretical analysis of the complexes between a ferric or a ferrous heme model and representative inhibitor ligands. Specifically, we evaluated the charge-transfer (CT) effect within these complexes by employing a series of theoretical methods based on density functional theory (DFT). Through a comprehensive analysis, we unveiled the relative significance of ligand-to-heme forward CT in the ferric and ferrous complexes of reversible inhibitors. In contrast, backward CT dominates over forward CT in the ferrous heme complexes of quasi-irreversible inhibitors. Further analysis using the compact frontier orbital method underscores the elevated electron-accepting abilities of quasi-irreversible inhibitors for π backdonation, which greatly amplifies their binding affinity for the ferrous heme. This study sheds light on the intricate mechanisms underlying P450 inhibition and provides valuable insights for future inhibitor design and development.

Through-Bond-Driven Through-Space Interactions in a Fullerene C60 Noncovalent Dyad: An Unusual Strong Binding between Spherical and Planar π Electron Clouds and Culmination of Dyadic Fractals

Hydrogen-bond-induced π-depletion as a criterion for π-stacking, a configurationally unique noncovalent strategy enabled an unconventional strong binding between the spherical N-fulleropyrrolidine (NFP) and the planar distributions of π electron clouds of three substituted pybates to form noncovalent fulleropyrrolidino-4-(pyrenyl) butanoate dyads of large computed interaction energies, varying between 37.49 and 44.93 kcal/mol. The geometrical distortion/bending of the alkyl tail of pybate in the noncovalent dyad was experimentally corroborated via UV-vis absorption spectroscopy, which translated into spectral broadening along with pronounced shifts in the n-π* transitions of the oxy-substituted pyrene in different solvents, ensuring through-bond interactions. Facile electron transfer through H-bond influenced the dynamic dispersive forces to be active, revealing the supremacy of through-bond over through-space interactions. The analyses of intermolecular forces using an absolutely localized molecular orbital-based energy decomposition analysis (ALMO-EDA) scheme revealed intricate insights into the intermolecular interactions and characteristic charge transfer; the dominance of forward electron transfer (pybate to NFP) over the reverse in offering stabilization was noted. Charge transfer was investigated further from natural bond orbital (NBO) and absolutely localized molecular orbital-based charge-transfer analysis (ALMO-CTA) methods, establishing the supremacy of donor-to-acceptor electron transfer over the reverse (acceptor-to-donor) one. The characteristic self-assembly of the noncovalent dyad in suitable solvents led to the formation of fractal networks via reaction-limited cluster aggregation with a fractal dimension of 2.37. Adoption of constrained molecular dynamics simulations indicated probable wrapping of pybates around NFP, leading to fractal-like assembly.

The stability of covalent dative bond significantly increases with increasing solvent polarity

It is generally expected that a solvent has only marginal effect on the stability of a covalent bond. In this work, we present a combined computational and experimental study showing a surprising stabilization of the covalent/dative bond in Me₃NBH₃ complex with increasing solvent polarity. The results show that for a given complex, its stability correlates with the strength of the bond. Notably, the trends in calculated changes of binding (free) energies, observed with increasing solvent polarity, match the differences in the solvation energies (ΔE^solv) of the complex and isolated fragments. Furthermore, the studies performed on the set of the dative complexes, with different atoms involved in the bond, show a linear correlation between the changes of binding free energies and ΔE^solv. The observed data indicate that the ionic part of the combined ionic-covalent character of the bond is responsible for the stabilizing effects of solvents.

Non covalent complexes are often considerably destabilized in the solvent. Here the authors combine vibrational Raman and NMR spectroscopy with a coupled-cluster computational investigation to show that the solvent polarity enhance the complex stability of a Me₃NBH₃ complex.

Anomeric effect, hyperconjugation and electrostatics: lessons from complexity in a classic stereoelectronic phenomenon

Understanding the interplay of multiple components (steric, electrostatic, stereoelectronic, dispersive, etc.) that define the overall energy, structure, and reactivity of organic molecules can be a daunting task. The task becomes even more difficult when multiple approaches based on different physical premises disagree in their analysis of a multicomponent molecular system. Herein, we will use a classic conformational "oddity", the anomeric effect, to discuss the value of identifying the key contributors to reactivity that can guide chemical predictions. After providing the background related to the relevant types of hyperconjugation and a brief historic outline of the origins of the anomeric effect, we outline variations of its patterns and provide illustrative examples for the role of the anomeric effect in structure, stability, and spectroscopic properties. We show that the complete hyperconjugative model remains superior in explaining the interplay between structure and reactivity. We will use recent controversies regarding the origin of the anomeric effect to start a deeper discussion relevant to any electronic effect. Why are such questions inherently controversial? How to describe a complex quantum system using a model that is "as simple as possible, but no simpler"? What is a fair test for such a model? Perhaps, instead of asking "who is right and who is wrong?" one should ask "why do we disagree?". Stereoelectronic thinking can reconcile quantum complexity with chemical intuition and build the conceptual bridge between structure and reactivity. Even when many factors contribute to the observed structural and conformational trends, electron delocalization is a dominating force when the electronic demand is high (i.e., bonds are breaking as molecules distort from their equilibrium geometries). In these situations, the role of orbital interactions increases to the extent where they can define reactivity. For example, negative hyperconjugation can unleash the "underutilized" stereoelectronic power of unshared electrons (i.e., the lone pairs) to stabilize a developing positive charge at an anomeric carbon. This analysis paves the way for the broader discussion of the omnipresent importance of negative hyperconjugation in oxygen-containing functional groups. From that point of view, the stereoelectronic component of the anomeric effect plays a unique role in guiding reaction design.

σ-Donation and π-Backdonation Effects in Dative Bonds of Main-Group Elements

The nature of donor-acceptor interactions is important for the understanding of dative bonding and can provide vital insights into many chemical processes. Here, we have performed a computational study to elucidate substantial differences between different types of dative interactions. For this purpose, a data set of 20 molecular complexes stabilized by dative bonds was developed (DAT20). A benchmark study that considers many popular density functionals with respect to accurate quantum chemical interaction energies and geometries revealed two different trends between the complexes of DAT20. This behavior was further explored by means of frontier molecular orbitals, extended-transition-state natural orbitals for chemical valence (ETS-NOCV), and natural energy decomposition analysis (NEDA). These methods revealed the extent of the forward and backdonation between the donor and acceptor molecules and how they influence the total interaction energies and molecular geometries. A new classification of dative bonds is suggested.

Anti-Electrostatic Main Group Metal-Metal Bonds That Activate CO2

There has been growing interest in the CO₂ capture and reduction by transition-metal-free catalysts. Here we performed a proof-of-concept study using an ab initio valence bond method called the block-localized wave function (BLW) method. The integrated BLW and density function theory (DFT) computations demonstrated that heterobimetallic Ae⁺/Al(I) (Ae represents alkaline earth metals Mg and Ca) Lewis acid/base combinations without transition metals can facilely capture and activate CO₂. There are two remarkable findings in this study. The first concerns the ionic nature of the metal-metal bonds. The experimentally synthesized low valent aluminum compound with a bidentate β-diketiminate (BDI) ligand, or (BDI)Al(I) in brief, is a Lewis base due to the lone pair on the aluminum cation though overall Al(I) is positively charged. Al(I) can form ionic metal-metal bonds with the alkaline earth metals of the positively charged Lewis acids (BDI)Ae⁺. This type of ionic metal-metal bonds is counterintuitive and antielectrostatic as both metals carry positive charges. The second finding is the CO₂ activation mechanism. (BDI)Al(I) can effectively bind and activate CO₂ by transferring one electron to CO₂, and the resulting complex can be best expressed as [(BDI)Al(I)]⁺[CO₂]^-. The participation of (BDI)Ae⁺ further enhances the capture and activation of CO₂ by (BDI)Al(I).

Nature and Strength of Lewis Acid/Base Interaction in Boron and Nitrogen Trihalides

We have quantum chemically investigated the bonding between archetypical Lewis acids and bases. Our state-of-the-art computations on the X3 B-NY3 Lewis pairs have revealed the origin behind the systematic increase in B-N bond strength as X and Y are varied from F to Cl, Br, I, H. For H3 B-NY3 , the bonding trend is driven by the commonly accepted mechanism of donor-acceptor [HOMO(base)-LUMO(acid)] interaction. Interestingly, for X3 B-NH3 , the bonding mechanism is determined by the energy required to deform the BX3 to the pyramidal geometry it adopts in the adduct. Thus, Lewis acids that can more easily pyramidalize form stronger bonds with Lewis bases. The decrease in the strain energy of pyramidalization on going from BF3 to BI3 is directly caused by the weakening of the B-X bond strength, which stems primarily from the bonding in the plane of the molecule (σ-like) and not in the π system, at variance with the currently accepted mechanism.

Ab initio valence bond theory: A brief history, recent developments, and near future

This Perspective presents a survey of several issues in ab initio valence bond (VB) theory with a primary focus on recent advances made by the Xiamen VB group, including a brief review of the earlier history of the ab initio VB methods, in-depth discussion of algorithms for nonorthogonal orbital optimization in the VB self-consistent field method and VB methods incorporating dynamic electron correlation, along with a concise overview of VB methods for complex systems and VB models for chemical bonding and reactivity, and an outlook of opportunities and challenges for the near future of the VB theory.

Energy Decomposition Analysis for Metal Surface-Adsorbate Interactions by Block Localized Wave Functions

The energy decomposition analysis based on block localized wave functions (BLW-EDA) allows one to gain physical insight into the nature of chemical bonding, decomposing the interaction energy in (1) a "frozen" term, accounting for the attraction due to electrostatic and dispersion interactions, modulated by Pauli repulsion, (2) the variationally assessed polarization energy, and (3) the charge transfer. This method has so far been applied to gas- and condensed-phase molecular systems. However, its standard version is not compatible with fractionally occupied orbitals (i.e., electronic smearing) and, as a consequence, cannot be applied to metallic surfaces. In this work, we propose a simple and practical extension of BLW-EDA to fractionally occupied orbitals, termed Ensemble BLW-EDA. As illustrative examples, we have applied the developed method to analyze the nature of the interaction of various adsorbates on Pt(111), ranging from physisorbed water to strongly chemisorbed ethylene. Our results show that polarization and charge transfer both contribute significantly at the adsorption minimum for all studied systems. The energy decomposition analysis provides details with respect to competing adsorption sites (e.g., CO on atop vs hollow sites) and elucidates the respective importance of polarization and charge transfer for the increased adsorption energy of H₂S compared to H₂O. Our development will enable a deeper understanding of the impact of charge transfer on catalytic processes in general.

Pleading for a Dual Molecular-Orbital/Valence-Bond Culture

Electron pairs through the looking glass might well discover that they can show two faces, one delocalized or the other localized, and that both are perfectly correct. Going back and forth between these two representations, according to which one is the most relevant and insightful for the case at hand, is easy and essential to get a complete understanding of electronic structure.

Diabatic Population Matrix Formalism for Performing Molecular Mechanics Style Simulations with Multiple Electronic States

An accurate description of nonbonded interactions is important in investigating dynamics of molecular systems. In many situations, fixed point charge models are successfully applied to explaining various chemical phenomena. However, these models with conventional formulations will not be appropriate in elucidating the detailed dynamics during nonadiabatic events. This is mainly because the chemical properties of any molecule, especially its electronic populations, significantly change with respect to molecular distortions in the vicinity of the surface crossing. To overcome this issue in molecular simulations yet within the framework of the fixed point charge model, we define a diabatic electronic population matrix and substitute it for the conventional adiabatic partial charges. We show that this matrix can be readily utilized toward attaining more reliable descriptions of Coulombic interactions, in combination with the interpolation formalism for obtaining the intramolecular interaction potential. We demonstrate how the mixed formalism with the diabatic charges and the interpolation can be applied to molecular simulations by conducting adiabatic and nonadiabatic molecular dynamics trajectory calculations of the green fluorescent protein chromophore anion in aqueous environment.

Charge-Shift Bonding Emerges as a Distinct Electron-Pair Bonding Family from Both Valence Bond and Molecular Orbital Theories

The charge-shift bonding (CSB) concept was originally discovered through valence bond (VB) calculations. Later, CSB was found to have signatures in atoms-in-molecules and electron-localization-function and in experimental electron density measurements. However, the CSB concept has never been derived from a molecular orbital (MO)-based theory. We now provide a proof of principle that an MO-based approach enables one to derive the CSB family alongside the distinctly different classical family of covalent bonds. In this bridging energy decomposition analysis, the covalent-ionic resonance energy, RECS, of a bond is extracted by cloning an MO-based purely covalent reference state, which is a constrained two-configuration wave function. The energy gap between this reference state and the variational TCSCF ground state yields numerical values for RECS, which correlate with the values obtained at the VBSCF level. This simple MO-based method, which only takes care of static electron correlation, is already sufficient for distinguishing the classical covalent or polar-covalent bonds from charge-shift bonds. The equivalence of the VB and MO-based methods is further demonstrated when both methods are augmented by dynamic correlation. Thus, it is shown from both MO and VB perspectives that the bonding in the CSB family does not arise from electron correlation. Considering that the existence of the CSB family is associated also with quite a few experimental observations that we already reviewed ( Shaik , S. , Danovich , D. , Wu , W. , and Hiberty , P. C. Nat. Chem. , 2009 , 1 , 443 - 449 ), the new bonding concept has passed by now two stringent tests. This derivation, on the one hand, supports the new concept and on the other, it creates bridges between the two main theories of electronic structure.

Energy decomposition analysis approaches and their evaluation on prototypical protein–drug interaction patterns

The partitioning of the interaction energy into chemical components such as electrostatics, polarization, and charge transfer is possible with energy decomposition analysis approaches. We review and evaluate these for biomolecular applications.

Large hyperconjugation in strained systems

Carbocations can appear as transient species, for instance, in elimination reactions and various rearrangements. Hyperconjugation (or conjugation) can then stabilize the cationic character and form a partial π bond. The effect of the electronic delocalization from strained substituents to a carbocation part was calculated. Very large hyperconjugation was found, sometimes more than 80 kcal mol(-1) , which is much larger than typical conjugation effects (56 kcal mol(-1) for the allyl cation).

On the origin of internal rotation in ammonia borane

The internal rotation in ammonia borane (AB) was studied on the basis of natural orbitals for chemical valence (NOCV) and eigenvectors for Pauli repulsion (NOPR). We found that the total hyperconjugation stabilization (ca. 5 kcal mol⁻¹), based on the charge transfer from the occupied σ (B–H) orbitals into the empty σ*(N–H), slightly favors the staggered conformation over the eclipsed one; however, the barrier to internal rotation in ammonia borane can be understood predominantly in a ‘classical’ way, as originating from the steric (Pauli) repulsion contributions (of the kinetic origin) that act solely between N–H and B–H bonds. Repulsion between the lone pair of ammonia and the adjacent B–H bonds was found to be dominant in absolute terms; however, it does not determine the rotational barrier. Similar conclusions on the role of CH↔HC repulsion appeared to be valid for isoelectronic ethane.

Figure

Design of Lewis acid-base complex: enhancing the stability and first hyperpolarizability of large excess electron compound

In the present paper, a new type of Lewis acid-base complex BX3...Li@Calix[4]pyrrole (X = H and F) was designed and assembled based on electride molecule Li@calix[4]pyrrole (as a Lewis base) and the electron deficient molecule BX3 (as a Lewis acid) by employing quantum mechanical calculation. The new Lewis acid-base complex offers an interesting push-excess electron-pull (P-e-P) framework to enhance the stability and nonlinear optical (NLO) response. To measure the nonlinear optical response, static first hyperpolarizabilities (β 0) are exhibited. Significantly, point-face assembled Lewis acid-base complex BF3...Li@Calix[4]pyrrole (II) has considerable first hyperpolarizabilities (β 0) value (1.4  ×  106 a.u.), which is about 117 times larger than reported 11,721 a.u. of electride Li@Calix[4]pyrrole. Further investigations show that, in BX3...Li@Calix[4]pyrrole with P-e-P framework, a strong charge-transfer transition from the ground state to the excited state contributes to the enhancement of first hyperpolarizability. Theory calculation of enthalpies of reaction (ΔrH0) at 298 K demonstrates that it is feasible to synthetize the complexes BX3...Li@Calix[4]pyrrole. In addition, compared with Li@Calix[4]pyrrole, the vertical ionization potential (VIP) and HOMO-LUMO gap of BX3...Li@Calix[4]pyrrole have obviously increased, due to the introduction of the Lewis acid molecule BX3. The novel Lewis acid-base NLO complex possesses not only a large nonlinear optical response but also higher stability.

Hyperconjugation in Carbocations, a BLW Study with DFT approximation

The geometry of ethyl cation is discussed, and the hyperconjugation effect in carbocations is evaluated at the B3LYP/6-311G(d) level. The Block Localized Wavefunction (BLW) method is used for all evaluations of the hyperconjugation, considered as the energy gained by the delocalization onto the C⁺ atom. This energy is defined as the energy difference between the delocalized (standard) calculation, where the electrons are freely delocalized, and a localized form where the positive charge sits on the carbon center. It is evaluated for 18 carbocations, including conjugated systems. In these cases we were particularly interested in the additional stabilization brought by hyperconjugative effects. Among other effects, the β-silicon effect is computed. Hyperconjugation amounts in several cases to an energy similar to conjugation effects.

Why Cyclooctatetraene Is Highly Stabilized: The Importance of "Two-Way" (Double) Hyperconjugation

Despite its highly nonplanar geometry, the tub-shaped D2d cyclooctatetraene (COT) minimum is far from being an unconjugated polyene model devoid of important π interactions. The warped skeleton of D2d COT results in the large stabilization (41.1 kcal/mol) revealed by its isodesmic bond separation energy (BSE). This originates largely from the "two-way" hyperconjugation, back and forth across the C-C single bonds, between the CC/CH σ(σ*) and the C═C (π*)π orbitals. These hyperconjugative effects compensate for the substantial, but not complete, loss of π conjugation upon ring puckering. C-C single bond rotation of 1,3-butadiene involves a similar interplay between π conjugation and "two-way" double hyperconjugation and serves as a simple model for the inversion of D2d to D4h COT. The perpendicular rotational transition states of many other systems, e.g., the allyl cation, styrene, biphenyl, and ethene, are stabilized similarly by "two-way" hyperconjugation.

A VALENCE BOND APPROACH BASED ON LEWIS STRUCTURES

In this paper, a valence bond method based on Lewis structures, called LVB, is presented. The method uses a Slater determinant expansion of doubly occupied orbitals for describing a Lewis structure, where two orbital sets, semi-localized orbitals, called bond orbitals, and localized hybrid atomic orbitals (HAOs), are employed. The levels of LVB method are fashioned as LVBS, LVBSD, etc. LVBS involves only the single bond orbital replacements with HAOs, while LVBSD involves also double replacements, and so on. Tests of three examples, methane, methylene, and benzene, show that the LVB method at both of LVBS and LVBSD levels gives results that match those of the VBSCF method very well, even though the form of LVB wave function is much compact.

Dispersion-corrected energy decomposition analysis for intermolecular interactions based on the BLW and dDXDM methods

As the simplest variant of the valence bond (VB) theory, the block-localized wave function (BLW) method defines the intermediate electron-localized state self-consistently at the DFT level and can be used to explore the nature of intermolecular interactions in terms of several physically intuitive energy components. Yet, it is unclear how the dispersion interaction affects such a kind of energy decomposition analysis (EDA) as standard density functional approximations neglect the long-range dispersion attractive interactions. Three electron densities corresponding to the initial electron-localized state, optimal electron-localized state, and final electron-delocalized state are involved in the BLW-ED approach; a density-dependent dispersion correction, such as the recently proposed dDXDM approach, can thus uniquely probe the impact of the long-range dispersion effect on EDA results computed at the DFT level. In this paper, we incorporate the dDXDM dispersion corrections into the BLW-ED approach and investigate a range of representative systems such as hydrogen-bonding systems, acid-base pairs, and van der Waals complexes. Results show that both the polarization and charge-transfer energies are little affected by the inclusion of the long-range dispersion effect, which thus can be regarded as an independent energy component in EDA.