Valence Bond Study Research Articles

Dichotomy of Delocalization/Localization and Charge-Shift Bonding in Germanazene and its Heavier Group 14 Analogues: a Valence Bond Study.

We present here a valence bond analysis of structure and π-delocalization in Ge3 (NH)3 , which models germanazene that was prepared by Power et al. To get a broader perspective, we explore the entire E3 (NH)3 series (E=C, Si, Ge, Sn, Pb). Thus, while (4n+2)π systems of carbon rings are aromatic with cyclic π-delocalization, the E3 (NH)3 rings are dominated by a nonbonded structure, wherein π-lone pairs are localized on the N atoms. Nevertheless, these molecules enjoy large covalent-ionic resonance energies of 153.0, 86.6, 74.2, 61.2, and 58.9 kcal/mol, respectively, for E=C, Si, Ge, Sn, Pb. The covalent-ionic mixing in E3 (NH)3 creates π-systems, which are stabilized by charge-shift bonding. Thus, unlike in benzene, in Ge3 (NH)3 delocalization of π-electron pairs of the N atoms is primarily confined to the domains of their adjacent Ge atoms. These features carry over to the substituted germanazene, Ge3 (NAr)3 (Ar=Ph).

The effect of immediate environment on bond strength of different bond types-A valence bond study.

The ability to design catalysis largely depends on our understanding of the electrostatic effect of the surrounding on the bonds participating in the reaction. Here, we used a simplistic model of point charges (PCs) to determine a set of rules guiding how to construct PC-bond arrangement that can strengthen or weaken different chemical bonds. Using valence bond theory to calculate the in situ bond energies, we show that the effect of the PC mainly depends on the bond's dipole moment irrespective of its type (being covalent or charge shift). That is, polar bonds are getting stronger or weaker depending on the signand location of the PC, whereas non- or weakly polar bonds become stronger or weaker depending only on the location of the PC and to a smaller extent compared with polar bonds. We also show that for polar bonds, the maximal bond strengthening and weakening effect can be achieved when the PC is placed along the bond axis, as close as possible to the more and less polarizable atom/fragment, respectively. Finally, due to the stabilizing effects of polarizability, we show that, overall, it is easier to cause bond strengthening compared with bond weakening. Particularly, for polar bonds, bond strengthening is larger than bond weakening obtained by an oppositely signed PC. These rules should be useful in the future design of catalysis in, e.g., enzyme active sites.

Hydride Abstraction as the Rate-Limiting Step of the Irreversible Inhibition of Monoamine Oxidase B by Rasagiline and Selegiline: A Computational Empirical Valence Bond Study.

Monoamine oxidases (MAOs) catalyze the degradation of a very broad range of biogenic and dietary amines including many neurotransmitters in the brain, whose imbalance is extensively linked with the biochemical pathology of various neurological disorders, and are, accordingly, used as primary pharmacological targets to treat these debilitating cognitive diseases. Still, despite this practical significance, the precise molecular mechanism underlying the irreversible MAO inhibition with clinically used propargylamine inhibitors rasagiline and selegiline is still not unambiguously determined, which hinders the rational design of improved inhibitors devoid of side effects current drugs are experiencing. To address this challenge, we present empirical valence bond QM/MM simulations of the rate-limiting step of the MAO inhibition involving the hydride anion transfer from the inhibitor α-carbon onto the N5 atom of the flavin adenin dinucleotide (FAD) cofactor. The proposed mechanism is strongly supported by the obtained free energy profiles, which confirm a higher reactivity of selegiline over rasagiline, while the calculated difference in the activation Gibbs energies of ΔΔG‡ = 3.1 kcal mol−1 is found to be in very good agreement with that from the measured literature kinact values that predict a 1.7 kcal mol−1 higher selegiline reactivity. Given the similarity with the hydride transfer mechanism during the MAO catalytic activity, these results verify that both rasagiline and selegiline are mechanism-based irreversible inhibitors and offer guidelines in designing new and improved inhibitors, which are all clinically employed in treating a variety of neuropsychiatric and neurodegenerative conditions.

Captodative Substitution Enhances the Diradical Character of Compounds, Reduces Aromaticity, and Controls Single-Molecule Conductivity Patterns: A Valence Bond Study.

