Energy and spectroscopic parameters of neutral and cations isomers of the CnH2 (n = 2-6) families using high-level ab-initio approaches.

Cationic species, previously detected from ion-induced desorption of solid methane by plasma desorption mass spectrometry (PDMS), and neutral species, are investigated using high-level ab-initio approaches. From a set of 25 cationic and 26 neutral structures belonging to CnH2 (n = 2-6) families, it was obtained the energy, rotational constants, harmonic vibrational frequency, charge distribution and excitation energies. The ZPVE-corrected energies, at CCSD(T)-F12; CCSD(T)-F12/RI/(cc-pVTZ-F12, cc-pVTZ-F12-CABS, cc-pVQZ/C) (n = 2-5) and CCSD(T)/cc-pVTZ (n = 6) levels, reveal that the topology of the most stable isomer vary with n and the charge. Out of 674 harmonic frequencies, those with maximum intensity are generally in the 3000-3500 cm-1 range. Analysis of 169 vertical transition energies calculated with the EOM-CCSD approach, suggest three C6H2 species as potential carriers of the diffuse interstellar bands (DIB). Systematic comparison of properties between neutral and cationic species can assist in the structural description of complex matrices.

Three-centre two-electron bonds from the quantum interference perspective.

The nature of the three-center two-electron (3c2e) chemical bond is investigated by the Interference Energy Analysis (IEA) method and using a SC(2, 3) (spin coupled wave function with two electrons and three orbitals) approach for describing 3c2e bonds. In this approach, each center involved in the bonding contributes with a one-electron state for the interference process. The species H3+, Li3+, LiH2+, C3H5+, C3H3+, R2CBeCR2 (R = H, CH3), C7H11+ and CH5+ are considered in the study. The results show that 3c2e bonds have the same features of a 2c2e bond: the stability of the chemical systems exhibiting 3c2e bonds derives from quantum interference among electronic states. Other (quasi-classical) factors are always overall destabilizing, mostly because of the nuclear repulsion. The interference energy of a 3c2e bond is about three times higher than that of a 2c2e bond involving atoms of the same elements. In particular, concerning Li3+ and C3H3+ we found no evidence that the 'aromatic' character attributed to those species imparts any special features to their chemical structures, compared to other 3c2e bonds. Therefore, these species exhibit multicenter bonds, essentially equivalent to the other studied cases.

Supported or unsupported three-center two-electron bonds? A criterion based on Interference Energy Analysis.

The classification of three-center two-electron (3c2e) bonds into supported (closed) or unsupported (open) was proposed by Lipscomb in his work on boranes and extended to transition metal complexes by Bau and co-workers. The species in which the interactions of the terminal atoms are negligible are called "unsupported bonds." Examples of chemical species that are said to exhibit such bonds are Li2H+, Na2H+, B2H7 -, Al2(CH3)7 -, and [(µ2-H)Cr2(CO)10]- although the general criterion for distinguishing these types of bonds is somewhat qualitative. Besides providing a unifying view of the nature of the chemical bond, in terms of quantum interference among electronic states, the Generalized Product Function Energy Partitioning method through the Interference Energy Analysis (IEA) is also potentially capable of providing a rigorous ground to the concept of supported bonds by looking at the specific interference energies between the orbital pairs associated with the bond. The IEA was performed in the species Li2H+, Na2H+, B2H7 -, C2H7 -, Al2H7 -, and [(µ2-H)Cr2(CO)10]-, as well as along the reaction path Li2H+ → Li2 + + H. The results shown that in all studied A-B-C bonds, the A-C interactions are as important as the A-B/B-C ones, leading to the conclusion that all studied 3c2e bonds are "supported," in the sense that the A-C interaction is not negligible. The particularity of those species in preferring linear geometry is completely explained by quasi-classical effects, more specifically, by minimization of the electron-electron and nucleus-nucleus repulsions.

Substituent Effects on the Quantum Interference of Two-Center One-Electron Bonds: [B2X6]- (X = H, F, Cl, CN, OH, CH3, and OCH3).

The interference energy analysis (IEA) provided by the generalized product function energy partitioning (GPF-EP) method was applied to investigate the influence of the neighboring atoms on the nature of the two-center one-electron (2c1e) bonds in the anion dimers of BX3 species (X = H, F, Cl, CN, OH, CH3, and OCH3). The species were studied at the GVB-PP(6/12).SC(1,2)/6-31**G++ level of calculation. The IEA has revealed that there is a balance between two main factors determining the chemical stability of the species. Quantum interference acts as the sole stabilizing effect in the formation of the chemical bonds, particularly as the result of the drop in kinetic energy, and the electronegativity of the substituent has a direct influence on the magnitude of this effect. The quasi-classical energy is responsible for the destabilizing factors, mainly the group bulkiness, and the "electron-withdrawing" effect in the case of the cyano group.

