Theoretical study of the oxidation mechanisms of naphthalene initiated by hydroxyl radicals: the O2 addition reaction pathways.

Atmospheric oxidation of the naphthalene-OH adduct [C10H8OH]˙ (R1) by molecular oxygen in its triplet electronic ground state has been studied using density functional theory along with the B3LYP, ωB97XD, UM05-2x and UM06-2x exchange-correlation functionals. From a thermodynamic viewpoint, the most favourable process is O2 addition at the C2 position in syn mode, followed by O2 addition at the C2 position in anti mode, O2 addition at the C4 position in syn mode, and O2 addition at the C4 position in anti mode, as the second, third and fourth most favourable processes. The syn modes of addition at these positions are thermodynamically favoured over the anti ones by the formation of an intramolecular hydrogen bond between the hydroxyl and peroxy substituents. Analysis of the computed structures, bond orders and free energy profiles demonstrate that the reaction steps involved in the oxidation of the naphthalene-OH adduct by O2 satisfy Hammond's principle. Kinetic rate constants and branching ratios under atmospheric pressure and in the fall-off regime have been supplied, using transition state and RRKM theories. By comparison with experiment, these data confirm the relevance of a two-step reaction mechanism. Whatever the addition mode, O2 addition in C4 position is kinetically favoured over O2 addition in C2 position, in contrast with the expectations drawn from thermodynamics and reaction energies. Under a kinetic control of the reaction, and in line with the computed reaction energy barriers, the most efficient process is O2 addition at the C4 position in syn mode, followed by O2 addition at the C2 position in syn mode, O2 addition at the C4 position in anti mode, and O2 addition at the C2 position in anti mode as the second, third and fourth most rapid processes. The computed branching ratios also indicate that the regioselectivity of the reaction decreases with increasing temperatures and decreasing pressures.

Benchmark theoretical study of the electric polarizabilities of naphthalene, anthracene, and tetracene.

By applying the principles of a focal point analysis to finite field calculations, the static electric dipole polarizabilities of naphthalene, anthracene, and tetracene have been estimated up to the level of coupled cluster theory including single, double, and perturbative triple excitations in the limit of an asymptotically complete basis set. Thermal corrections to polarizabilities have been empirically estimated according to Born-Oppenheimer molecular dynamical simulations at 298 K employing density functional theory. Confrontation with the latest experimental values for isotropic polarizabilities ultimately indicates relative accuracies of the order of 1.2% (for naphthalene) to 3.5% (for anthracene), provided diffuse functions are incorporated in the focal point analysis.

Half-metallicity and spin-contamination of the electronic ground state of graphene nanoribbons and related systems: an impossible compromise?

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 S(2) 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.

Quantum chemical study of conformational fingerprints in the photoelectron spectra and (e, 2e) electron momentum distributions of n-hexane.

The main purpose of the present work is to simulate from many-body quantum mechanical calculations the results of experimental studies of the valence electronic structure of n-hexane employing photoelectron spectroscopy (PES) and electron momentum spectroscopy (EMS). This study is based on calculations of the valence ionization spectra and spherically averaged (e, 2e) electron momentum distributions for each known conformer by means of one-particle Green's function [1p-GF] theory along with the third-order algebraic diagrammatic construction [ADC(3)] scheme and using Kohn-Sham orbitals derived from DFT calculations employing the Becke 3-parameters Lee-Yang-Parr (B3LYP) functional as approximations to Dyson orbitals. A first thermostatistical analysis of these spectra and momentum distributions employs recent estimations at the W1h level of conformational energy differences, by Gruzman et al. [J. Phys. Chem. A 2009, 113, 11974], and of correspondingly obtained conformer weights using MP2 geometrical, vibrational, and rotational data in thermostatistical calculations of partition functions beyond the level of the rigid rotor-harmonic oscillator approximation. Comparison is made with the results of a focal point analysis of these energy differences using this time B3LYP geometries and the corresponding vibrational and rotational partition functions in the thermostatistical analysis. Large differences are observed between these two thermochemical models, especially because of strong variations in the contributions of hindered rotations to relative entropies. In contrast, the individual ionization spectra or momentum profiles are almost insensitive to the employed geometry. This study confirms the great sensitivity of valence ionization bands and (e, 2e) momentum distributions on the molecular conformation and sheds further light on spectral fingerprints of through-space methylenic hyperconjugation, in both PES and EMS experiments.

