Diboron Bonds Between BX3 (X=H,F, CH3) and BYZ2 (Y=H,F; Z=CO,N2, CNH)

The ability of B atoms on two different molecules to engage with one another in a noncovalent diboron bond is studied by ab initio calculations. Due to electron donation from its substituents, the trivalent B atom of BYZ₂ (Z=CO, N₂, and CNH; Y=H and F) has the ability to in turn donate charge to the B of a BX₃ molecule (X=H, F, and CH₃), thus forming a B⋅⋅⋅B diboron bond. These bonds are of two different strengths and character. BH(CO)₂ and BH(CNH)₂, and their fluorosubstituted analogues BF(CO)₂ and BF(CNH)₂, engage in a typical noncovalent bond with B(CH₃)₃ and BF₃, with interaction energies in the 3–8 kcal/mol range. Certain other combinations result in a much stronger diboron bond, in the 26–44 kcal/mol range, and with a high degree of covalent character. Bonds of this type occur when BH₃ is added to BH(CO)₂, BH(CNH)₂, BH(N₂)₂, and BF(CO)₂, or in the complexes of BH(N₂)₂ with B(CH₃)₃ and BF₃. The weaker noncovalent bonds are held together by roughly equal electrostatic and dispersion components, complemented by smaller polarization energy, while polarization is primarily responsible for the stronger ones.     

Bonding analysis of metal-metal multiple bonds in R3M-M'R3 (M, M' = Cr, Mo, W; R = Cl, NMe2)

The bonding situation of homonuclear and heteronuclear metal-metal multiple bonds in R(3)M-M'R(3) (M, M' = Cr, Mo, W; R = Cl, NMe(2)) is investigated by density functional theory (DFT) calculations, with the help of energy decomposition analysis (EDA). The M-M' bond strength increases as M and M' become heavier. The strongest bond is predicted for the 5d-5d tungsten complexes (NMe(2))(3)W-W(NMe(2))(3) (D(e) = 103.6 kcal/mol) and Cl(3)W-WCl(3) (D(e) = 99.8 kcal/mol). Although the heteronuclear molecules with polar M-M' bonds are not known experimentally, the predicted bond dissociation energies of up to 94.1 kcal/mol for (NMe(2))(3)Mo-W(NMe(2))(3) indicate that they are stable enough to be isolated in the condensed phase. The results of the EDA show that the stronger R(3)M-M'R(3) bonds for heavier metal atoms can be ascribed to the larger electrostatic interaction caused by effective attraction between the expanding valence orbitals in one metal atom and the more positively charged nucleus in the other metal atom. The orbital interaction reveal that the covalency of the homonuclear and heteronuclear R(3)M-M'R(3) bonds is due to genuine triple bonds with one σ- and one degenerate π-symmetric component. The metal-metal bonds may be classified as triple bonds where π-bonding is much stronger than σ-bonding; however, the largest attraction comes from the quasiclassical contribution to the metal-metal bonding. The heterodimetallic species show only moderate polarity and their properties and stabilities are intermediate between the corresponding homodimetallic species, a fact which should allow for the experimental isolation of heterodinuclear species. CASPT2 calculations of Cl(3)M-MCl(3) (M = Cr, Mo, W) support the assignment of the molecules as triply bonded systems.     

Direct estimate of the strength of conjugation and hyperconjugation by the energy decomposition analysis method

The intrinsic strength of pi interactions in conjugated and hyperconjugated molecules has been calculated using density functional theory by energy decomposition analysis (EDA) of the interaction energy between the conjugating fragments. The results of the EDA of the trans-polyenes H2C=CH-(HC=CH)n-CH=CH2 (n = 1-3) show that the strength of pi conjugation for each C=C moiety is higher than in trans-1,3-butadiene. The absolute values for the conjugation between Si=Si pi bonds are around two-thirds of the conjugation between C=C bonds but the relative contributions of DeltaE pi to DeltaE orb in the all-silicon systems are higher than in the carbon compounds. The pi conjugation between C=C and C=O or C=NH bonds in H2C=CH--C(H)=O and H2C=CH-C(H)=NH is comparable to the strength of the conjugation between C=C bonds. The pi conjugation in H2C=CH-C(R)=O decreases when R = Me, OH, and NH2 while it increases when R = halogen. The hyperconjugation in ethane is around a quarter as strong as the pi conjugation in ethyne. Very strong hyperconjugation is found in the central C-C bonds in cubylcubane and tetrahedranyltetrahedrane. The hyperconjugation in substituted ethanes X3C-CY3 (X,Y = Me, SiH3, F, Cl) is stronger than in the parent compound particularly when X,Y = SiH3 and Cl. The hyperconjugation in donor-acceptor-substituted ethanes may be very strong; the largest DeltaE pi value was calculated for (SiH3)3C-CCl3 in which the hyperconjugation is stronger than the conjugation in ethene. The breakdown of the hyperconjugation in X3C-CY3 shows that donation of the donor-substituted moiety to the acceptor group is as expected the most important contribution but the reverse interaction is not negligible. The relative strengths of the pi interactions between two C=C double bonds, one C=C double bond and CH3 or CMe3 substituents, and between two CH3 or CMe3 groups, which are separated by one C-C single bond, are in a ratio of 4:2:1. Very strong hyperconjugation is found in HC[triple bond]C-C(SiH3)3 and HC[triple bond]C-CCl3. The extra stabilization of alkenes and alkynes with central multiple bonds over their terminal isomers coming from hyperconjugation is bigger than the total energy difference between the isomeric species. The hyperconjugation in Me-C(R)=O is half as strong as the conjugation in H2C=CH-C(R)=O and shows the same trend for different substituents R. Bond energies and lengths should not be used as indicators of the strength of hyperconjugation because the effect of sigma interactions and electrostatic forces may compensate for the hyperconjugative effect.     

