Theoretical study of the oxidation mechanisms of thiophene initiated by hydroxyl radicals.

The mechanisms for the oxidation of thiophene by OH radicals under inert conditions (Ar) have been studied using density functional theory in conjunction with various exchange-correlation functionals. These results were compared with benchmark CBS-QB3 theoretical results. Kinetic rate constants were estimated by means of variational transition state theory (VTST) and the statistical Rice-Ramsperger-Kassel-Marcus (RRKM) theory. Effective rate constants were calculated via a steady-state analysis based upon a two-step model reaction mechanism. In line with experimental results, the computed branching ratios indicate that the most kinetically efficient process involves OH addition to a carbon atom adjacent to the sulfur atom. Due to the presence of negative activation energies, pressures larger than 10(4) bar are required to reach the high-pressure limit. Nucleus-independent chemical shift indices and natural bond orbital analysis show that the computed activation energies are dictated by changes in aromaticity and charge-transfer effects due to the delocalization of lone pairs from sulfur to empty π(*) orbitals. Graphical Abstract CBS-QB3 energy profiles for the reaction pathways 1-3 characterizing the oxidation of thiophene by hydroxyl radicals into the related products.

Theoretical study of molecular vibrations in electron momentum spectroscopy experiments on furan: an analytical versus a molecular dynamical approach.

The influence of thermally induced nuclear dynamics (molecular vibrations) in the initial electronic ground state on the valence orbital momentum profiles of furan has been theoretically investigated using two different approaches. The first of these approaches employs the principles of Born-Oppenheimer molecular dynamics, whereas the so-called harmonic analytical quantum mechanical approach resorts to an analytical decomposition of contributions arising from quantized harmonic vibrational eigenstates. In spite of their intrinsic differences, the two approaches enable consistent insights into the electron momentum distributions inferred from new measurements employing electron momentum spectroscopy and an electron impact energy of 1.2 keV. Both approaches point out in particular an appreciable influence of a few specific molecular vibrations of A1 symmetry on the 9a1 momentum profile, which can be unravelled from considerations on the symmetry characteristics of orbitals and their energy spacing.

Theoretical study of the oxidation mechanisms of naphthalene initiated by hydroxyl radicals: the OH-addition pathway.

The oxidation mechanisms of naphthalene by OH radicals under inert (He) conditions have been studied using density functional theory along with various exchange-correlation functionals. Comparison has been made with benchmark CBS-QB3 theoretical results. Kinetic rate constants were correspondingly estimated by means of transition state theory and statistical Rice-Ramsperger-Kassel-Marcus (RRKM) theory. Comparison with experiment confirms that, on the OH-addition reaction pathway leading to 1-naphthol, the first bimolecular reaction step has an effective negative activation energy around -1.5 kcal mol(-1), whereas this step is characterized by an activation energy around 1 kcal mol(-1) on the OH-addition reaction pathway leading to 2-naphthol. Effective rate constants have been calculated according to a steady state analysis upon a two-step model reaction mechanism. In line with experiment, the correspondingly obtained branching ratios indicate that, at temperatures lower than 410 K, the most abundant product resulting from the oxidation of naphthalene by OH radicals must be 1-naphthol. The regioselectivity of the OH(•)-addition onto naphthalene decreases with increasing temperatures and decreasing pressures. Because of slightly positive or even negative activation energies, the RRKM calculations demonstrate that the transition state approximation breaks down at ambient pressure (1 bar) for the first bimolecular reaction steps. Overwhelmingly high pressures, larger than 10(5) bar, would be required for restoring to some extent (within ∼5% accuracy) the validity of this approximation for all the reaction channels that are involved in the OH-addition pathway. Analysis of the computed structures, bond orders, and free energy profiles demonstrate that all reaction steps involved in the oxidation of naphthalene by OH radicals satisfy Leffler-Hammond's principle. Nucleus independent chemical shift indices and natural bond orbital analysis also show that the computed activation and reaction energies are largely dictated by alterations of aromaticity, and, to a lesser extent, by anomeric and hyperconjugative effects.

Theoretical study of the oxidation mechanisms of naphthalene initiated by hydroxyl radicals: the H abstraction pathway.

Reaction mechanisms for the initial stages of naphthalene oxidation at high temperatures (T ≥ 600 K) have been studied theoretically using density functional theory along with various exchange-correlation functionals, as well as the benchmark CBS-QB3 quantum chemical approach. These stages correspond to the removal of hydrogen atoms by hydroxyl radical and the formation thereby of 1- and 2-naphthyl radicals. Bimolecular kinetic rate constants were estimated by means of transition state theory. The excellent agreement with the available experimental kinetic rate constants demonstrates that a two-step reaction scheme prevails. Comparison with results obtained with density functional theory in conjunction with various exchange-correlation functionals also shows that DFT remains unsuited for quantitative insights into kinetic rate constants. Analysis of the computed structures, bond orders, and free energy profiles demonstrates that the reaction steps involved in the removal of hydrogen atoms by OH radicals satisfy Hammond's principle. Computations of branching ratios also show that these reactions do not exhibit a particularly pronounced site-selectivity.