The present contribution uses a valence bond (VB) perspective to consider the captodative substitution strategy, a method to enhance the diradical character of (potentially aromatic) compounds. We confirm the qualitative reasoning that has generally been used to rationalize the diradical-character-enhancing effect of captodative substitution: this type of substitution scheme disproportionally stabilizes specific Dewar/diradical(oid) VB structures, thus increasing their weight in the full ground-state wave function. Furthermore, we assess the effect of captodative substitution on the aromaticity of the considered compound. We observe a clear trade-off between diradical character and aromaticity for our model systems: as one of these properties increases, the other decreases. This finding is especially significant within the field of single-molecule electronics because it enables unification of the previously observed inverse proportionality between the aromaticity of a compound and the magnitude of conductance through that molecule, with the observed proportionality between diradical character and the magnitude of conductance associated with a compound. To some extent, both properties, i.e., aromaticity and diradical character, appear to be the flip-sides of the same coin.

A valence bond study of the activation of methyl halides bonds by electric fields

In the present work, the activation of methyl halides bonds under experience of an external electric field (EEF) is explained from the Valence Bond theory perspective. The dissociation mechanism of C–X bonds (X [Formula: see text] Cl, Br, I) influenced by a homogeneous and a heterogeneous field placed parallel to the bond axis is presented. For all examples, an increase in the electric field strength have similar consequences: (i) the decrease of the energy depth along the dissociation path, (ii) an increase of the equilibrium interatomic distance (at high EEFs), and (iii) the transition from a homolytic to a heterolytic dissociation after some field magnitude. These general behaviors are explained through the curve crossing between the ionic and the covalent structure at some field strength.

Capturing the Role of Explicit Solvent in the Dimerization of RuV(bda) Water Oxidation Catalysts

AbstractA ground‐breaking empirical valence bond study for a soluble transition‐metal complex is presented. The full reaction of catalyst monomers approaching and reacting in the RuV oxidation state were studied. Analysis of the solvation shell in the reactant and along the reaction coordinate revealed that the oxo itself is hydrophobic, which adds a significant driving force to form the dimer. The effect of the solvent on the reaction between the prereactive dimer and the product was small. The solvent seems to lower the barrier for the isoquinoline (isoq) complex while it is increased for pyridines. By comparing the reaction in the gas phase and solution, the proposed π‐stacking interaction of the isoq ligands is found to be entirely driven by the water medium.

Capturing the Role of Explicit Solvent in the Dimerization of RuV (bda) Water Oxidation Catalysts.

A ground-breaking empirical valence bond study for a soluble transition-metal complex is presented. The full reaction of catalyst monomers approaching and reacting in the RuV oxidation state were studied. Analysis of the solvation shell in the reactant and along the reaction coordinate revealed that the oxo itself is hydrophobic, which adds a significant driving force to form the dimer. The effect of the solvent on the reaction between the prereactive dimer and the product was small. The solvent seems to lower the barrier for the isoquinoline (isoq) complex while it is increased for pyridines. By comparing the reaction in the gas phase and solution, the proposed π-stacking interaction of the isoq ligands is found to be entirely driven by the water medium.

Bonding in B2 and B2+: Insights from full configuration interaction and valence bond studies

Full Configuration Interaction (Full-CI) and Valence Bond Self-Consistent Field (VBSCF) methods have been used to study the electronic structure and bonding in B2 and B2+ molecules. The bonding analysis based on these calculations shows that the B2 molecule is stabilised due to the formation of a double σ bond, one strong σ-bond together with one second weaker σ-bond, and two weak π bonds. Upon ionization one π electron is removed from the system and B2+ is formed, which has a one electron σ bond, instead of a π bond. It has been shown that a few carefully chosen VB configurations are enough to describe the bonding; with these structures, geometrical parameters as well as dissociation energies of these unusual molecular species are in agreement with full-CI results.

The driving force for Π-bond localization and bond alternation in trisannelated benzenes.