Kinetics and thermodynamics of H. transfer from (eta5-C5R5)Cr(CO)3H (R = Ph,Me, H) to methyl methacrylate and styrene

The rates of H/D exchange have been measured between (a) the activated olefins methyl methacrylate-d(5) and styrene-d(8), and (b) the Cr hydrides (eta(5)-C(5)Ph(5))Cr(CO)(3)H (2a), (eta(5)-C(5)Me(5))Cr(CO)(3)H (2b), and (eta(5)-C(5)H(5))Cr(CO)(3)H (2c). With a large excess of the deuterated olefin the first exchange goes to completion before subsequent exchanges begin, at a rate first order in olefin and in hydride. (Hydrogenation is insignificant except with styrene and CpCr(CO)(3)H; in most cases, the radicals arising from the first H. transfer are too hindered to abstract another H. .) Statistical corrections give the rate constants k(reinit) for H. transfer to the olefin from the hydride. With MMA, k(reinit) decreases substantially as the steric bulk of the hydride increases; with styrene, the steric bulk of the hydride has little effect. At longer times, the reaction of MMA or styrene with 2a gives the corresponding metalloradical 1a as termination depletes the concentration of the methyl isobutyryl radical 3 or the alpha-methylbenzyl radical 4; computer simulation of [1a] as f(t) gives an estimate of k(tr), the rate constant for H. transfer from 3 or 4 back to Cr. These rate constants imply a DeltaG (50 degrees C) of +11 kcal/mol for H. transfer from 2a to MMA, and a DeltaG (50 degrees C) of +10 kcal/mol for H. transfer from 2a to styrene. The CH(3)CN pK(a) of 2a, 11.7, implies a BDE for its Cr-H bond of 59.6 kcal/mol, and DFT calculations give 58.2 kcal/mol for the Cr-H bond in 2c. In combination the kinetic DeltaG values, the experimental BDE for 2a, and the calculated DeltaS values for H. transfer imply a C-H BDE of 45.6 kcal/mol for the methyl isobutyryl radical 3 (close to the DFT-calculated 49.5 kcal/mol), and a C-H BDE of 47.9 kcal/mol for the alpha-methylbenzyl radical 4 (close to the DFT-calculated 49.9 kcal/mol). A solvent cage model suggests 46.1 kcal/mol as the C-H BDE for the chain-carrying radical in MMA polymerization.     

Heats of formation of boron hydride anions and dianions and their ammonium salts [BnHmy-][NH4+]y with y=1-2

Thermochemical parameters of the closo boron hydride BnHn2- dianions, with n=5-12, the B3H8- and B11H14- anions, and the B5H9 and B10H14 neutral species were predicted by high-level ab initio electronic structure calculations. Total atomization energies obtained from coupled-cluster CCSD(T)/complete basis set (CBS) extrapolated energies, plus additional corrections were used to predict the heats of formation of the simplest BnHmy- species in the gas phase in kcal/mol at 298 K: DeltaHf(B3H8-)=-23.1+/-1.0; DeltaHf(B5H52-)=119.4+/-1.5; DeltaHf(B6H62-)=64.1+/-1.5; and DeltaHf(B5H9)=24.1+/-1.5. The heats of formation of the larger species were evaluated by the G3 method from hydrogenation reactions (values at 298 K, in kcal/mol with estimated error bars of+/-3 kcal/mol): DeltaHf(B7H72-)=51.8; DeltaHf(B8H82-)=46.1; DeltaHf(B9H92-)=24.4; DeltaHf(B10H102-)=-12.5; DeltaHf(B11H112-)=-11.8; DeltaHf(B12H122-)=-86.3; DeltaHf(B11H14-)=-57.3; and DeltaHf(B10H14)=18.7. A linear correlation between atomization energies of the dianions and energies of the BH units was found. The heats of formation of the ammonium salts of the anions and dianions were predicted using lattice energies (UL) calculated from an empirical expression based on ionic volumes. The UL values (0 K) of the BnHn2- dianions range from 319 to 372 kcal/mol. The values of UL for the B3H8- and B11H14- anions are 113 and 135 kcal/mol, respectively. The calculated lattice energies and gas-phase heats of formation of the constituent ions were used to predict the heats of formation of the ammonium crystal salts [BnHmy-][NH4+]y. These results were used to evaluate the thermodynamics of the H2 release reactions from the ammonium hydro-borate salts.     