Electron momentum spectroscopy of 1-butene: a theoretical analysis using molecular dynamics and molecular quantum similarity.

The results of experimental studies of the valence electronic structure of 1-butene by means of electron momentum spectroscopy (EMS) have been reinterpreted on the basis of molecular dynamical simulations in conjunction with the classical MM3 force field. The computed atomic trajectories demonstrate the importance of thermally induced nuclear dynamics in the electronic neutral ground state, in the form of significant deviations from stationary points on the potential energy surface and considerable variations of the C-C-C-C dihedral angle. These motions are found to have a considerable influence on the computed spectral bands and outer-valence electron momentum distributions. Euclidean distances between spherically averaged electron momentum densities confirm that thermally induced nuclear motions need to be fully taken into account for a consistent interpretation of the results of EMS experiments on conformationally flexible molecules.

Photoelectron and electron momentum spectroscopy of tetrahydrofuran from a molecular dynamical perspective.

The results of experimental studies of the valence electronic structure of tetrahydrofuran employing He I photoelectron spectroscopy as well as Electron Momentum Spectroscopy (EMS) have been reinterpreted on the basis of Molecular Dynamical simulations employing the classical MM3 force field and large-scale quantum mechanical simulations employing Born-Oppenheimer Molecular Dynamics in conjunction with the dispersion corrected ωB97XD exchange-correlation functional. Analysis of the produced atomic trajectories demonstrates the importance of thermal deviations from the lowest energy path for pseudorotation, in the form of considerable variations of the ring-puckering amplitude. These deviations are found to have a significant influence on several outer-valence electron momentum distributions, as well as on the He I photoelectron spectrum.

Half-metallicity of graphene nanoribbons and related systems: a new quantum mechanical El Dorado for nanotechnologies... or a hype for materials scientists?

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.

Focal point analysis of the singlet-triplet energy gap of octacene and larger acenes.

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.

Bonding in negative ions: the role of d orbitals in the heavy analogues of pyridine and furan radical anions.

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.

Electron momentum spectroscopy of norbornadiene at the benchmark ADC(3) level.

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.

Probing nuclear dynamics in momentum space: a new interpretation of (e, 2e) electron impact ionization experiments on ethanol.

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.

Benchmark Dyson orbital study of the ionization spectrum and electron momentum distributions of ethanol in conformational equilibrium.

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.

Study of the photoelectron and electron momentum spectra of cyclopentene using benchmark Dyson orbital theories.

A complete study of the valence electronic structure and related electronic excitation properties of cyclopentene in its C(s) ground state geometry is presented. Ionization spectra obtained from this compound by means of photoelectron spectroscopy (He I and He II) and electron momentum spectroscopy have been analyzed in details up to electron binding energies of 30 eV using one-particle Green's function (1p-GF) theory along with the outer-valence (OVGF) and the third-order algebraic diagrammatic construction [ADC(3)] schemes. The employed geometries derive from DFT/B3LYP calculations in conjunction with the aug-cc-pVTZ basis set, and closely approach the structures inferred from experiments employing microwave spectroscopy or electron diffraction in the gas phase. The 1p-GF/ADC(3) calculations indicate that the orbital picture of ionization breaks down at electron binding energies larger than approximately 17 eV in the inner-valence region, and that the outer-valence 7a' orbital is also subject to a significant dispersion of the ionization intensity over shake-up states. This study confirms further the rule that OVGF pole strengths smaller than 0.85 foretell a breakdown of the orbital picture of ionization at the ADC(3) level. Spherically averaged (e, 2e) electron momentum distributions at an electron impact energy of 1200 eV that were experimentally inferred from an angular analysis of EMS intensities have been interpreted by comparison with accurate simulations employing ADC(3) Dyson orbitals. Very significant discrepancies were observed with momentum distributions obtained from several outer-valence ionization bands using standard Kohn-Sham orbitals.

Theoretical study of the fragmentation pathways of norbornane in its doubly ionized ground state.

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.

Aromaticity of giant polycyclic aromatic hydrocarbons with hollow sites: super ring currents in super-rings.

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 fate of dicationic states in molecular clusters of benzene and related compounds.