The Valence-Bond (VB) Model and Its Intimate Relationship to the Symmetric or Permutation Group.

VB and molecular orbital (MO) models are normally distinguished by the fact the first looks at molecules as a collection of atoms held together by chemical bonds while the latter adopts the view that each molecule should be regarded as an independent entity built up of electrons and nuclei and characterized by its molecular structure. Nevertheless, there is a much more fundamental difference between these two models which is only revealed when the symmetries of the many-electron Hamiltonian are fully taken into account: while the VB and MO wave functions exhibit the point-group symmetry, whenever present in the many-electron Hamiltonian, only VB wave functions exhibit the permutation symmetry, which is always present in the many-electron Hamiltonian. Practically all the conflicts among the practitioners of the two models can be traced down to the lack of permutation symmetry in the MO wave functions. Moreover, when examined from the permutation group perspective, it becomes clear that the concepts introduced by Pauling to deal with molecules can be equally applied to the study of the atomic structure. In other words, as strange as it may sound, VB can be extended to the study of atoms and, therefore, is a much more general model than MO.

Mechanistic Insights into the Formation of Lithium Fluoride Nanotubes.

Born-Oppenheimer molecular dynamics (BOMD) and periodic density functional theory (DFT) calculations have been applied for describing the mechanism of formation of lithium fluoride (LiF) nanotubes with cubic, hexagonal, octagonal, decagonal, dodecagonal, and tetradecagonal cross-sections. It has been shown that high energy structures, such as nanowires, nanorings, nanosheets, and nanopolyhedra are transient species for the formation of stable nanotubes. Unprecedented (LiF)n clusters (n≤12) were also identified, some of them lying less than 10 kcal mol-[1] above their respective global minima. Such findings indicate that stochastic synthetic techniques, such as laser ablation and chemical vapor deposition, should be combined with a template-driven procedure in order to generate the nanotubes with adequate efficiency. Apart from the stepwise growth of LiF units, the formation of nanotubes was also studied by rolling up a planar square sheet monolayer, which could be hypothetically produced from the exfoliation of the FCC crystal structure. It was shown that both pathways could lead to the formation of alkali halide nanotubes, a still unprecedented set of one-dimensional materials.

One-electron bonds are not "half-bonds".

Despite the success of the molecular orbital (MO) and valence-bond (VB) models to describe the electronic structure and properties of molecules, neither MO nor VB provides an explanation for the nature of the chemical bond. The first to address this problem was Ruedenberg, who showed that chemical bonds result from quantum interference. He developed a method to calculate the interference contribution to the total electronic energy and density and applied it to molecules containing typical two-centre two-electron (2c-2e) covalent bonds. To test the generality of Ruedenberg's hypothesis, we developed a powerful Interference Energy Analysis (IEA) method to calculate the interference contributions of individual chemical bonds to the total energy of diatomic and polyatomic molecules, and showed that any two-electron bond, despite its polarity, results from quantum interference. Nevertheless, many stable molecules are experimentally known whose chemical structures clearly indicate the existence of two-centre one-electron bonds (2c-1e). Therefore, the question remains if quantum interference will be the dominant effect for these systems. This work describes the extension of the IEA for treating two-centre one-electron bonds, making use of a Generalised Product Function (GPF) built from spin coupled wave functions of N electrons in M orbitals, SC(N,M). Several diatomic and polyatomic molecules were analysed and whenever possible the results were compared with the analogous case of a two-electron bond. The results indicate that interference is the dominant effect for the one-electron bonds, which reinforces the role of quantum interference as the central element in chemical bonding theory.

Unexpected reversal of stability in strained systems containing one-electron bonds.