Why does perfluorination render bicyclo[2.2.0]hex-1(4)-ene stable toward dimerization? Calculations provide the answers

B3LYP calculations with two different basis sets have been performed to understand why bicyclo[2.2.0]hex-1(4)-ene (1a) undergoes dimerization with DeltaH(++) = 11.5 kcal/mol, but dimerization of perfluorobicyclo[2.2.0]hex-1(4)-ene (1b) has never been observed. The former reaction is computed to be exothermic by 37.2 kcal/mol, whereas the latter is calculated to be endothermic by 7.4 kcal/mol. The 44.6 kcal/mol difference between the enthalpies of these two reactions can be dissected into contributions of 24.5 kcal/mol for the difference between the enthalpies for forming diradical intermediates 2a and 2b and 20.1 kcal/mol for cyclization of 2a and 2b to, respectively, 3a and 3b. The latter enthalpy difference is largely attributable to repulsions between the endo-fluorines in the dimer, although the exo-fluorines also are found to contribute. The former enthalpy difference is attributable to the difference between the dissociation enthalpies of the pi bonds in 1a and 1b, which is shown to amount to 16 +/- 1 kcal/mol. About 25% of the stronger pi bond in fluoroalkene 1b is found to be due to hyperconjugation of the eight C-F bonds in 1b with the filled pi orbital. However, the major contributor to the stronger pi bond in 1b is shown to be the unfavorable interaction that results when a pyramidalized radical center is syn to a C-F bond. Both of these effects, which contribute to the greater strength of the pi bond in 1b, relative to that in 1a, are analyzed and discussed.     

Complexes of HArF and AuX (X = F,Cl, Br,I). Comparison of H-bonds,halogen bonds,F-shared bonds and covalent bonds

The various sorts of complexes in which HArF and AuX (X = F, Cl, Br, I) can engage are probed by MP2/aug-cc-pVTZ calculations. The most weakly bound are those containing a halogen bond (XB) of the AuX⋯FArH sort, with binding energies less than 8 kcal/mol. H-bonded dimers FArH⋯XAu are a little stronger, held together by some 12 kcal/mol. Being the most strongly bound places the F atom of HArF roughly midway between Ar and Au in an F-shaped structure, bound by some 43–54 kcal/mol. The last sort of product involves atomic rearrangements wherein the H atom migrates from Ar to Au, followed by formation of a covalent Ar–Au bond. The resulting molecular unit is stabilized by 30–40 kcal/mol relative to the original HArF and AuX reactants. The H-bonded dimers are held together by an unusually large polarization component, surpassing electrostatic attraction, while dispersion predominates for the halogen bonds. Perturbations of the geometries and stretching frequencies offer a ready means of distinguishing the different types of complexes by spectroscopic techniques.     

Electronic structure of germanium monohydrides Ge(n)H, n = 1-3

Quantum chemical calculations were applied to investigate the electronic structure of germanium hydrides, Ge(n)H (n = 1, 2, 3), their cations, and anions. Computations using a multiconfigurational quasi-degenerate perturbation approach (MCQDPT2) based on complete active space wave functions (CASSCF), multireference perturbation theory (MRMP2), and density functional theory reveal that Ge(2)H has a (2)B(1) ground state with a doublet-quartet gap of approximately 39 kcal/mol. A quasidegenerate (2)A(1) state has been derived to be 2 kcal/mol above the ground state (MCQDPT2/aug-cc-pVTZ). In the case of the cation Ge(3)H(+) and anion Ge(3)H(-), singlet low-lying electronic states are derived, that is, (1)A' and (1)A(1), respectively. The singlet-triplet energy gap is estimated to 6 kcal/mol for the cation. An "Atoms in Molecules" (AIM) analysis shows a certain positive charge on the Ge(n) (n = 1, 2, 3) unit in its hydrides, in accordance with the NBO analysis. The topologies of the electron density of the germanium hydrides are different from that of the lithium-doped counterparts. On the basis of our electron localization function (ELF) analysis, the Ge-H bond in Ge(2)H is characterized as a three-center-two-electron bond. Some key thermochemical parameters of Ge(n)H have also been derived.     

Orbital overlap and chemical bonding