To investigate the factors that may cause bond length alternation and π-bond localization in annelated benzenes, ab initio valence bond calculations were performed. The results reveal that the bond length alternation of annelated benzene is determined by the strain-induced hybridization change from the partially optimized geometries, in which the central benzene ring is constrained to a regular hexagon, to the equilibrium geometries rather than the previously recognized re-hybridization effect that the carbon atoms in the central ring are deviated from sp2 hybridization. Meanwhile, the π-π interaction also provides a sort of driving force, which facilitates bond length alternation, which in turn magnifies π-bond localization. A subsequent potential energy curve study shows that σ-strain and π-π interactions have different mechanisms for the effect.

Multistate empirical valence bond study of temperature and confinement effects on proton transfer in water inside hydrophobic nanochannels.

Microscopic characteristics of an aqueous excess proton in a wide range of thermodynamic states, from low density amorphous ices (down to 100 K) to high temperature liquids under the critical point (up to 600 K), placed inside hydrophobic graphene slabs at the nanometric scale (with interplate distances between 3.1 and 0.7 nm wide) have been analyzed by means of molecular dynamics simulations. Water-proton and carbon-proton forces were modeled with a multistate empirical valence bond method. Densities between 0.07 and 0.02 Å(-3) have been considered. As a general trend, we observed a competition between effects of confinement and temperature on structure and dynamical properties of the lone proton. Confinement has strong influence on the local structure of the proton, whereas the main effect of temperature on proton properties is observed on its dynamics, with significant variation of proton transfer rates, proton diffusion coefficients, and characteristic frequencies of vibrational motions. Proton transfer is an activated process with energy barriers between 1 and 10 kJ/mol for both proton transfer and diffusion, depending of the temperature range considered and also on the interplate distance. Arrhenius-like behavior of the transfer rates and of proton diffusion are clearly observed for states above 100 K. Spectral densities of proton species indicated that in all states Zundel-like and Eigen-like complexes survive at some extent. © 2016 Wiley Periodicals, Inc.

Can Aromaticity Coexist with Diradical Character? An Ab Initio Valence Bond Study of S2N2 and Related 6π‐Electron Four‐Membered Rings E2N2 and E42+ (E=S, Se, Te)

A series of 6π-electron 4-center species, E(2)N(2) and E(4)(2+) (E=S, Se, Te) is studied by means of ab initio valence bond methods with the aims of settling some controversies on 1) the diradical character of these molecules and 2) the radical sites, E or N, of the preferred diradical structure. It was found that for all molecules, the cumulated weights of the two possible diradical structures are always important and close to 50 %, making these molecules comparable to ozone in terms of diradical character. While the two diradical structures are degenerate in the E(4)(2+) dications, they have on the contrary strongly unequal weights in the E(2)N(2) neutral molecules. In these three molecules, the electronic structure is dominated by one diradical structure, in which the radical sites are the two nitrogen atoms, while the other diradical structure is much less important. The ordering of the various VB structures in terms of their calculated weights is confirmed by the relative energies of individual VB structures. In all cases, the major diradical structure (or both diradical structures when they are degenerate) is (are) the lowest one(s), while the covalent VB structures lie higher in energy. The vertical resonance energies are considerable in S(2)N(2) and S(4)(2+), about 80 % of the estimated value for benzene, and diminish as one goes down the periodic table (S→Se→Te). This confirms the aromatic character of these species, as already demonstrated for S(2)N(2) on the basis of magnetic criteria. This and the high weights and stabilities of one or both diradical structures in all systems indicates that aromaticity and diradical character do not exclude each other, contrary to what is usually claimed. Furthermore, it is shown that the diradical structures find their place in a collective electron flow responsible for the ring currents in the π system of these species.

The Nature of the Idealized Triple Bonds Between Principal Elements and the σ Origins of Trans-Bent Geometries-A Valence Bond Study.