A benchmark theoretical study of the electronic ground state and of the singlet-triplet split of benzene and linear acenes.

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 (C(4n+2)H(2n+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 (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 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-pVinfinityZ), the S(0)-T(1) 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-pVinfinityZ level for the S(0)-T(1) 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 S(0)-T(1) energy gap of 0 to approximately 4 kcal/mol (approximately 0.17 eV) in the limit of an infinitely large polyacene.

Study of the valence wave function of thiophene with high resolution electron momentum spectroscopy and advanced Dyson orbital theories.

Results of an exhaustive experimental study of the valence electronic structure of thiophene using high resolution electron momentum spectroscopy at impact energies of 1200 and 2400 eV are presented. The measurements were performed using an electron momentum spectrometer of the third generation at Tsinghua University, which enables energy, polar and azimuthal angular resolutions of the order of DeltaE = 0.8 eV, Deltatheta = +/-0.53 degrees and Deltaphi = +/-0.84 degrees . These measurements were interpreted by comparison with Green's function calculations of one-electron and shake-up ionization energies as well as of the related Dyson orbital electron momentum distributions, using the so-called third-order algebraic diagrammatic construction scheme (ADC(3)). Comparison of spherically averaged theoretical electron momentum distributions with experimental results very convincingly confirms the presence of two rather intense pi-2 pi*+1 shake-up lines at electron binding energies of 13.8 and 15.5 eV, with pole strengths equal to 0.18 and 0.13, respectively. Analysis of the electron momentum distributions associated with the two lowest 2A2 (pi3-1) and 2B1 (pi2-1) cationic states provides indirect evidence for a symmetry lowering and nuclear dynamical effects due to vibronic coupling interactions between these two states. ADC(3) Dyson orbital momentum distributions are systematically compared with distributions derived from Kohn-Sham (B3LYP) orbitals, and found to provide most generally superior insights into experiment.

A benchmark theoretical study of the electron affinities of benzene and linear acenes.

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 Moller-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 (C(2+4n)H(2+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 ( approximately -1.53 and approximately -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.

Imaging momentum orbital densities of conformationally versatile molecules: a benchmark theoretical study of the molecular and electronic structures of dimethoxymethane.

The main purpose of the present work is to predict from benchmark many-body quantum mechanical calculations the results of experimental studies of the valence electronic structure of dimethoxymethane employing electron momentum spectroscopy, and to establish once and for all the guidelines that should systematically be followed in order to reliably interpret the results of such experiments on conformationally versatile molecules. In a first step, accurate calculations of the energy differences between stationary points on the potential energy surface of this molecule are performed using Hartree-Fock (HF) theory and post-HF treatments of improving quality (MP2, MP3, CCSD, CCSD(T), along with basis sets of increasing size. This study focuses on the four conformers of this molecule, namely the trans-trans (TT), trans-gauche (TG), gauche-gauche (G+G+), and gauche-gauche (G+G-) structures, belonging to the C2v, C1, C2, and Cs symmetry point groups, respectively. A focal point analysis supplemented by suited extrapolations to the limit of asymptotically complete basis sets is carried out to determine how the conformational energy differences at 0 K approach the full CI limit. In a second step, statistical thermodynamics accounting for hindered rotations is used to calculate Gibbs free energy corrections to the above energy differences, and to evaluate the abundance of each conformer in the gas phase. It is found that, at room temperature, the G+G+ species accounts for 96% of the conformational mixture characterizing dimethoxymethane. In a third step, the valence one-electron and shake-up ionization spectrum of dimethoxymethane is analyzed according to calculations on the G+G+ conformer alone by means of one-particle Green's function [1p-GF] theory along with the benchmark third-order algebraic diagrammatic construction [ADC(3)] scheme. A complete breakdown of the orbital picture of ionization is noted at electron binding energies above 22 eV. A comparison with available (e,2e) ionization spectra enables us to identify specific fingerprints of through-space orbital interactions associated with the anomeric effect. At last, based on our 1p-GF/ADC(3) assignment of spectral bands, accurate and spherically averaged (e,2e) electron momentum distributions at an electron impact energy of 1200 eV are computed from the related Dyson orbitals. Very significant discrepancies are observed with momentum distributions obtained for several outer-valence levels using standard Kohn-Sham orbitals.

Probing molecular conformations in momentum space: the case of n-pentane.

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 degrees to 10 degrees 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.

Study of the molecular structure, ionization spectrum, and electronic wave function of 1,3-butadiene using electron momentum spectroscopy and benchmark Dyson orbital theories.

The scope of the present work is to reconcile electron momentum spectroscopy with elementary thermodynamics, and refute conclusions drawn by Saha et al. in J. Chem. Phys. 123, 124315 (2005) regarding fingerprints of the gauche conformational isomer of 1,3-butadiene in electron momentum distributions that were experimentally inferred from gas phase (e,2e) measurements on this compound [M. J. Brunger et al., J. Chem. Phys. 108, 1859 (1998)]. Our analysis is based on thorough calculations of one-electron and shake-up ionization spectra employing one-particle Green's function theory along with the benchmark third-order algebraic diagrammatic construction [ADC(3)] scheme. Accurate spherically averaged electron momentum distributions are correspondingly computed from the related Dyson orbitals. The ionization spectra and Dyson orbital momentum distributions that were computed for the trans-conformer of 1,3-butadiene alone are amply sufficient to quantitatively unravel the shape of all available experimental (e,2e) electron momentum distributions. A comparison of theoretical ADC(3) spectra for the s-trans and gauche energy minima with inner- and outer-valence high-resolution photoelectron measurements employing a synchrotron radiation beam [D. M. P. Holland et al., J. Phys. B 29, 3091 (1996)] demonstrates that the gauche structure is incompatible with ionization experiments in high-vacuum conditions and at standard temperatures. On the other hand, outer-valence Green's function calculations on the s-trans energy minimum form and approaching basis set completeness provide highly quantitative insights, within approximately 0.2 eV accuracy, into the available experimental one-electron ionization energies. At last, analysis of the angular dependence of relative (e,2e) ionization intensities nicely confirms the presence of one rather intense pi(-2) pi(*+1) satellite at approximately 13.1 eV in the ionization spectrum of the s-trans conformer.

Green's function study of the one-electron and shake-up ionization spectra of unsaturated hydrocarbon cage compounds.

The valence one-electron and shake-up ionization spectra of stella-2,6-diene, stella-2,6-dione, bicyclo-[2.2.2]-octane-2,5-dione, and bicyclo-[2.2.1]-heptane-2,5-dione have been exhaustively studied, up to the double ionization threshold and beyond, by means of one-particle Green's function theory. This study is based on calculations employing the outer-valence Green's function and the third-order algebraic diagrammatic construction schemes, along with a variety of basis sets. A comparison is made with available ultraviolet (He I) photoelectron and (e, 2e) electron-impact ionization spectra, with main focus on the identification of spectral fingerprints for cyclic strains and through-bond pi-conjugation. As a byproduct, our results demonstrate that it is impossible to reliably assign complex (e, 2e) ionization spectra by resorting only to Hartree-Fock or Kohn-Sham orbital energies and to the related electron momentum distributions.

The band 12 issue in the electron momentum spectra of norbornane: a comparison with additional Green's Function calculations and ultraviolet photoemission measurements.

In continuation of a recent study of the electronic structure of norbornane [J. Chem. Phys., 2004, 121, 10525] by means of electron momentum spectroscopy (EMS), we present Green's Function calculations of the ionization spectrum of this compound at the ADC(3) level using basis sets of varying quality, along with accurate evaluations at the CCSD(T) level of the vertical (26.5 eV) and adiabatic (22.1 eV) double ionization thresholds under C(2v) symmetry. The obtained results are compared with newly recorded ultraviolet photoemission spectra (UPS), up to binding energies of 40 eV. The theoretical predictions are entirely consistent with experiment and indicate that, in a vertical depiction of ionization, shake-up states at binding energies larger than approximately 26.5 eV tend to decay via emission of a second electron in the continuum. A band of s-type symmetry that has been previously seen at approximately 25 eV in the electron impact ionization spectra of norbornane is entirely missing in the UPS measurements and theoretical ADC(3) spectra. With regard to these results and to the time scales characterizing electron-electron interactions in EMS (10(-17) s) as compared with that (10(-13) s) of photon-electron interactions in UPS, and considering the p-type symmetry of the electron momentum distributions for the nearest 1b(1) and 1b(2) orbitals, this additional band can certainly not be due to adiabatic double ionization processes starting from the ground electronic state of norbornane, or to exceptionally strong vibronic coupling interactions between cationic states derived from ionization of the latter orbitals. It is therefore tentatively ascribed to autoionization processes via electronically excited and possibly dissociating states.

Investigation into the valence electronic structure of norbornene using electron momentum spectroscopy, Green's function, and density functional theories.

Results of a study of the valence electronic structure of norbornene (C(7)H(10)), up to binding energies of 30 eV, are reported. Experimental electron momentum spectroscopy (EMS) and theoretical Green's function and density functional theory approaches were utilized in this investigation. A stringent comparison between the electron momentum spectroscopy and theoretical orbital momentum distributions found that, among the tested models, the combination of the Becke-Perdew functional and a polarized valence basis set of triple-zeta quality provides the best representation of the electron momentum distributions for all 19 valence orbitals of norbornene. This experimentally validated model was then used to extract other molecular properties of norbornene (geometry, infrared spectrum). When these calculated properties are compared to corresponding results from independent measurements, reasonable agreement is typically found. Due to the improved energy resolution, EMS is now at a stage to very finely image the effective topology of molecular orbitals at varying distances from the molecular center, and the way the individual atomic components interact with each other, often in excellent agreement with theory. This will be demonstrated here. Green's Function calculations employing the third-order algebraic diagrammatic construction scheme indicate that the orbital picture of ionization breaks down at binding energies larger than about 22 eV. Despite this complication, they enable insights within 0.2 eV accuracy into the available ultraviolet emission and newly presented (e,2e) ionization spectra. Finally, limitations inherent to calculations of momentum distributions based on Kohn-Sham orbitals and employing the vertical depiction of ionization processes are emphasized, in a formal discussion of EMS cross sections employing Dyson orbitals.

Theoretical study of the conversion of sulfonyl precursors into chains of poly(p-phenylene vinylene).

The elimination and side reactions involved in the thermal conversion of sulfonyl precursor chains into poly(p-phenylene vinylene) (PPV) have been studied in detail, using Density Functional theory, along with the MPW1K functional. The performance of the MPW1K functional for describing radical dissociation and internal conversion reactions of sulfonyl precursors has been assessed against the results of benchmark CCSD(T) calculations. Enthalpies as well as entropies are calculated at different temperatures at the level of the rigid rotor-harmonic oscillator approximation. Entropy effects on internal elimination reactions are very limited. In sharp contrast, at the temperatures under which the conversion is usually performed (550 K), entropy contributions to the activation energies are found to be very significant and to strongly favor direct radical dissociations of the precursors. Further radical side reactions following an E(i) conversion through an alkyl substituent may also significantly contribute to the formation of sp(3) defects and/or cross-linked structures in the polymer-an advantageous feature for the making of materials with improved photoluminescence efficiencies.

Theoretical study of the internal elimination reactions of xanthate precursors.

The gas-phase internal elimination (E(i)) reaction of ethyl xanthate (CH(3)-CH(2)-S-CS-O-CH(3)) has been investigated by means of Hartree-Fock, second-order Møller-Plesset, and density functional theory (DFT) using the Becke three-parameter Lee-Yang-Parr (B3LYP) functional and the modified Perdew-Wang one-parameter model for kinetics (MPW1K). Considerable differences between the ground- and transition-state geometries and the calculated activation energies are observed from one approach to the other, which justifies first a careful calibration of the methods against the results of benchmark CCSD(T) calculations. Compared with these, DFT calculations along with the MPW1K functional are found to be an appropriate choice for describing the E(i) reaction of xanthate precursors. The precursor conformation and the transition states involved in the internal conversion of xanthate precursors of cyano derivatives of ethylene, and of cis- and trans-stilbene, are then characterized in detail by means of this functional.

Probing molecular conformations with electron momentum spectroscopy: the case of n-butane.

High-resolution (e,2e) measurements of the valence electronic structure and momentum-space electron density distributions of n-butane have been exhaustively reanalyzed in order to cope with the presence of two stable structures in the gas phase, namely the all-staggered and gauche conformers. The measurements are compared to a series of Boltzmann-weighted simulations based on the momentum-space form of Kohn-Sham (B3LYP) orbital densities, and to ionization spectra obtained from high-level [ADC(3)] one-particle Green's Function calculations. Indubitable improvements in the quality of the simulated (e,2e) ionization spectra and electron momentum profiles are seen when the contributions of the gauche form of n-butane are included. Both the one-electron binding energies and momentum distributions consistently image the distortions and topological changes that molecular orbitals undergo due to torsion of the carbon backbone, and thereby exhibit variations which can be traced experimentally. With regard to the intimate relation of (e,2e) cross sections with orbital densities, electron momentum spectroscopy can therefore be viewed as a very powerful, but up to now largely unexploited, conformational probe. The study also emphasizes the influence of thermal agitation in photoionization experiments of all kind.

Norbornane: an investigation into its valence electronic structure using electron momentum spectroscopy, and density functional and Green's function theories.

We report on the results of an exhaustive study of the valence electronic structure of norbornane (C(7)H(12)), up to binding energies of 29 eV. Experimental electron momentum spectroscopy and theoretical Green's function and density functional theory approaches were all utilized in this investigation. A stringent comparison between the electron momentum spectroscopy and theoretical orbital momentum distributions found that, among all the tested models, the combination of the Becke-Perdew functional and a polarized valence basis set of triple-zeta quality provides the best representation of the electron momentum distributions for all of the 20 valence orbitals of norbornane. This experimentally validated quantum chemistry model was then used to extract some chemically important properties of norbornane. When these calculated properties are compared to corresponding results from other independent measurements, generally good agreement is found. Green's function calculations with the aid of the third-order algebraic diagrammatic construction scheme indicate that the orbital picture of ionization breaks down at binding energies larger than 22.5 eV. Despite this complication, they enable insights within 0.2 eV accuracy into the available ultraviolet photoemission and newly presented (e,2e) ionization spectra, except for the band associated with the 1a(2) (-1) one-hole state, which is probably subject to rather significant vibronic coupling effects, and a band at approximately 25 eV characterized by a momentum distribution of "s-type" symmetry, which Green's function calculations fail to reproduce. We note the vicinity of the vertical double ionization threshold at approximately 26 eV.