The chemical bonds in the diatomic molecules Li(2)-F(2) and Na(2)-Cl(2) at different bond lengths have been analyzed by the energy decomposition analysis (EDA) method using DFT calculations at the BP86/TZ2P level. The interatomic interactions are discussed in terms of quasiclassical electrostatic interactions DeltaE(elstat), Pauli repulsion DeltaE(Pauli) and attractive orbital interactions DeltaE(orb). The energy terms are compared with the orbital overlaps at different interatomic distances. The quasiclassical electrostatic interactions between two electrons occupying 1s, 2s, 2p(sigma), and 2p(pi) orbitals have been calculated and the results are analyzed and discussed. It is shown that the equilibrium distances of the covalent bonds are not determined by the maximum overlap of the sigma valence orbitals, which nearly always has its largest value at clearly shorter distances than the equilibrium bond length. The crucial interaction that prevents shorter bonds is not the loss of attractive interactions, but a sharp increase in the Pauli repulsion between electrons in valence orbitals. The attractive interactions of DeltaE(orb) and the repulsive interactions of DeltaE(Pauli) are both determined by the orbital overlap. The net effect of the two terms depends on the occupation of the valence orbitals, but the onset of attractive orbital interactions occurs at longer distances than Pauli repulsion, because overlap of occupied orbitals with vacant orbitals starts earlier than overlap between occupied orbitals. The contribution of DeltaE(elstat) in most nonpolar covalent bonds is strongly attractive. This comes from the deviation of quasiclassical electron-electron repulsion and nuclear-electron attraction from Coulomb's law for point charges. The actual strength of DeltaE(elstat) depends on the size and shape of the occupied valence orbitals. The attractive electrostatic contributions in the diatomic molecules Li(2)-F(2) come from the s and p(sigma) electrons, while the p(pi) electrons do not compensate for nuclear-nuclear repulsion. It is the interplay of the three terms DeltaE(orb), DeltaE(Pauli), and DeltaE(elstat) that determines the bond energies and equilibrium distances of covalently bonded molecules. Molecules like N(2) and O(2), which are usually considered as covalently bonded, would not be bonded without the quasiclassical attraction DeltaE(elstat).     

On the viability of small endohedral hydrocarbon cage complexes: X@C4H4, X@C8H8, X@C8H14, X@C10H16, X@C12H12, AND X@C16H16

Small hydrocarbon complexes (X@cage) incorporating cage-centered endohedral atoms and ions (X = H(+), H, He, Ne, Ar, Li(0,+), Be(0,+,2+), Na(0,+), Mg(0,+,2+)) have been studied at the B3LYP/6-31G(d) hybrid HF/DFT level of theory. No tetrahedrane (C(4)H(4), T(d)()) endohedral complexes are minima, not even with the very small hydrogen atom or beryllium dication. Cubane (C(8)H(8), O(h)()) and bicyclo[2.2.2]octane (C(8)H(14), D(3)(h)()) minima are limited to encapsulating species smaller than Ne and Na(+). Despite its intermediate size, adamantane (C(10)H(16), T(d)()) can enclose a wide variety of endohedral atoms and ions including H, He, Ne, Li(0,+), Be(0,+,2+), Na(0,+), and Mg(2+). In contrast, the truncated tetrahedrane (C(12)H(12), T(d)()) encapsulates fewer species, while the D(4)(d)() symmetric C(16)H(16) hydrocarbon cage (see Table of Contents graphic) encapsulates all but the larger Be, Mg, and Mg(+) species. The host cages have more compact geometries when metal atoms, rather than cations, are inside. This is due to electron donation from the endohedral metals into C-C bonding and C-H antibonding cage molecular orbitals. The relative stabilities of endohedral minima are evaluated by comparing their energies (E(endo)) to the sum of their isolated components (E(inc) = E(endo) - E(cage) - E(x)) and to their exohedral isomer energies (E(isom) = E(endo) - E(exo)). Although exohedral binding is preferred to endohedral encapsulation without exception (i.e., E(isom) is always exothermic), Be(2+)@C(10)H(16) (T(d)(); -235.5 kcal/mol), Li(+)@C(12)H(12) (T(d)(); 50.2 kcal/mol), Be(2+)@C(12)H(12) (T(d)(); -181.2 kcal/mol), Mg(2+)@C(12)H(12) (T(d)(); -45.0 kcal/mol), Li(+)@C(16)H(16) (D(4)(d)(); 13.3 kcal/mol), Be(+)@C(16)H(16) (C(4)(v)(); 31.8 kcal/mol), Be(2+)@C(16)H(16) (D(4)(d)(); -239.2 kcal/mol), and Mg(2+)@C(16)H(16) (D(4)(d)(); -37.7 kcal/mol) are relatively stable as compared to experimentally known He@C(20)H(20) (I(h)()), which has an E(inc) = 37.9 kcal/mol and E(isom) = -35.4 kcal/mol. Overall, endohedral cage complexes with low parent cage strain energies, large cage internal cavity volumes, and a small, highly charged guest species are the most viable synthetic targets.     