Ring strain energy is a very well documented feature of neutral cycloalkanes, and influences their structural, thermochemical and reactivity properties. In this work, we apply density functional theory and high-level coupled cluster calculations to describe the geometry and relative stability of C6H12+˙ radical cations, whose cyclic isomers are prototypes of singly-charged cycloalkanes. Molecular ions with the mentioned stoichiometry were produced via electron impact experiments using a gaseous cyclohexane sample (20-2000 eV). From our calculations, in addition to structures that resemble linear and branched alkenes as well as distinct conformers of cyclohexane, we have found low-lying species containing three-, four- and five-membered rings with the presence of an elongated C-C bond. Remarkably, the stability trend of these ring-bearing radical cations is anomalous, and the three-membered species are up to 11.3 kcal mol-1 more stable than the six-membered chair structure. Generalized Valence Bond calculations and the Spin Coupled theory with N electrons and M orbitals were used in conjunction with the Generalized Product Function Energy Partitioning (GPF-EP) method and Interference Energy Analysis (IEA) to describe the chemical bonding in such moieties. Our results confirm that these elongated C-C motifs are one-electron sigma bonds. Our calculations also reveal the effects that drive thermochemical preference of strained systems over their strained-free isomers, and the origin of the unusual stability trend observed for cycloalkane radical cations.

Soft X-ray Chlorine Photolysis on Chlorobenzene Ice: An Experimental and Theoretical Study.

An experimental and theoretical study of the photoinduced homolysis of the carbon-chlorine bond in an ice matrix of chlorobenzene is presented. A condensed chlorobenzene film has been grown in situ and near edge X-ray fine structure (NEXAFS) spectra were collected after exposing the condensed film to a monochromatic photon beam centered at the 2822 eV resonant excitation of chlorine and at 2850 eV. The photoabsorption to the Cl 1s → σ* and Cl 1s → π* states has been measured and the hypothesis of free radical coupling reactions was investigated via time-dependent density functional theory (TD-DFT) and complete active space self-consistent field (CASSCF) calculations. Also, potential energy pathways to the C-Cl cleavage have been obtained at the CASSCF level to the Cl 1s → σ*, 1s → π*, and 1s → ∞ states. A strong dissociative character was only found for the Cl 1s → σ* resonance.

Theoretical study of the absolute inner-shell photoionization cross sections of the formic acid and some of its hydrogen-bonded clusters.

Inner-shell absolute photoabsorption and photoionization cross sections of the formic acid, HCOOH, and its small hydrogen-bonded clusters, i.e., (HCOOH)2, HCOOH2 +, HCOHOH+, and HCOOH·H3O+, were calculated at the time-dependent density functional theory (TDDFT) level, and the results were used to analyze the effect of the formic acid clustering on the carbon and oxygen K-edge photoionization cross sections. The discrete electronic pseudospectra obtained with square-integrable (L2) basis set calculations were used in an analytic continuation procedure based on continued fraction functions to obtain the photoabsorption cross sections. Symmetry adapted cluster configuration interaction calculations on the small formic acid clusters have also been performed at the oxygen K-edge to assign the discrete transitions and ionization potentials in support to the TDDFT results.

Molecular inner-shell photoabsorption/photoionization cross sections at core-valence-separated coupled cluster level: Theory and examples.

Oxygen, nitrogen, and carbon K-shell photoabsorption and photoionization cross sections have been calculated within core-valence-separated coupled cluster (CC) linear response theory for a number of molecular systems, namely, water, ammonia, ethylene, carbon dioxide, acetaldehyde, furan, and pyrrole. The cross sections below and above the K-edge core ionization thresholds were obtained, on the same footing, from L2 basis set calculations of the discrete electronic pseudospectrum yielded by an asymmetric-Lanczos-based formulation of CC linear response theory at the CC singles and doubles (CCSD) and CC singles and approximate doubles (CC2) levels. An analytic continuation procedure for both discrete and continuum cross sections as well as a Stieltjes imaging procedure for the photoionization cross section were applied and the results critically compared.

Doubly and Triply Charged Species Formed from Chlorobenzene Reveal Unusual C-Cl Multiple Bonding.

In free-radical halogenation of aromatics, singly charged ions are usually formed as intermediates. These stable species can be easily observed by time-of-flight mass spectrometry (TOF-MS). Here we used electron and proton beams to ionize chlorobenzene (C6H5Cl) and investigate the ions stability by TOF-MS. Additionally to the singly charged parent ion and its fragments, we find a significant yield of doubly and triply charged parent ions not previously reported. In order to characterize these species, we used high-level theoretical methods based on density functional theory (DFT), coupled-cluster (CC), and generalized valence bond (GVB) to calculate the structure, relative stabilities, and bonding of these dications and trications. The most stable isomers exhibit unusual carbon-chlorine multiple bonding: a terminal C═Cl double bond in a formyl-like CHCl moiety (1, rC-Cl = 1.621 Å) and a ketene-like C═C═Cl cumulated species (2, rC-Cl = 1.542 Å). The calculations suggest that an excited state of 2 has a nitrile-like C≡Cl triple bond structure.