Calculations employing density functional theory indicate that, rather than undergoing fragmentation, dicationic clusters of benzene, hexafluorobenzene, and naphthalene produced by sequential one-electron or sudden double-ionization experiments on the neutrals can relax via the formation of inter-ring covalent C-C bonds, along with a series of proton transfers that enable a substantial reduction of inter- and intramolecular Coulomb repulsions. The theoretically predicted chemically bound structures correspond to deep local energy minima on the potential energy surface pertaining to the lowest electronic state of the dications and can therefore be regarded as metastable (kinetically long-lived) species. This discovery invalidates on theoretical grounds the liquid-droplet model of multiply charged clusters and sheds very unexpected light on possible consequences in chemistry of the intermolecular Coulombic decay (ICD) mechanism [Cederbaum, L. S.; et al. Phys. Rev. Lett. 1997, 79, 4778; Jahnke, T.; et al. Phys. Rev. Lett. 2004, 93, 163401] for deep inner-valence ionized states. Propagation of charge rearrangement reactions and proton transfers to several monomers may eventually lead to the formation of rather extended dicationic assemblies.

Probing the shape and stereochemistry of molecular orbitals in locally flexible aromatic chains: a penning ionization electron spectroscopy and Green's function study of the electronic structure of biphenyl.

We report on the results of an exhaustive study of the interplay between the valence electronic structure, the topology and reactivity of orbitals, and the molecular structure of biphenyl by means of Penning ionization electron spectroscopy in the gas phase upon collision with metastable He*(2(3)S) atoms. The measurements are compared with one-particle Green's function calculations of one-electron and shake-up valence ionization spectra employing the third-order algebraic diagrammatic construction scheme [ADC(3)]. Penning ionization intensities are also analyzed by means of the exterior electron-density model and comparison with photoelectron spectra: in contrast with the lines originating from sigma orbitals, ionization lines belonging to the pi-band system have large Penning ionization cross sections due to their greater extent outside the molecular van der Waals surface. The involved chemi-ionization processes are further experimentally investigated using collision-energy-resolved Penning ionization electron spectroscopy. The cross sections of pi-ionization bands exhibit a markedly negative collision-energy dependence and indicate that the interaction potential that prevails between the molecule and the He*(2(3)S) atom is strongly attractive in the pi-orbital region. On the other hand, the partial ionization cross sections pertaining to sigma-ionization channels are characterized by more limited collision-energy dependencies, as a consequence of rather repulsive interactions within the sigma-orbital region. A comparison of ADC(3) simulations with the Penning ionization electron spectra and UV photoelectron spectra measured by Kubota et al. [Chem. Phys. Lett. 1980, 74, 409] on thin films of biphenyl deposited at 170 and 109 K on copper demonstrates that biphenyl molecules lying at the surface of polycrystalline layers adopt predominantly a planar configuration, whereas within an amorphous sample most molecules have twisted structures similar to those prevailing in the gas phase.

Collision-energy-resolved Penning ionization electron spectroscopy of bromomethanes (CH3Br, CH2Br2, and CHBr3) by collision with He*(2(3)S) metastable atoms.

Ionization of bromomethanes (CH3Br, CH2Br2, and CHBr3) upon collision with metastable He*(2(3)S) atoms has been studied by means of collision-energy-resolved Penning ionization electron spectroscopy. Lone-pair (nBr) orbitals of Br4p characters have larger ionization cross sections than sigma(C-Br) orbitals. The collision-energy dependence of the partial ionization cross sections shows that the interaction potential between the molecule and the He*(2(3)S) atom is highly anisotropic around CH3Br or CH2Br2, while isotropic attractive interactions are found for CHBr3. Bands observed at electron energies of approximately 2 eV in the He*(2(3)S) Penning ionization electron spectra (PIES) of CH2Br2 and CHBr3 have no counterpart in ultraviolet (He I) photoionization spectra and theoretical (third-order algebraic diagrammatic construction) one-electron and shake-up ionization spectra. Energy analysis of the processes involved demonstrates that these bands and further bands overlapping with sigma(C-Br) or piCH2 levels are related to autoionization of dissociating (He+ - Br-) pairs. Similarly, a band at an electron energy of approximately 1 eV in the He*(2(3)S) PIES spectra of CH3Br has been ascribed to autoionizing Br** atoms released by dissociation of (unidentified) excited states of the target molecule. A further autoionization (S) band can be discerned at approximately 1 eV below the lone-pair nBr bands in the He*(2(3)S) PIES spectrum of CHBr3. This band has been ascribed to the decay of autoionizing Rydberg states of the target molecule (M**) into vibrationally excited states of the molecular ion. It was found that for this transition, the interaction potential that prevails in the entrance channel is merely attractive.

Theoretical study of the internal conversion of sulfoxide precursors of poly-isothianaphthene and related polymers.

In the present contribution, we theoretically investigate the suitability of the sulfoxide route for the synthesis of conjugated polymers of relevance for the fabrication of low-band gap materials with improved characteristics. The study focuses specifically on the internal elimination (E(i)) reactions of sulfoxide precursors of model oligomers of trans- and cis-poly-isothianaphtene (PITN), trans-poly-isothianaphtene vinylene (PITNV), and trans-poly-(ethylene dioxythiophene vinylene) (PEDOTV). These reactions have been characterized in detail by means of Density Functional Theory, along with the MPW1K functional (Modified Perdew-Wang 1-parameter model for kinetics).