Aromaticity in metallabenzenes

The electronic structure and bonding situation in 21 metallabenzenes (metal=Os, Ru, Ir, Rh, Pt, and Pd) were investigated at the DFT level (BP86/TZ2P) by using an energy decomposition analysis (EDA) of the interaction energy between various fragments. The aim of the work is to estimate the strength of the pi bonding and the aromatic character of the metallacyclic compounds. Analysis of the electronic structure shows that the metallacyclic moiety has five occupied pi orbitals, two with b1 symmetry and three with a2 symmetry, which describe the pi-bonding interactions. The metallabenzenes are thus 10 pi-electron systems. This holds for 16-electron and for 18-electron complexes. The pi bonding in the metallabenzenes results mainly from the b1 contribution, but the a2 contribution is not negligible. Comparison of the pi-bonding strength in the metallacyclic compounds with acylic reference molecules indicates that metallabenzenes should be considered as aromatic compounds whose extra stabilization due to aromatic conjugation is weaker than in benzene. The calculated aromatic stabilization energies (ASEs) are between 8.7 kcal mol(-1) for 13 and 37.6 kcal mol(-1) for 16 which is nearly as aromatic as benzene (ASE=42.5 kcal mol(-1)). The classical metallabenzene model compounds 1 and 4 exhibit intermediate aromaticity with ASE values of 33.4 and 17.6 kcal mol(-1). The greater stability of the 5d complexes compared with the 4d species appears not to be related to the strength of pi conjugation. From the data reported here there is no apparent trend or pattern which indicates a correlation between aromatic stabilization and particular ligands, metals, coordination numbers or charge. The lower metal-C5H5 binding energy of the 4d complexes correlates rather with weaker sigma-orbital interactions.     

Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5‐Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1:1 5‐hydroxytryptamine‐water (5‐HT‐H₂O) complexes are studied at the MP2/6‐311++G(d,p) level. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) analyses and the localized molecular orbital energy decomposition analysis (LMO‐EDA) were performed to explore the nature of the hydrogen‐bonding interactions in these complexes. Various types of hydrogen bonds (H‐bonds) are formed in these 5‐HT‐H₂O complexes. The intermolecular C4H5^5‐HT···O^w H‐bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5‐HT‐H₂O complexes. H‐bond in which nitrogen atom of amino in 5‐HT acted as proton donors was stronger than other H‐bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5‐HT‐H₂O complexes.     

Theoretical prediction of the heats of formation of C2H5O* radicals derived from ethanol and of the kinetics of beta-C-C scission in the ethoxy radical

Thermochemical parameters of three C(2)H(5)O* radicals derived from ethanol were reevaluated using coupled-cluster theory CCSD(T) calculations, with the aug-cc-pVnZ (n = D, T, Q) basis sets, that allow the CC energies to be extrapolated at the CBS limit. Theoretical results obtained for methanol and two CH(3)O* radicals were found to agree within +/-0.5 kcal/mol with the experiment values. A set of consistent values was determined for ethanol and its radicals: (a) heats of formation (298 K) DeltaHf(C(2)H(5)OH) = -56.4 +/- 0.8 kcal/mol (exptl: -56.21 +/- 0.12 kcal/mol), DeltaHf(CH(3)C*HOH) = -13.1 +/- 0.8 kcal/mol, DeltaHf(C*H(2)CH(2)OH) = -6.2 +/- 0.8 kcal/mol, and DeltaHf(CH(3)CH(2)O*) = -2.7 +/- 0.8 kcal/mol; (b) bond dissociation energies (BDEs) of ethanol (0 K) BDE(CH(3)CHOH-H) = 93.9 +/- 0.8 kcal/mol, BDE(CH(2)CH(2)OH-H) = 100.6 +/- 0.8 kcal/mol, and BDE(CH(3)CH(2)O-H) = 104.5 +/- 0.8 kcal/mol. The present results support the experimental ionization energies and electron affinities of the radicals, and appearance energy of (CH(3)CHOH+) cation. Beta-C-C bond scission in the ethoxy radical, CH(3)CH2O*, leading to the formation of C*H3 and CH(2)=O, is characterized by a C-C bond energy of 9.6 kcal/mol at 0 K, a zero-point-corrected energy barrier of E0++ = 17.2 kcal/mol, an activation energy of Ea = 18.0 kcal/mol and a high-pressure thermal rate coefficient of k(infinity)(298 K) = 3.9 s(-1), including a tunneling correction. The latter value is in excellent agreement with the value of 5.2 s(-1) from the most recent experimental kinetic data. Using RRKM theory, we obtain a general rate expression of k(T,p) = 1.26 x 10(9)p(0.793) exp(-15.5/RT) s(-1) in the temperature range (T) from 198 to 1998 K and pressure range (p) from 0.1 to 8360.1 Torr with N2 as the collision partners, where k(298 K, 760 Torr) = 2.7 s(-1), without tunneling and k = 3.2 s(-1) with the tunneling correction. Evidence is provided that heavy atom tunneling can play a role in the rate constant for beta-C-C bond scission in alkoxy radicals.     