Are One-Electron Bonds Any Different from Standard Two-Electron Covalent Bonds?

The nature of the chemical bond is perhaps the central subject in theoretical chemistry. Our understanding of the behavior of molecules developed amazingly in the last century, mostly with the rise of quantum mechanics (QM) and QM-based theories such as valence bond theory and molecular orbital theory. Such theories are very successful in describing molecular properties, but they are not able to explain the origin of the chemical bond. This problem was first analyzed in the 1960s by Ruedenberg, who showed that covalent bonds are the direct result of quantum interference between one-electron states. The generality of this result and its quantification were made possible through the recent development of the generalized product function energy partitioning (GPF-EP) method by our group, which allows the partitioning of the electronic density and energy into their interference and quasi-classical (noninterference) contributions. Furthermore, with GPF wave functions these effects can be analyzed separately for each bond of a molecule. This interference energy analysis has been applied to a large variety of molecules, including diatomics and polyatomics, molecules with single, double, and triple bonds, molecules with different degrees of polarity, linear or branched molecules, cyclic or acyclic molecules, conjugated molecules, and aromatics, in order to verify the role played by quantum interference. In all cases the conclusion is exactly the same: for each bond in each of the molecules considered, the main contribution to its stability comes from the interference term. Two-center one-electron (2c1e) bonds are the simplest kind of chemical bonds. Yet they are often viewed as odd or unconventional cases of bonding. Are they any different from conventional (2c2e) bonds? If so, what differences can we expect in the nature of (2c1e) bonds relative to electron-pair bonds? In this Account, we extend the GPF-EP method to describe bonds involving N electrons in M orbitals (N < M) and show its application to (2c1e) bonds. As examples we chose the molecules H2+, H3C·CH3+, B2H4-, [Cu·BH3(PH3)3], and an alkali-metal cation dimer, and we evaluated the components of the electronic energy and density, which account for the formation of the bond, and compared the results with those for the respective analogous molecules exhibiting the "conventional" two-electron bond. In all cases, it was verified that interference is the dominant effect for the one-electron bonds. The GPF-EP results clearly indicate that molecules exhibiting (2c1e) bonds should not be considered as special systems, since one- and two-electron bonds result from quantum interference and therefore there is no conceptual difference between them. Moreover, these results show that quantum interference provides a way to unify the chemical bond concept.

Diboryne Nanostructures Stabilized by Multitopic N-Heterocyclic Carbenes: A Computational Study.

Different families of nanomaterials produced from the stabilization of diboryne (B≡B) units by multitopic N-heterocyclic carbenes (NHCs), such as nanowires, nanorings, and nanotents, were studied by computational methods. Density functional theory calculations with and without periodic boundary conditions were applied to estimate the dependence of the electronic and thermochemical properties of different diboryne macromolecules with respect to the nature of the bridging ligand. Our results show that all diboryne nanostructures studied herein are viable candidates for synthesis. The Janus-type multitopic naphthobis(imidazolylidene) (5), anthrobis(imidazolylidene) (10), and pyracenetetrakis(imidazolylidene) (16) compounds are the best candidates for generating diboryne nanowires. A path to covalent organic frameworks, nanocages, and nanotubes from the optimized diboryne nanostructures is also described. Rather than just scientific curiosity, diboryne nanostructures emerge as interesting targets for the synthesis of main-group nanomaterials.

Quantum Interference Contribution to the Dipole Moment of Diatomic Molecules.

The interference energy partitioning analysis method developed by our group and used to study the nature of the chemical bond was extended to partition the electric dipole moment in quasi-classical and interference contributions. Our results show that interference participates in charge displacement in polar molecules, providing, directly or indirectly, a relevant contribution for the total dipole moment. A linear correlation was found between the interference contribution of the dipole moment from the bond electron group, µINT(bond), and the difference of electronegativity of the atoms which form the bond, ΔXAB. This interesting result reinforces the fact that electronegativity is not a property of an atom alone, but rather a property of the atom in the molecule and that ΔXAB can only be associated with that part of the total charge displacement resulting from the formation of the chemical bond. The partitioning of the total dipole moment into quasi-classical and interference contributions provides new insights about the reasons for the failure of the ΔXAB criterion in predicting the correct orientation of the dipole moment in several molecules. The results of the present work also bring additional evidence for the previously proposed mechanism of formation of polar bonds.

Time-dependent density functional theory description of total photoabsorption cross sections.

