How the addition of atomic hydrogen to a multiple bond can be catalyzed by water molecules.

Observational data show complex organic molecules in the interstellar medium (ISM). Hydrogenation of small unsaturated carbon double bond could be one way for molecular complexification. It is important to understand how such reactivity occurs in the very cold and low-pressure ISM. Yet, there is water ice in the ISM, either as grain or as mantle around grains. Therefore, the addition of atomic hydrogen on double-bonded carbon in a series of seven molecules have been studied and it was found that water catalyzes this reaction. The origin of the catalysis is a weak charge transfer between the π MO of the unsaturated molecule and H atom, allowing a stabilizing interaction with H2O. This mechanism is rationalized using the non-covalent interaction and the quantum theory of atoms in molecules approaches.

The significant role of water in reactions occurring on the surface of interstellar ice grains: Hydrogenation of pure ketene H2CCO ice versus hydrogenation of mixed H2CCO/H2O ice at 10 K.

Water ice plays an important role in reactions taking place on the surface of interstellar ice grains, ranging from catalytic effects that reduce reaction barrier heights to effects that stabilize the reaction products and intermediates formed, or that favor one reaction pathway over another, passing through water-involvement in the reaction to produce more complex molecules that cannot be formed without water or water-derived fragments H, O and OH. In this context, we have combined experimental and theoretical studies to investigate ketene (CH2CO) + H solid-state reaction at 10 K in the presence and absence of water molecules under interstellar conditions, through H-bombardment of CH2CO and CH2CO/H2O ices. We show in the present study that with or without water, the ketene molecule reacts with H atoms to form four reaction products, namely CO, H2CO, CH4 and CH3CHO. Based on the amounts of CH2CO consumed during the hydrogenation processes, the CH2CO + 2H reaction appears to be more efficient in the presence of water. This underlines the catalytic role of water ice in reactions occurring on the surface of interstellar ice grains. However, if we refer to the yields of reaction products formed during the hydrogenation of CH2CO and CH2CO/H2O ices, we find that water molecules favor the reaction pathway to form CH3CHO and deactivate that leading to CH4 and H2CO. These experimental results are in good agreements with the theoretical predictions that highlight the catalytic effect of H2O on the CH2CO + H reaction, whose potential energy barrier drops from 4.6 kcal mol-1 (without water) to 3.8 and 3.6 kcal mol-1 with one and two water molecules respectively.

Conformational preference analysis in C2H6 using orbital forces and non-covalent interactions; comparison with related systems.

Dynamic Orbital Forces (DOF) and Non-Covalent Interactions (NCIs) are used to analyze the attractive/repulsive interactions responsible of the conformational preference of ethane and some related compounds. In ethane, it is found that the stabilization of the staggered conformation with respect to the adiabatic eclipsed one arises from both attractive and repulsive interactions in CH3⋯CH3. Attractive ones are predominant in a ratio 2 : 1, with an important role of a σ MO. On the contrary, the stabilization of the staggered conformation with respect to the vertical eclipsed one arises almost only from repulsive π interactions. Weak long-range H⋯H repulsions also favour the staggered conformation. From the sum of DOFs, yielding intrinsic bond energies, the rotation barrier can be decomposed into a weakening of the C-C bond (ca. 7 kcal mol-1), moderated by a strengthening of C-H ones (ca. 4 kcal mol-1). This evidences the decrease of hyperconjugation in the eclipsed conformation with respect to the staggered one. In the compounds CH3-SiH3, SiH3-SiH3, CH3-CF3 and CF3-CF3, the conformational preference is predominantly or exclusively due to repulsive interactions, with respect as well to adiabatic as to vertical eclipsed structures.

Formation of CO, CH4, H2CO and CH3CHO through the H2CCO + H surface reaction under interstellar conditions.

The reaction of ketene (H2CCO) with hydrogen atoms has been studied under interstellar conditions through two different experimental methods, occurring on the surface and in the bulk of H2CCO ice. We show that ketene interaction with H-atoms at 10 K leads mainly to four reaction products, carbon monoxide (CO), methane (CH4), formaldehyde (H2CO) and acetaldehyde (CH3CHO). A part of these results shows a chemical link between a simple organic molecule such as H2CCO and a complex one such as CH3CHO, through H-addition reactions taking place in dense molecular clouds. The H-addition processes are very often proposed by astrophysical models as mechanisms for the formation of complex organic molecules based on the abundance of species already detected in the interstellar medium. However, the present study shows that the hydrogenation of ketene under non-energetic conditions may also lead efficiently to fragmentation processes and the formation of small species such as CO, CH4 and H2CO, without supplying external energy such as UV photons or high energy particles. Such fragmentation pathways should be included in the astrophysical modeling of H2CCO + H in the molecular clouds of the interstellar medium. To support these results, theoretical calculations have explicitly showed that, under our experimental conditions, H-atom interactions with the CC bond of ketene lead mainly to CH3CHO, CH4 and CO. By investigating the formation and reactivity of the reaction intermediate H3C-CO radical, our calculations demonstrate that the H3C-CO + H reaction evolves through two barrierless pathways to form either CH3CHO or CH4 and CO fragments.