We describe herein a valence bond (VB) study of 27 triply bonded molecules of the general type X≡Y, where X and Y are main element atoms/fragments from groups 13-15 in the periodic table. The following conclusions were derived from the computational data: (a) Single π-bond and double π-bond energies for the entire set correlate with the "molecular electronegativity", which is the sum of the X and Y electronegativites for X≡Y. The correlation with the molecular electronegativity establishes a simple rule of periodicity: π-bonding strength generally increases from left to right in a period and decreases down a column in the periodic table. (b) The σ frame invariably prefers trans bending, while π-bonding gets destabilized and opposes the trans distortion. In HC≡CH, the π-bonding destabilization overrides the propensity of the σ frame to distort, while in the higher row molecules, the σ frame wins out and establishes trans-bent molecules with 2(1)/2 bonds, in accord with recent experimental evidence based on solid state (29)Si NMR of the Sekiguchi compound. Thus, in the trans-bent molecules "less bonds pay more". (c) All of the π bonds show significant bonding contributions from the resonance energy due to covalent-ionic mixing. This quantity is shown to correlate linearly with the corresponding "molecular electronegativity" and to reflect the mechanism required to satisfy the equilibrium condition for the bond. The π bonds for molecules possessing high molecular electronegativity are charge-shift bonds, wherein bonding is dominated by the resonance energy of the covalent and ionic forms, rather than by either form by itself.

Trends in Aromatic Oxidation Reactions Catalyzed by Cytochrome P450 Enzymes: A Valence Bond Modeling.

The mixed density functional theory (DFT) and valence bond study described herein focuses on the activation of 17 benzene derivatives by the active species of Cytochrome P450, so-called Compound I (Cpd I), as well as by the methoxy radical, as a potentially simple model of Cpd I (Jones, J. P.; Mysinger, M.; Korzekwa, K. R. Drug Metab. Dispos. 2002, 30, 7-12). Valence bond modeling is employed to rationalize the P450 mechanism and its spin-state selectivity from first principles of electronic structure and to predict activation energies independently, using easily accessible properties of the reactants: the singlet-triplet excitation energies, the ionization potentials of the aromatics, and the electron affinity of Cpd I and/or of the methoxy radical. It is shown that the valence bond model rationalizes all the mechanistic aspects and predicts activation barriers (for 35 reactions) with reasonable accuracy compared to the DFT barriers with an average deviation of ±1.0 kcal·mol(-1) (for DFT barriers, see: Bathelt, C. M.; Ridder, L.; Mulholland, A. J.; Harvey, J. N. Org. Biomol. Chem. 2004, 2, 2998-3005). The valence bond modeling also reveals the mechanistic similarities between the P450 Cpd I and methoxy reactions and enables one to make predictions of barriers for reactions from other studies.

Bonding Conundrums in the C2 Molecule: A Valence Bond Study.

The ab initio VB study for the electronic structure of the C2 molecule in the ground state is presented in this work. VB calculations involving 78 chemically relevant VB structures can predict the bonding energy of C2 quite well. Sequentially, a VBCIS calculation provides spectroscopic parameters that are very close to full CI calculated values in the same basis set. Furthermore, the analysis of the bonding scheme shows that a triply bonded structure is the major one in terms of weights, and the lowest in energy at the equilibrium distance. The second structure in terms of weights is an ethylene-like structure, displaying a σ + π double bond. The structure with two suspended π bonds but no σ bond contributes only marginally to the ground state. This ordering of weights for the VB structures describing the C2 molecule is shown to be consistent with the shape of the molecular orbitals and with the multireference character of the ground state. With the triply bonded bonding scheme, the natures of the π and σ bonds are investigated, and then the corresponding "in situ" bond strengths are estimated. The contribution of the covalent-ionic resonance energy to π and σ bonding is revealed and discussed.

The Inverted Bond in [1.1.1]Propellane is a Charge‐Shift Bond

Bonded or not bonded? An ab initio valence bond study of [1.1.1]propellane shows that the two bridgehead carbons are linked by a strong and direct sigma bond that is neither classically covalent nor classically ionic, but rather a charge-shift bond, in which the covalent-ionic resonance energy plays the major role. As such, the central bond of [1.1.1]propellane closely resembles the single bond of difluorine.

Covalent Excited States of Polyenes C2nH2n+2 (n = 2−8) and Polyenyl Radicals C2n-1H2n+1 (n = 2−8): An Ab Initio Valence Bond Study

Covalent Excited States of Polyenes C2nH2n+2 (n = 2−8) and Polyenyl Radicals C2n-1H2n+1 (n = 2−8): An Ab Initio Valence Bond Study

A Valence Bond Study of the Low‐Lying States of the NF Molecule