Characteristics of antiaromatic ring pi multi-hydrogen bonds in (H2O)n-C4H4 (n = 1, 2) complexes

By counterpoise-corrected optimization method, the six antiaromatic ring pi multi-hydrogen bond structures with diversiform shapes for (H2O)n-C4H4 (n = 1,2) have been obtained at the MP2/aug-cc-pVDZ level. At the CCSD(T)/aug-cc-pVDZ level, the interaction energy obtained mainly depends on the numbers of H2O and fold numbers of the pi multi-hydrogen bond. The interaction energy order is -2.342 (1a with pi mono-hydrogen) < -2.777 (1b with pi bi-hydrogen) < -4.683 (2a with pi bi-hydrogen) < -4.734 (2b with pi tri-hydrogen) < -4.782 (2c with pi tri-hydrogen) < -5.009 kcal/mol (2d with pi tetra-hydrogen bond). Strangely, why is the interaction energy of the pi bi-hydrogen bond in 1b close to that of the pi mono-hydrogen bond in 1a (their difference is only 15.7%)? The reason is that a pi-type H-bond (as an accompanying interaction) between two lone pairs of the O-atom and a near pair of H-atoms of C4H4 exists shoulder by shoulder in structures 1a, 2a, 2b, and 2c and contributes to the interaction energy. Another accompanying interaction, a repulsive interaction between the pi H-bond (using the H-atom(s) of H2O) and the near pair of H-atoms of C4H4, is also found. For the structures and interaction energies, the pi-type H-bond produces four effects: bending the strong pi H-bond, attracting the pair of H-atoms of C4H4 so that they deviate from the C4 ring plane, showing the interaction energy contribution, and bringing the larger electron correlation contribution. The repulsive interaction also produces four effects: pushing the pair of H-atoms of C4H4 so that they deviate from its ring plane, elongating the distance of the pi H-bond, promoting the formation of pi-type H-bond, and slightly influencing the interaction energy. In the present paper, one C=C bond with two H2O (over and below the ring plane) forms a pi H-bond link in two ways: a strong-weak pi H-bond link and a strong-strong pi H-bond link. The stability contribution of the former is more favorable than the latter. One H2O forms a pi H-bond with C4H4 in two ways. One strong pi H-bond part (over or below the ring plane) always is accompanied by another H-bond part. The accompanying part is either a weak pi H-bond or pi-type H-bond.     

Theoretical investigation of carbon-sulfur triple bonds

By means of first principle calculations we have investigated a set of molecules that are presumed to contain carbon-sulfur triple bonds, namely HCSOH, H(3)SCH, cis-FCSF, F(3)CCSF(3), and F(5)SCSF(3). For HCSOH, FCSF, and H(3)SCH we used the CCSD(T) methodology and the correlation-consistent basis sets. On the other hand, F(3)CCSF(3) and F(5)SCSF(3) were studied at the B3LYP, M06-2X, MP2, and G3 levels of theory. We found that none of these molecules display a carbon-sulfur adiabatic bond dissociation energy (ABDE) as strong as diatomic CS (170.5 kcal mol(-1)), or a diabatic bond dissociation energy (DBDE) larger than the one found in SCO (212.0 kcal mol(-1)), although the DBDE of FCSF comes quite close at 208.3 kcal mol(-1). The CS ABDEs of F(3)CCSF(3), F(5)SCSF(3), and H(3)SCH are comparable to that of a single C-S bond. In contrast with the experimental results, F(3)CCSF(3) and F(5)SCSF(3) are predicted to be linear with C(3v) and C(s) symmetry, respectively, at the B3LYP/6-311+G(3df,2p) level. MP2/6-311+G(2df,2p) calculations support the C(3v) symmetry for F(3)CCSF(3), despite F(5)SCSF(3) not having a perfect linear structure; the CSC angle is 174.6°, which is nearly 20° larger than the experimental value. The analysis of the carbene structures of HCSOH and H(3)SCH revealed that they are not significant, because the triplet state is dissociative in these cases. However, for F(3)CCSF(3) and F(5)SCSF(3) , the carbene triplet states lie 0.81 and 0.77 eV above the singlet state, respectively. In the same vein, our investigation supports the presence of a strong double bond for HCSOH. The conflicting evidence available for F(3)CCSF(3) and F(5)SCSF(3) makes it very difficult to determine the nature of the CS bonds. However, the bond dissociation energies and the singlet-triplet splittings clearly suggest that these compounds should be considered as masked sulfinylcarbenes. The analysis of the bond dissociation energies challenges the existence of a triple bond in these five molecules, but from a strictly thermodynamic standpoint, cis-FCSF is found to be the candidate most likely to exhibit triple-bond character.     