Orbital energies and nuclear forces in DFT: Interpretation and validation.

The bonding and antibonding character of individual molecular orbitals has been previously shown to be related to their orbital energy derivatives with respect to nuclear coordinates, known as dynamical orbital forces. Albeit usually derived from Koopmans' theorem, in this work we show a more general derivation from conceptual DFT, which justifies application in a broader context. The consistency of the approach is validated numerically for valence orbitals in Kohn-Sham DFT. Then, we illustrate its usefulness by showcasing applications in aromatic and antiaromatic systems and in excited state chemistry. Overall, dynamical orbital forces can be used to interpret the results of routine ab initio calculations, be it wavefunction or density based, in terms of forces and occupations.

Ligand effects on coordination properties of organolithium compounds: insights from computational experiments on a "weakened" Li.

A simple numerical experiment is presented which allows tuning the lithium electrophilicity, a parameter strongly affected by the solvent and additives coordination. A series of coordination of Li+ to carbanions or polydentate nucleophiles is examined showing the potential and the limits of this approach. The results suggest that such a simple trick can be remarkably helpful to model and decipher the effects of solvation on the structure and properties of lithiated organometallic species.

The "Inverted Bonds" Revisited: Analysis of "In Silico" Models and of [1.1.1]Propellane by Using Orbital Forces.

This article dwells on the nature of "inverted bonds", which refer to the σ interaction between two sp hybrids by their smaller lobes, and their presence in [1.1.1]propellane. Firstly, we study H3 C-C models of C-C bonds with frozen H-C-C angles reproducing the constraints of various degrees of "inversion". Secondly, the molecular orbital (MO) properties of [1.1.1]propellane and [1.1.1]bicyclopentane are analyzed with the help of orbital forces as a criterion of bonding/antibonding character and as a basis to evaluate bond energies. Triplet and cationic states of [1.1.1]propellane species are also considered to confirm the bonding/antibonding character of MOs in the parent molecule. These approaches show an essentially non-bonding character of the σ central C-C interaction in propellane. Within the MO theory, this bonding is thus only due to π-type MOs (also called "banana" MOs or "bridge" MOs) and its total energy is evaluated to approximately 50 kcal mol-1 . In bicyclopentane, despite a strong σ-type repulsion, a weak bonding (15-20 kcal mol-1 ) exists between both central C-C bonds, also due to π-type interactions, though no bond is present in the Lewis structure. Overall, the so-called "inverted" bond, as resulting from a σ overlap of the two sp hybrids by their smaller lobes, appears highly questionable.

Analysis of Reaction Processes On the Basis of the Evolution of Dynamic Orbital Forces: Examples of Cycloadditions, SN 2 Substitution, Nucleophilic Addition, and Hydrogen Transposition.

The derivative of the energy of a canonical molecular orbital (MO) [or dynamical orbital forces (DOFs)] with respect to a bond length provides a reliable index of the bonding/antibonding character of this MO on this bond. The DOFs of selected MOs as a function of the reaction coordinate were computed for a panel of model reaction mechanisms: [2+4] (Diels-Alder) cycloaddition, [2+2] cycloaddition, second-order nucleophilic substitution (SN 2), nucleophilic addition to a carbonyl group, and [1,2] hydrogen transposition. The results highlight the nature of the reorganization of the main MOs and the stage of the reaction coordinate (RC) at which it occurs. For instance, in the Diels-Alder reaction, one can identify a part of the reaction that is dominated by repulsive four-electron interactions and another part dominated by attractive two-electron interactions. Also, the shape of the DOF as a function of the reaction coordinate reveals the existence of avoided MO crossings and their location on the RC. Even for spontaneous reactions with monotonic variation in the potential energy, extrema of the MO energy and sudden electron rearrangements can be put into evidence. This study provides quantitative support to classical MO analyses of reactivity such as correlation diagrams and frontier approximation.

