Search

An analysis using the formalism of crystalline orbitals for extended systems with periodicity in one dimension demonstrates that any antiferromagnetic and half-metallic spin-polarization of the edge states in n-acenes, and more generally in zigzag graphene nanoislands and nanoribbons of finite width, would imply a spin contamination that increases proportionally to system size, in sharp and clear contradiction with the implications of Lieb's theorem for compensated bipartite lattices and the expected value for a singlet (S = 0) electronic ground state. Verifications on naphthalene, larger n-acenes (n = 3-10) and rectangular nanographene islands of increasing size, as well as a comparison using unrestricted Hartree-Fock theory along with basis sets of improving quality against various many-body treatments demonstrate altogether that antiferromagnetism and half-metallicity in extended graphene nanoribbons will be quenched by an exact treatment of electron correlation, at the confines of non-relativistic many-body quantum mechanics. Indeed, for singlet states, symmetry-breakings in spin-densities are necessarily the outcome of a too approximate treatment of static and dynamic electron correlation in single-determinantal approaches, such as unrestricted Hartree-Fock or Density Functional Theory. In this context, such as the size-extensive spin-contamination to which it relates, half-metallicity is thus nothing else than a methodological artefact. [ABSTRACT FROM AUTHOR]

A benchmark theoretical study of the electronic ground state and of the vertical and adiabatic singlet-triplet (ST) excitation energies of benzene (n=1) and n-acenes (C4n+2H2n+4) ranging from naphthalene (n=2) to heptacene (n=7) is presented, on the ground of single- and multireference calculations based on restricted or unrestricted zero-order wave functions. High-level and large scale treatments of electronic correlation in the ground state are found to be necessary for compensating giant but unphysical symmetry-breaking effects in unrestricted single-reference treatments. The composition of multiconfigurational wave functions, the topologies of natural orbitals in symmetry-unrestricted CASSCF calculations, the T1 diagnostics of coupled cluster theory, and further energy-based criteria demonstrate that all investigated systems exhibit a 1Ag singlet closed-shell electronic ground state. Singlet-triplet (S0-T1) energy gaps can therefore be very accurately determined by applying the principles of a focal point analysis onto the results of a series of single-point and symmetry-restricted calculations employing correlation consistent cc-pVXZ basis sets (X=D, T, Q, 5) and single-reference methods [HF, MP2, MP3, MP4SDQ, CCSD, CCSD(T)] of improving quality. According to our best estimates, which amount to a dual extrapolation of energy differences to the level of coupled cluster theory including single, double, and perturbative estimates of connected triple excitations [CCSD(T)] in the limit of an asymptotically complete basis set (cc-pV∞Z), the S0-T1 vertical excitation energies of benzene (n=1) and n-acenes (n=2–7) amount to 100.79, 76.28, 56.97, 40.69, 31.51, 22.96, and 18.16 kcal/mol, respectively. Values of 87.02, 62.87, 46.22, 32.23, 24.19, 16.79, and 12.56 kcal/mol are correspondingly obtained at the CCSD(T)/cc-pV∞Z level for the S0-T1 adiabatic excitation energies, upon including B3LYP/cc-PVTZ corrections for zero-point vibrational energies. In line with the absence of Peierls distortions, extrapolations of results indicate a vanishingly small S0-T1 energy gap of 0 to ∼4 kcal/mol (∼0.17 eV) in the limit of an infinitely large polyacene. [ABSTRACT FROM AUTHOR]

A benchmark theoretical determination of the electron affinities of benzene and linear oligoacenes ranging from naphthalene to hexacene is presented, using the principles of a focal point analysis. These energy differences have been obtained from a series of single-point calculations at the Hartree–Fock, second-, third-, and partial fourth-order Mo\ller–Plesset (MP2, MP3, and MP4SDQ) levels and from coupled cluster calculations including single and double excitations (CCSD) as well as perturbative estimates of connected triple excitations [CCSD(T)], using basis sets of improving quality, containing up to 1386, 1350, 1824, 1992, 1630, and 1910 basis functions in the computations, respectively. Studies of the convergence properties of these energy differences as a function of the size of the basis set and order attained in electronic correlation enable a determination of the vertical electron affinities of the four larger terms of the oligoacene (C2+4nH2+2n) series within chemical accuracy (0.04 eV). According to our best estimates, these amount to +0.28, +0.82, +1.21, and +1.47 eV when n=3, 4, 5, and 6. Adiabatic electron affinities have been further calculated by incorporating corrections for zero-point vibrational energies and for geometrical relaxations. The same procedure was applied to determine the vertical electron affinities of benzene and naphthalene, which are found to be markedly negative (∼-1.53 and ∼-0.48 eV, respectively). Highly quantitative insights into experiments employing electron transmission spectroscopy on these compounds were also amenable from such an approach, provided diffuse atomic functions are deliberately removed from the basis set, in order to enforce confinement in the molecular region and enable a determination of pseudoadiabatic electron affinities (with respect to the timescale of nuclear motions). Comparison was made with calculations employing density functional theory and especially designed models that exploit the integer discontinuity in the potential or incorporate a potential wall in the unrestricted Kohn–Sham orbital equation for the anion. [ABSTRACT FROM AUTHOR]