The electronic structures of the three lowest-lying states of NF are investigated by means of modern valence bond (VB) methods such as the VB self-consistent field (VBSCF), breathing orbital VB (BOVB), and VB configuration interaction (VBCI) methods. The wave functions for the three states are expressed in terms of 9-12 VB structures, which can be further condensed into three or four classical Lewis structures, whose weights are quantitatively estimated. Despite the compactness of the wave functions, the BOVB and VBCI methods reproduce the spectroscopic properties and dipole moments of the three states well, in good agreement with previous computational studies and experimental values. By analogy to the isoelectronic O(2) molecule, the ground state (3)Sigma(-) possesses both a sigma bond and 3-electron pi bonds. However, here the polar sigma bond contributes the most to the overall bonding. It is augmented by a fractional (19%) contribution of three-electron pi bonding that arises from pi charge transfer from fluorine to nitrogen. In the singlet (1)Delta and (1)Sigma(+) excited states the pi-bonding component is classically covalent, and it contributes 28% and 37% to the overall bonding picture for the two states, respectively. The resonance energies are calculated and reveal that pi bonding contributes at least 24, 35 and 42 kcal mol(-1) to the total bonding energies of the (3)Sigma(-), (1)Delta and (1)Sigma(+) states, respectively. Some unusual properties of the NF molecule, like the equilibrium distance shortening and bonding energy increasing upon excitation, the counterintuitive values of the dipole moments and the reversal of the dipole moments as the bond is stretched, are interpreted in the light of the simple valence bond picture. The overall polarity of the molecule is very small in the ground state, and is opposite to the relative electronegativity of N vs F in the singlet excited states. The values of the dipole moments in the three states are quantitatively accounted for by the calculated weights of the VB structures.

Testing the Validity of the Conventional Resonance Model for Protonated Carbonyl, Imine and Thiocarbonyl Compounds. An Ab Initio Valence Bond Study

The conventional resonance model describes protonated carbonyls, imines, and thiocarbonyls by a superposition of two structures, one pi polar-covalent and the other of carbenium type. The validity of this model is clearly supported by high level valence bond calculations, giving a 32% weight for the carbenium form in protonated carbonyl, 19% in protonated formamine and thioformaldehyde. The carbenium form is further stabilized by pi-donating substituents. Solvation effects do not fundamentally change the gas-phase picture.

Electronic polarization effects on charge carriers in anthracene: A valence bond study

A semiempirical quantum-chemical model based on a fragment orbital formalism is presented to assess molecular parameters relevant to charge transport in organic crystals. The mixed valence bond/Hartree--Fock approach provides an efficient integrated framework to evaluate the electronic polarization effects induced by localized charge carriers and the associated impact on the matrix elements mediating electron migration in the hopping regime. This formalism, applied here to anthracene clusters of increasing sizes and dimensionalities, yields the electrostatic and polarization contributions to the total interaction energy of the neutral and charged aggregates and leads to a reduction in the effective bandwidth by $\ensuremath{\sim}10%--20%$ as a result of the polarization cloud.

The Menshutkin Reaction in the Gas Phase and in Aqueous Solution: A Valence Bond Study

The recently developed (L. Song, W. Wu, Q. Zhang, S. Shaik, J. Phys. Chem. A 2004, 108, 6017-6024) valence bond method coupled to a polarized continuum model (VBPCM) is applied to the Menshutkin reaction, NH3+CH3Cl-->CH3NH3(+)+Cl-, in the gas phase and in aqueous solution. The computed barriers and reaction energies at the level of the breathing orbital VB method (P. C. Hiberty, J. P. Flament, E. Noizet, Chem. Phys. Lett. 1992, 189, 259), BOVB and VBPCM//BOVB, are comparable to CCSD(T) and CCSD(T)//PCM results and to experimental values in solution. The gas-phase reaction is endothermic and leads to an ion-pair complex via a late transition state. By contrast, the reaction in the aqueous phase is exothermic and leads to separate solvated ions as reaction products, via an early transition state. The VB calculations provide also the reactivity parameters needed to apply the valence bond state correlation diagram method, VBSCD (S. Shaik, A. Shurki, Angew. Chem. Int. Ed. 1999, 38, 586). It is shown that the reactivity parameters along with their semiempirical derivations provide together a satisfactory qualitative and quantitative account of the barriers.