The time-dependent version of the density functional theory (TDDFT) has been used to calculate the total photoabsorption cross section of a number of molecules, namely, benzene, pyridine, furan, pyrrole, thiophene, phenol, naphthalene, and anthracene. The discrete electronic pseudo-spectra, obtained in a L2 basis set calculation were used in an analytic continuation procedure to obtain the photoabsorption cross sections. The ammonia molecule was chosen as a model system to compare the results obtained with TDDFT to those obtained with the linear response coupled cluster approach in order to make a link with our previous work and establish benchmarks.

Nature of the chemical bond and origin of the inverted dipole moment in boron fluoride: a generalized valence bond approach.

The generalized product function energy partitioning (GPF-EP) method has been applied to investigate the nature of the chemical bond and the origin of the inverted dipole moment of the BF molecule. The calculations were carried out with GPF wave functions treating all of the core electrons as a single Hartree-Fock group and the valence electrons at the generalized valence bond perfect-pairing (GVB-PP) or full GVB levels, with the cc-pVTZ basis set. The results show that the chemical structure of both X (1)Σ(+) and a (3)Π states is composed of a single bond. The lower dissociation energy of the excited state is attributed to a stabilizing intraatomic singlet coupling involving the B 2sp-like lobe orbitals after bond dissociation. An increase of electron density on the B atom caused by the reorientation of the boron 2sp-like lobe orbitals is identified as the main responsible effect for the electric dipole inversion in the ground state of BF. Finally, it is shown that π back-bonding from fluorine to boron plays a minor role in the electron density displacement to the bonding region in both states. Moreover, this effect is associated with changes in the quasi-classical component of the electron density only and does not contribute to covalency in either of the states. Therefore, at least for the case of the BF molecule, the term back-bonding is misleading, since it does not contribute to the bond formation.

Coupled-Cluster Study of the Lower Energy Region of the Ground Electronic State of the HSO2 Potential Energy Surface.

This work reports CCSD(T)/aug-cc-pV(T+d)Z ab initio calculations for the lower energy region of the ground electronic state of the HSO2 system. Optimized geometries, total energies, zero-point vibrational energies, frequencies, complete basis set extrapolations, and reaction paths are reported at the same level of calculation. The connection of the two minima (synperiplanar HOSO and HSO2) with the dissociation limit H + SO2 through the van der Waals minimum H···SO2 was established. An important quantitative discrepancy with previous works is the fact that at the present level of calculation the energy difference between transition states connecting the global minimum synperiplanar HOSO to the HSO2 minimum (TS5) and to the van der Waals minimum H···SO2 (TS6) is negligible, implying that the forward barriers after the synperiplanar HOSO global minimum have practically the same height. This result suggests that these two transition states may be involved in the path of the global minimum toward the exit channel H + SO2. As a consequence, trajectories for the OH + SO collisions could evolve through the well formed by the HSO2 minimum, therefore opening two competitive channels for the OH + SO → H + SO2 reaction, a fact never reported in trajectory calculations.

The non-covalent nature of the molecular structure of the benzene molecule.

The benzene molecule is one of the most emblematic systems in chemistry, with its structural features being present in numerous different compounds. We have carried out an analysis of the influence of quantum mechanical interference on the geometric features of the benzene molecule, showing that many of the characteristics of its equilibrium geometry are a consequence of non-covalent contributions to the energy. This result implies that quasi-classical reasoning should be sufficient to predict the defining aspects of the benzene structure such as its planarity and equivalence of its bond lengths.

Interference energy in C-H and C-C bonds of saturated hydrocarbons: dependence on the type of chain and relationship to bond dissociation energy.

Interference energy for C-H and C-C bonds of a set of saturated hydrocarbons is calculated by the generalized product function energy partitioning (GPF-EP) method in order to investigate its sensitivity to the type of chain and also its contribution to the bond dissociation energy. All GPF groups corresponding to chemical bonds are calculated by use of GVB-PP wave functions to ensure the correct description of bond dissociation. The results show that the interference energies are practically the same for all the C-H bonds, presenting only small variations (0.5 kcal.mol(-1)) due to the structural changes in going from linear to branched and cyclic chains. A similar trend is verified for the C-C bonds, the sole exception being the cyclopropane molecule, for which only the C-C bond exhibits a more significant variation. On the other hand, although the interference energy is quantitatively the most important contribution to the bond dissociation energy (DE), one cannot predict DE only from the bond interference energy. Differences in the dissociation energies of C-C and C-H bonds due to structural changes in the saturated hydrocarbons can be mainly attributed to quasi-classical effects.