A comprehensive study, throughout the valence region, of the electronic structure and electron momentum density distributions of the four conformational isomers of n-pentane is presented. Theoretical (e,2e) valence ionization spectra at high electron impact energies (1200 eV+electron binding energy) and at azimuthal angles ranging from 0° to 10° in a noncoplanar symmetric kinematical setup are generated according to the results of large scale one-particle Green’s function calculations of Dyson orbitals and related electron binding energies, using the third-order algebraic-diagrammatic construction [ADC(3)] scheme. The results of a focal point analysis (FPA) of relative conformer energies [A. Salam and M. S. Deleuze, J. Chem. Phys. 116, 1296 (2002)] and improved thermodynamical calculations accounting for hindered rotations are also employed in order to quantitatively evaluate the abundance of each conformer in the gas phase at room temperature and reliably predict the outcome of experiments on n-pentane employing high resolution electron momentum spectroscopy. Comparison with available photoelectron measurements confirms the suggestion that, due to entropy effects, the trans-gauche (tg) conformer strongly dominates the conformational mixture characterizing n-pentane at room temperature. Our simulations demonstrate therefore that experimental measurements of (e,2e) valence ionization spectra and electron momentum distributions would very consistently and straightforwardly image the topological changes and energy variations that molecular orbitals undergo due to torsion of the carbon backbone. The strongest fingerprints for the most stable conformer (tt) are found for the electron momentum distributions associated with ionization channels at the top of the inner-valence region, which sensitively image the development of methylenic hyperconjugation in all-staggered n-alkane chains. [ABSTRACT FROM AUTHOR]

In this work we discuss in some computational and analytical details the issue of half-metallicity in zig-zag graphene nanoribbons and nanoislands of finite width, i.e. the coexistence of metallic nature for electrons with one spin orientation and insulating nature for the electrons of opposite spin, which has been recently predicted from so-called first-principle calculations employing Density Functional Theory. It is mathematically demonstrated and computationally verified that, within the framework of non-relativistic and time-independent quantum mechanics, like the size-extensive spin-contamination to which it relates, half-metallicity is nothing else than a methodological artefact, due to a too approximate treatment of electron correlation in the electronic ground state.

Due to its unique mechanical properties, graphene can be applied for reinforcement in nanocomposites. We analyse the Young's modulus of graphene at the semi-empirical PM6 level of theory. The internal forces are calculated and the Young's modulus is predicted for a finite graphene sheet when external strain is applied on the system. These results are in a good agreement with theoretical and experimental results from the literature giving values of about 1 TPa for the Young's modulus. Stress-strain curves are computed for elongation up to 20%. In addition, the influence of the presence of a single vacancy, as well as for oxygenation of a vacancy, on the mechanical properties of graphene has been analysed. Our results indicate that when applying the deformation locally onto the system, higher local stress can be induced, as confirmed by Finite Element Analysis. Also, the presence of structural defects in the system will stiffen the system upon low strain, but reduces the elastic limit from more than 20% strain for pristine graphene to less than 10% strain when defects are present.