Organohelium compounds: structures, stabilities and chemical bonding analyses.

This paper deals with the possibility of forming short and relatively strong carbon-helium bonds in small typical organic molecules through substitution of one or several H atoms by He(+). A structural and energetics study (based on high-level calculations) of this unusual bonding, as well as a topological characterization of the resulting cations, is undertaken. Stable species generally requires substitution of about half of the hydrogen atoms for formation. Under these conditions, the number of such species appears to be potentially unlimited. "True" C-He bonds exhibit equilibrium distances ranging from 1.327 (C2H2He2(2+)) to 1.129 Å (He2CO(2+)). The energies of neutral He releasing range from approximately 5 kcal mol(-1) [He2CO(2+), (Z)-C2H2He2(2+)] to 25 kcal mol(-1) (C2HHe3(3+)), but remain most frequently around 10 kcal mol(-1). However, most of He(+)-substituted hydrocarbons are metastable with respect to C-C cleavage, except derivatives of ethene. Atoms in molecules (AIM) and electron localization function (ELF) topological descriptors classify the C-He bond as a weak charge-shift interaction [S. Shaik, D. Danovich, B. Silvi, D. L. Lauvergnat, P. C. Hiberty, Chem. Eur. J. 2005, 11, 6358-6371] in agreement with a recent publication by Rzepa [S. H. Rzepa, Nat. Chem. 2010, 2, 390-393]. He2CO(2+) is the only investigated compound that presents a C-He bonding ELF basin, which indicates a non-negligible covalent contribution to the bond. Other modifications in the electronic structure, such as the breaking of the triple bond in ethyne derivatives or the loss of aromaticity in C6H3He3(3+), are also nicely revealed by the ELF topology.

Theoretical design of strong neutral radical-boron adducts: trisubstituted boranes as potential radical scavengers.

The conditions of formation of strong two-center one-electron bonds in neutral compounds are discussed. Both molecular orbital and valence bond analyses show that good candidates are adducts of radicals ˙AR3 (A=C, Si, Ge) of low ionization energy (IE) with boranes BX3 of high electron affinity (EA). This is confirmed by ab initio calculations. The bond energy of adducts of B(CF3)3 with various radicals ranges from 18 kcal mol(-1) for ˙CH3 to approximately 40 kcal mol(-1) for Me3Si˙, and a clear correlation with IE-EA difference is found. This allows one to expect B(CF3)3, among other fluoroboranes, to be an efficient radical scavenger.

Inductive effects on proton affinity of benzene derivatives: analysis using fictitious hydrogen atoms.

Pure inductive effects on the gas-phase basicity of seven benzene derivatives (3- and 4-substitution) are monitored in a continuous way using fictitious hydrogen atoms bearing an adjustable nuclear charge Z*. This approach (H* method) affords three main advantages over existing treatments: such entities are by definition purely inductive (without any underlying assumptions), use of empirical parameters is circumvented, and yet the method has been designed to remain particularly easy to use. We directly establish the linear dependence of proton affinities on inductive effects and, more quantitatively, measure accurate sensitivities rho(I)* analogous to Taft's coefficients. Functional centers exhibit contrasted values, up to a factor of 3, which finds an interpretation within the framework of the HSAB theory. The sensitivities rho(I)* for 3- and 4-substitution are quantified. The associated para/meta rho(I)* ratio ranges from 1.02 to 1.16 according to the functional center. These values, always slightly superior to unity, denote a contribution of pi electrons in the transmission of the inductive effect. This effect, first identified by Exner, is shown to account for ca. 30% of the basicity of benzoic acid, which is taken as an example.

Electron pair localization function: a practical tool to visualize electron localization in molecules from quantum Monte Carlo data.

In this work we introduce an electron localization function describing the pairing of electrons in a molecular system. This function, called "electron pair localization function," is constructed to be particularly simple to evaluate within a quantum Monte Carlo framework. Two major advantages of this function are the following: (i) the simplicity and generality of its definition; and (ii) the possibility of calculating it with quantum Monte Carlo at various levels of accuracy (Hartree-Fock, multiconfigurational wave functions, valence bond, density functional theory, variational Monte Carlo with explicitly correlated trial wave functions, fixed-node diffusion Monte Carlo, etc). A number of applications of the electron pair localization function to simple atomic and molecular systems are presented and systematic comparisons with the more standard electron localization function of Becke and Edgecombe are done. Results illustrate that the electron pair localization function is a simple and practical tool for visualizing electronic localization in molecular systems.