Electronic structure and chain-length effects in diplatinum polyynediyl complexes trans,trans-[(X)(R3P)2Pt(C triple bond C)(n)Pt(PR3)2(X)]: a computational investigation

Structure and bonding in the title complexes are studied using model compounds trans,trans-[(C6H5)(H3P)2Pt(C triple bond C)(n)Pt(PH3)2(C6H5)] (PtCxPt; x = 2n = 4-26) at the B3LYP/LACVP* level of density functional theory. Conformations in which the platinum square planes are parallel are very slightly more stable than those in which they are perpendicular (DeltaE = 0.12 kcal mol(-1) for PtC8Pt). As the carbon-chain length increases, progressively longer C triple bond C triple bonds and shorter triple bond C-C triple bond single bonds are found. Whereas the triple bonds in HCxH become longer (and the single bonds shorter) as the interior of the chain is approached, the PtC triple bond C triple bonds in PtCxPt are longer than the neighboring triple bond. Also, the Pt-C bonds are shorter at longer chain lengths, but not the H-C bonds. Accordingly, natural bond orbital charge distributions show that the platinum atoms become more positively charged, and the carbon chain more negatively charged, as the chain is lengthened. Furthermore, the negative charge is localized at the two terminal C triple bond C atoms, elongating this triple bond. Charge decomposition analyses show no significant d-pi* backbonding. The HOMOs of PtCxPt can be viewed as antibonding combinations of the highest occupied pi orbital of the sp-carbon chain and filled in-plane platinum d orbitals. The platinum character is roughly proportional to the Pt/Cx/Pt composition (e.g., x = 4, 31 %; x = 20, 6 %). The HOMO and LUMO energies monotonically decrease with chain length, the latter somewhat more rapidly so that the HOMO-LUMO gap also decreases. In contrast, the HOMO energies of HCxH increase with chain length; the origin of this dichotomy is analyzed. The electronic spectra of PtC4Pt to PtC10Pt are simulated. These consist of two pi-pi* bands that redshift with increasing chain length and are closely paralleled by real systems. A finite HOMO-LUMO gap is predicted for PtCinfinityPt. The structures of PtCxPt are not strictly linear (average bond angles 179.7 degrees -178.8 degrees ), and the carbon chains give low-frequency fundamental vibrations (x = 4, 146 cm(-1); x = 26, 4 cm(-1)). When the bond angles in PtC12Pt are constrained to 174 degrees in a bow conformation, similar to a crystal structure, the energy increase is only 2 kcal mol(-1). The above conclusions should extrapolate to (C triple bond C)(n) systems with other metal endgroups.     

Hyperconjugative stabilization in alkyl carbocations: direct estimate of the beta-effect of group-14 elements

DFT calculations at the BP86/TZ2P level have been carried out for the primary, secondary, and tertiary carbenium ions [H2C-CH(EH3)2](+) (1a-e), [HC{CH(EH3)2}2](+), (2a-e), and [C{CH(EH3)2}3](+) (3a-e) for E = C, Si, Ge, Sn, Pb. The nature of the interaction between the carbenium center H(2-n)C(+) and the substituents {CH(EH3)2}m has been investigated with an energy decomposition analysis (EDA) aiming at estimating the strength of the pi hyperconjugation which electronically stabilizes the carbenium ions. The results of the EDA show that the calculated DeltaEpi values can be used as a measure for the strength of hyperconjugation in carbenium ions arising from the interactions of saturated groups possessing pi orbitals. The theoretical data suggest that the ability of sigma C-E bonds to stabilize positive charges by hyperconjugation follow the order C < Si < Ge < Sn < Pb. Hyperconjugation of C-Si bonds is much stronger than hyperconjugation of C-C bonds while the further rising from silicon to lead is smaller and has about the same step size for each element. The strength of the hyperconjugation in primary, secondary, and tertiary alkyl carbenium ions does not increase linearly with the number of hyperconjugating groups; the incremental stabilization becomes smaller from primary to secondary to tertiary cations. The effect of hyperconjugation is reflected in the shortening of the C-C bond distances and in the lengthening of the C-E bonds, which exhibits a highly linear relationship between the calculated C-C and C-E distances in carbocations 1-3 and the hyperconjugation estimated by the DeltaEpi values.     