A benchmark theoretical study of the electronic ground state and of the vertical and adiabatic singlet-triplet (ST) excitation energies of n-acenes (C(4n+2)H(2n+4)) ranging from octacene (n = 8) to undecacene (n = 11) is presented. The T1 diagnostics of coupled cluster theory and further energy-based criteria demonstrate that all investigated systems exhibit predominantly a (1)A(g) singlet closed-shell electronic ground state. Singlet-triplet (S(0)-T(1)) energy gaps can therefore be very accurately determined by applying the principle of a focal point analysis (FPA) onto the results of a series of single-point and symmetry-restricted calculations employing correlation consistent cc-pVXZ basis sets (X = D, T, Q, 5) and single-reference methods [HF, MP2, MP3, MP4SDQ, CCSD, and CCSD(T)] of improving quality. According to our best estimates, which amount to a dual extrapolation of energy differences to the level of coupled cluster theory including single, double, and perturbative estimates of connected triple excitations [CCSD(T)] in the limit of an asymptotically complete basis set (cc-pV∞Z), the S(0)-T(1) vertical (adiabatic) excitation energies of these compounds amount to 13.40 (8.21), 10.72 (6.05), 8.05 (3.67), and 7.10 (2.58) kcal/mol, respectively. In line with the absence of Peierls distortions (bond length alternations), extrapolations of results obtained at this level for benzene (n = 1) and all studied n-acenes so far (n = 2-11) indicate a vanishing S(0)-T(1) energy gap, in the limit of an infinitely large polyacene, within an uncertainty of 1.5 kcal/mol (0.06 eV). Lacking experimental values for the S(0)-T(1) energy gaps of n-acenes larger than hexacene, comparison is made with recent optical and electrochemical determinations of the HOMO-LUMO band gap. Further issues such as scalar relativistic, core correlation, and diagonal Born-Oppenheimer corrections (DBOCs) are tentatively examined., (© 2011 American Chemical Society)

We have used a potential wall method to investigate the role of d orbitals in the a(2) singly-occupied molecular orbitals of (2)A(2) negative ion states of two molecular series: pyridine, phosphabenzene, arsabenzene, stibabenzene (C(5)H(5)X, X = {N, P, As, Sb}), and furan, thiophene, selenophene, tellurophene (C(4)H(4)X, X = {O, S, Se, Te}). Unlike for the lower lying doubly occupied orbitals, heteroatom d-carbon p in-phase (bonding) interactions in these a(2) orbitals are clearly identified and explain the 0.5 eV stabilization of the (2)A(2) radical anion state in those compounds where the heteroatoms have d orbitals in the valence shell, compared to compounds where d orbitals are missing in the valence shell of the heteroatoms. The performance of both the potential wall approach and the approximate expression of Tozer and De Proft for calculating negative electron affinities has been also investigated, through a comparison with results obtained using electron-transmission spectroscopy experiments.

An extensive study, throughout the valence region, of the electronic structure, ionization spectrum, and electron momentum distributions of norbornadiene is presented, on the ground of accurate calculations of valence one-electron and shake-up ionization energies and of the related Dyson orbitals, using one-particle Green's function (1p-GF) theory in conjunction with the so-called third-order algebraic diagrammatic construction scheme [ADC(3)]. Comparison is made with results obtained from standard (B3LYP) Kohn-Sham orbitals and measurements employing electron momentum spectroscopy, taking into account the contamination of inner- and outer-valence spectral bands by numerous shake-up states. Four relatively intense shake-up lines at 12.1, 16.4, 17.6, and 17.8 eV are found to yield recognizable spectral fingerprints in the EMS experiments. Valence bands at electron binding energies larger than 20 eV are subject to a complete breakdown of the orbital picture of ionization.

The [3 + 2] and [3 + 3] cyclodimerisation processes of small nitrile oxides, XCNO (X = F, Cl, Br, CN, CH(3)) are investigated by ab initio coupled cluster theory at the CCSD, CCSD(T) and MR-AQCC levels for the first time. The favoured dimerisation process is a multi-step reaction to furoxans (1,2,5-oxadiazole-2-oxides) involving dinitrosoalkene-like intermediates with diradical character. The rate determining step for all but the F-species is the first, corresponding to the C-C bond formation. The kinetic energy barrier depends on the nature of the substituent X, generally increasing with decreasing electronegativity and increasing pi-donor ability of the substituent: F (DeltaG(298) = 0 kJ mol(-1)) < Cl (72) < Br (90) < CH(3) (104) < CN (114) (MR-AQCC(2,2)//UB3LYP/cc-pVTZ). Following initial C-C bond formation, three possible dinitrosoethylene diradical pathways are explored. Two of them are new, and one of them is a low-energy three-step path with implications for cycloreversion, tautomerism and detection of dinitrosoethylene intermediates. Alternative one-step, concerted [3 + 2] and [3 + 3] cyclodimerisation processes leading to 1,2,4-oxadiazole-4-oxides and 1,4,2,5-dioxadiazines have kinetic energy barriers around 100-240 kJ mol(-1) (CCSD//B3LYP), some 1.6 to 2.5 times higher than those leading to furoxans, supporting the experimental observations of furoxan formation as nitrile oxide loss channels during storage, trapping/re-vaporisation and reactions of nitrile oxides. Potential polymerisation initiation processes for NCCNO, involving the 1,2-dipolar NC substituent are also explored.

Calculations of electron momentum distributions for equilibrium geometries, employing advanced Dyson orbital theories and statistical thermodynamics beyond the RRHO approximation, fail to quantitatively reproduce the outermost momentum distribution profile inferred from (e, 2e) electron impact ionization experiments on ethanol employing high-resolution electron momentum spectroscopy. A very detailed study of the influence on this momentum distribution of nuclear dynamics in the initial ground state and in the final ionized state is presented according to a thermal averaging over exceedingly large sets of model structures as well as Born-Oppenheimer molecular dynamical simulations on the potential energy surface of the radical cation. Our results give very convincing albeit qualitative indications that the strong turn-up of the (e, 2e) ionization intensities characterizing the highest occupied molecular orbital (HOMO) of ethanol at low electron momenta is the combined result of (1) the extraordinarily flat nature of the conformational energy map of ethanol, which enables significant departures from energy minima in the ground electronic state, (2) strong anomeric interactions between an oxygen lone pair and the central C-C bond for the minor but significant fraction of conformers exhibiting a hydroxyl torsion angle (alpha) at around 90 degrees, and, last but not least, (3) the possibility to observe with this minor conformer fraction ultrafast and highly significant extensions of the central C-C bond, resulting, in turn, in an enhanced delocalization of the HOMO from the oxygen lone pair region onto the methyl group, immediately after the sudden removal of an electron. This charge transfer appears to occur at the very first stages, that is, within an effective time scale on the order of approximately 10 fs, of an ultrafast dissociation of the ethanol radical cation into a methyl radical and a protonated form of formaldehyde.

The parent furoxan (1,2,5-oxadiazole 2-oxide), synthesized from glyoxime and NO(2)(g), has been investigated in the gas phase for the first time by mid-infrared and He I photoelectron spectroscopy, and in the liquid phase by Raman spectroscopy. The ground-state geometry has been obtained from quantum-chemical calculations at the B3LYP, MPn (n = 2-4), CISD, QCISD, CCSD, CCSD(T), RSPTn (n = 2,3), MRCI, and MR-AQCC levels using 6-311++G(2d,2p), cc-pVTZ, aug-cc-pVTZ, cc-pCVTZ, and cc-pVQZ basis sets. Furoxan is predicted to be planar, with a strong exocyclic and a relatively weak endocyclic N-O bond. The furoxan moiety is electron rich, indicated e.g. by a large negative NPA charge (-0.46 e). According to various aromaticity indices, furoxan is nearly as aromatic as furan and furazan. Unlike alkyl- and cyano-substituted furoxans, the parent furoxan, upon thermolysis, does not cleave to the monomer nitrile oxide, yielding only HNCO, HCN, CO(2), CO, NO, and H(2)O decomposition products.

An extensive study, throughout the valence region, of the electronic structure, ionization spectrum, and electron momentum distributions of ethanol is presented, on the ground of a model that focuses on a mixture of the gauche and anti conformers in their energy minimum form, using weight coefficients obtained from thermostatistical calculations that account for the influence of hindered rotations. The analysis is based on accurate calculations of valence one-electron and shakeup ionization energies and of the related Dyson orbitals, using one-particle Green's Function (1p-GF) theory in conjunction with the so-called third-order Algebraic Diagrammatic Construction scheme [ADC(3)]. The confrontation against available UPS (HeI) measurements indicates the presence in the spectral bands of significant conformational fingerprints at outer-valence ionization energies ranging from approximately 14 to approximately 18 eV. The shakeup onset is located at approximately 24 eV, and a shoulder at approximately 14.5 eV in the He I spectrum can be specifically ascribed to the minor anti (C(s)) conformer fraction. Thermally and spherically averaged Dyson orbital momentum distributions are computed for seven resolvable bands in model (e, 2e) ionization spectra at an electron impact energy of 1.2 keV. A comparison is made with results obtained from standard (B3LYP) Kohn-Sham orbitals and EMS measurements employing a high-resolution spectrometer of the third generation. The analysis is qualitatively in line with experiment and reveals a tremendously strong influence of the molecular conformation on the outermost electron momentum distributions. Quantitatively significant discrepancies with experiment can nonetheless be tentatively ascribed to strong dynamical disorder in the gas phase molecular structure.

The potential energy surface of norbornane in its dicationic singlet ground state has been investigated in detail using density functional theory along with the nonlocal hybrid and gradient-corrected Becke three-parameter Lee-Yang-Parr functional (B3LYP) and the cc-pVDZ basis set. For the sake of more quantitative insight into the chemical reactions induced by double ionization of norbornane, this study was supplemented by a calculation of basic thermodynamic state functions coupled to a focal point analysis of energy differences obtained using correlation treatments and basis sets of improving quality, enabling an extrapolation of these energy differences at the CCSD(T) level in the limit of an asymptotically complete (cc-pV infinity Z) basis set. Our results demonstrate the likelihood of an ultrafast intramolecular rearrangement of the saturated hydrocarbon cage after a sudden removal of two electrons into a kinetically metastable five-membered cyclic C5H8+-CH+-CH3 intermediate, prior to a Coulomb explosion into C5H7+=CH2 and CH3+ fragments, which might explain a tremendous rise of electron-impact (e, 2e) ionization cross sections at electron binding energies around the double-ionization threshold. The first step is straightforward and strongly exothermic (DeltaH298 = -114.0 kcal mol-1). The second step is also exothermic (DeltaH298 = -10.2 kcal mol-1) but requires an activation enthalpy (DeltaH298) of 39.7 kcal/mol. The various factors governing the structure of this intermediate, such as electrostatic interactions, inductive effects, cyclic strains, and methylenic hyperconjugation interactions, are discussed in detail.

We present a systematic theoretical study based on semi-empirical, Hartree-Fock (HF), and density functional theory (DFT) models of a series of polycyclic aromatic hydrocarbons (PAHs) that exhibit hollow sites. In this study we focus particularly on the magnetic criteria of aromaticity, namely (1)H NMR and nucleus-independent chemical shifts (NICS), and on their relationships with other electronic properties. The computed shifts and NICS indices indicate that an external magnetic field induces exceptionally strong ring currents in even-layered PAH doughnuts, in particular in the layer directly adjacent to the central hole of double-layered compounds. These exceptionally strong ring currents also correlate with particularly small HOMO-LUMO gaps and electronic excitation energies and to abnormally high polarizabilities, indicating in turn that these compounds have a more pronounced metallic character. Comparison is made with further depictions of aromaticity in these systems and in [18]-[66]annulene rings by employing topological, structural, and energetic criteria.

The mechanism of 1,2-addition reactions of HF and HCl to Si=Si, Si=C, and C=C bonds has been investigated by ab initio quantum chemical methods. Geometries and relative energies of the stationary points and all the transition states were determined by using the MP2/6-311++G(d,p), B3LYP/6-311++G(d,p), and CBS-Q levels of theory. The investigated reactions can be characterized by two main thermodynamic profiles. The type in which the reagent molecule attacks a carbon atom is moderately exothermic with a high activation barrier. The second type in which a hydrogen halide attacks a silicon is strongly exothermic with a low activation energy. At the early stage of all the reactions a weakly bonded initial complex is found which indicates that the initial step of all the reactions is an electrophilic attack of hydrogen halide. The geometry and charge distribution of the transition state of the reactions indicate two main types of mechanism. If silicon is attacked, the halogen-silicon bond formation precedes the H-Y bond breaking. If, however, carbon is attacked, the first step is always an ionic dissociation of the hydrogen halide and a carbenium ion formation, which is stabilized by the C-Y bond formation in the final step of the reaction. The reaction diagrams and proposed mechanisms explain the experimentally found regioselectivity well.

The mechanism of 1,2-addition reactions of water, methanol, and trifluoromethanol to Si=Si, Si=C, and C=C bonds has been investigated by ab initio quantum chemical methods. Geometries and relative energies of the stationary points and all the transition states were determined using the MP2/6-311++G(d,p), B3LYP/6-311++G(d,p), and CBS-Q levels of theory. The investigated reactions can be characterized by two main thermodynamical profiles. The type in which the reagent molecule attacks a carbon atom is moderately exothermic, with a high activation barrier. The second type, in which water or alcohol attacks a silicon, is strongly exothermic, with a small activation energy. At the early stage of all the reactions, a weakly bonded initial complex is found which determines the further mechanism of the reaction. On the basis of the HOMO, LUMO, and Laplacian of electron distribution of disilene and silene, several mechanisms have been assumed, depending on the substrate (disilene, substituted disilene, silene, or ethene) and the reagent (water, methanol, or trifluoromethanol). The reaction diagrams and proposed mechanisms explain the experimentally found regioselectivity and diastereoselectivity well.