Computational study of analogues of the uranyl ion containing the -N=U=N- unit: density functional theory calculations on UO2(2+), UON+, UN2, UO(NPH3)3+, U(NPH3)2(4+), [UCl4[NPR3]2] (R = H, Me), and [UOCl4[NP(C6H5)3]

The electronic and geometric structures of the title species have been studied computationally using quasi-relativistic gradient-corrected density functional theory. The valence molecular orbital ordering of UO2(2+) is found to be pi g < pi u < sigma g < sigma u (highest occupied orbital), in agreement with previous experimental conclusions. The significant energy gap between the sigma g and sigma u orbitals is traced to the "pushing from below" mechanism: a filled-filled interaction between the semi-core uranium 6p atomic orbitals and the sigma u valence level. The U-N bonding in UON+ and UN2 is significantly more covalent than the U-O bonding in UON+ and UO2(2+). UO(NPH3)3+ and U(NPH3)2(4+) are similar to UO2(2+), UON+, and UN2 in having two valence molecular orbitals of metal-ligand sigma character and two of pi character, although they have additional orbitals not present in the triatomic systems, and the U-N sigma levels are more stable than the U-N pi orbitals. The inversion of U-N sigma/pi orbital ordering is traced to significant N-P (and P-H) sigma character in the U-N sigma levels. The pushing from below mechanism is found to destabilize the U-N f sigma molecular orbital with respect to the U-N d sigma level in U(NPH3)2(4+). The uranium f atomic orbitals play a greater role in metal-ligand bonding in UO2(2+), UN2, and U(NPH3)2(4+) than do the d atomic orbitals, although, while the relative roles of the uranium d and f atomic orbitals are similar in UO2(2+) and U(NPH3)2(4+), the metal d atomic orbitals have a more important role in the bonding in UN2. The preferred UNP angle in [UCl4(NPR3)2] (R = H, Me) and [UOCl4(NP(C6H5)3)]- is found to be close to 180 degrees in all cases. This preference for linearity decreases in the order R = Ph > R = Me > R = H and is traced to steric effects which in all cases overcome an electronic preference for bending at the nitrogen atom. Comparison of the present iminato (UNPR3) calculations with previous extended Hückel work on d block imido (MNR) systems reveals that in all cases there is little or no preference for linearity over bending at the nitrogen when R is (a) only sigma-bound to the nitrogen and (b) sterically unhindered. The U/N bond order in iminato complexes is best described as 3.     

Spectroscopic detection and theoretical confirmation of the role of Cr2(CO)5(C5R5)2 and .Cr(CO)2(ketene)(C5R5) as intermediates in carbonylation of N=N=CHSiMe3 to O=C=CHSiMe3 by .Cr(CO)3(C5R5) (R = H, CH3)

Conversion of N=N=CHSiMe3 to O=C=CHSiMe3 by the radical complexes .Cr(CO)3C5R5 (R = H, CH3) derived from dissociation of [Cr(CO)3(C5R5)]2 have been investigated under CO, Ar, and N2 atmospheres. Under an Ar or N2 atmosphere the reaction is stoichiometric and produces the Cr[triple bond]Cr triply bonded complex [Cr(CO)2(C5R5)]2. Under a CO atmosphere regeneration of [Cr(CO)3(C5R5)]2 (R = H, CH3) occurs competitively and conversion of diazo to ketene occurs catalytically as well as stoichiometrically. Two key intermediates in the reaction, .Cr(CO)2(ketene)(C5R5) and Cr2(CO)5(C5R5)2 have been detected spectroscopically. The complex .Cr(13CO)2(O=13C=CHSiMe3)(C5Me5) has been studied by electron spin resonance spectroscopy in toluene solution: g(iso) = 2.007; A(53Cr) = 125 MHz; A(13CO) = 22.5 MHz; A(O=13C=CHSiMe3) = 12.0 MHz. The complex Cr2(CO)5(C5H5)2, generated in situ, does not show a signal in its 1H NMR and reacts relatively slowly with CO. It is proposed to be a ground-state triplet in keeping with predictions based on high level density functional theory (DFT) studies. Computed vibrational frequencies are also in good agreement with experimental data. The rates of CO loss from 3Cr2(CO)5(C5H5)2 producing 1[Cr(CO)2(C5H5)]2 and CO addition to 3Cr2(CO)5(C5H5)2 producing 1[Cr(CO)3(C5H5)]2 have been measured by kinetics and show DeltaH approximately equal 23 kcal mol(-1) for both processes. Enthalpies of reduction by Na/Hg under CO atmosphere of [Cr(CO)n(C5H5)]2 (n = 2,3) have been measured by solution calorimetry and provide data for estimation of the Cr[triple bond]Cr bond strength in [Cr(CO)2(C5H5)]2 as 72 kcal mol(-1). The complex [Cr(CO)2(C5H5)]2 does not readily undergo 13CO exchange at room temperature or 50 degrees C implying that 3Cr2(CO)5(C5H5)2 is not readily accessed from the thermodynamically stable complex [Cr(CO)2(C5H5)]2. A detailed mechanism for metalloradical based conversion of diazo and CO to ketene and N2 is proposed on the basis of a combination of experimental and theoretical data.