Search

International audience; A periodic ab initio Hartme-Fock method has been used to evaluate a number of electronic and structural properties of MgHs. The calculated quantities are the equilibrium lattice parameters, the cohesive energy, the elastic constants, the energy-band struo ture, the density of' states, the electronic charge distribution and the Compton profiles. For the two former properties correlation effects are taken into account and improve the agreement with the experimental data.The strongly but not fitlly ionic character of MgHs is confirmed by the present study.

This study investigates the electronic configuration of titanium nitride-nickel (TiN-Ni) nanocomposites, in order to explain the high stability and activity of this hydrogenolysis catalyst. TiN-Ni is compared to a titanium oxide-nickel reference (TiO2-Ni). Strong metal support interactions are observed between the TiN and Ni. Scanning transmission electron microscopy (STEM) and energy-dispersive X-ray spectroscopy (EDS) illustrate that the Ni distributes more homogeneously on the nitride support. Computational comparison of TiN Ni and TiO2-Ni provides evidence of preferential Ni adsorption onto nitrogen sites of the nitride support. DFT calculations also predict a charge polarization between Ti and Ni atoms. X-ray photoelectron spectroscopy (XPS) corroborates computational analysis by revealing a suppression of surface nitride species upon deposition of Ni onto TiN. Shifts in Ti 2p and Ni 2p binding energy positions are also evident, which indicate an electronic perturbation between TiN and Ni. We conclude that the presence of nitrogen in the nitride support influences the electronic and structural properties of TiN-Ni and is partly responsible for the beneficial catalytic properties reported for this nanocomposite.

The structural, electronic and vibrational properties of two common defects in diamond, CHN and CHB, describing the case in which a carbon C1 atom is substituted by a nitrogen atom, or by a boron atom, breaking a C1–C2 bond, followed by the saturation of the dangling bond of C2 by a hydrogen atom, are investigated at the quantum mechanical level, by using a periodic supercell approach, hybrid DFT functionals and a local Gaussian-type basis set as implemented in the CRYSTAL code. The effect of concentration of the defects has been explored, by considering two supercells containing 64 and 216 atoms (S64 and S216). Formation and hydrogenation energies, geometries, Mulliken charges and the band structure of both defects are reported. The vibrational features of the defects have been investigated, by generating the IR and Raman spectra, and by analyzing graphically and through the isotopic substitution (H → D, 11B → 10B and 14N → 15N) the nature of the most relevant modes related to the defects. The computed C–H stretching mode of CHN, once corrected for anharmonicity (3408 cm $$^{-1}$$ ), falls to wavenumbers very close to the experimental peak observed at 3394 cm $$^{-1}$$ , which can then be reasonably attributed to this specific defect. The present manuscript is included in a special volume in honor and memory of Janos Angyan. Although he did not study in particular the kind of defects discussed in the present manuscript, the many methodological contributions he introduced in computational science have inspired many of the tools we have been using here. One of the present authors, RD, in particular, is grateful to Janos for the illuminating discussions they had in Paris, Nancy and Torino.

cited By 0; International audience; The VN3H defect in diamond (a vacancy surrounded by three nitrogen and one carbon atoms, the latter being saturated by a hydrogen atom) is investigated quantum-mechanically by use of a periodic supercell approach, an all-electron Gaussian-type basis set, "hybrid" functionals of density functional theory, and the Crystal program. Three fully optimized structural models (supercells containing 32, 64, and 128 atoms) are considered to investigate the effect of defect concentration. The electronic configuration of the defect is reported along with a description of its structural features. In particular, the influence of the lone-pair electrons of the three nitrogen atoms on the C-H bond is discussed. A thorough characterization of the vibrational spectroscopic features of the VN3H defect is also presented, where the anharmonicity of the most relevant normal modes is discussed. The infrared and Raman spectra show specific peaks, which allow for the identification of this particular defect among the many defects that are commonly present in both natural and irradiation-damaged diamonds. In particular, the main feature of the spectral fingerprint of the defect (i.e. the C-H stretching mode), experimentally observed at 3107 cm-1, is here computed at 3094 cm-1 with the B3LYP "hybrid" functional (with an anharmonic redshift of 157 cm-1 with respect to its harmonic value). The role played by the three nitrogen atoms on the spectral features of the defect is clearly identified through the redshift due to the 14N → 15N isotopic substitution. © 2017 the Owner Societies.

The adsorption of TiCl4 on the surfaces of MgCl2 crystals has been investigated by means of state-of-the-art periodic hybrid DFT methods, as the first step of a comprehensive study aiming to elucidate the structure of the active species in industrial MgCl2-supported Ziegler-Natta catalysts for ethene and propene polymerization. A first distinctive feature of the approach was the thorough evaluation of dispersion forces, crucial because the binding of TiCl4 on MgCl2 surfaces turned out to be essentially dispersion-driven. Also important was a careful investigation of the effects of different choices on basis set and density functional (DF) on the quantitative aspects of the results; this allowed us to trace the unusually large disagreement in the previous literature and identify unambiguous trends. In particular, three full sets of calculations were run adopting the B3LYP(-D), PBE0(-D) and M06 approximations: to the best of our knowledge, the last represents the first case of M06 functional implementation in a periodic code (CRYSTAL) of widespread use. The results consistently indicated that the adsorption of TiCl4 on well-formed MgCl2 crystals under conditions relevant for catalysis can only occur on MgCl2(1 1 0) or equivalent lateral faces, whereas the interaction with MgCl2(1 0 4) - for decades claimed as the most important surface in stereoselective catalysts - is too weak for the formation of stable adducts. The implications of these findings for catalysis are discussed. (C) 2011 Elsevier Inc. All rights reserved.

The present paper studies MX crystals in rock-salt structure (M: Li, Na, K; X: F, Cl, Br, I). They are often described as being formed by ions. Pictures based on quantum mechanical calculations sustain and quantify it. The tools used are (i) the Quantum Theory of Atoms in Molecules, (ii) the Electron Localization Function, and (iii) the maximization of the probability to find in a spatial domain a number of electrons equal to that of the ion under consideration. The present paper shows that the images provided by these three different tools to analyze the quantum mechanical calculations yield, for these systems, very similar results, in the sense that the spatial domains and probability distributions are close. While results for the first two methods are already present in the literature, the last of the methods is applied for the first time to these systems, and details about the method of calculation and program are also given.

The relative stability of nine different well defined CeO2 surfaces has been studied by periodic density functional calculations using GGA + U and B3LYP exchange-correlation functional. Both methods consistently predict that CeO2(111) is the most stable surface and also provide a consistent picture of the most stable surfaces which indeed are in agreement with previous studies based on empirical interatomic potentials. The facility of ceria surfaces to undergo a redox process has been investigated by forcing spin-polarized solutions, which lead to the occupancy of Ce 4f orbitals. These calculations provide evidence that surfaces with low-coordinated Ce cations are likely to be reduced more easily than regular low-index Miller surfaces. Fil: Branda, Maria Marta. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Instituto de Física del Sur. Universidad Nacional del Sur. Departamento de Física. Instituto de Física del Sur; Argentina Fil: Ferullo, Ricardo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Instituto de Química del Sur. Universidad Nacional del Sur. Departamento de Química. Instituto de Química del Sur; Argentina Fil: Causá, Mauro. Università degli Studi di Napoli Federico II; Italia Fil: Illas, Francesc. Universidad de Barcelona; España

PAH-based models, with an even or odd number of unsaturated carbon atoms and pi electrons (even and odd PAHs for short), are selected to investigate, by molecular and periodic methods, their electron distribution and border reactivity toward ozone, and also to represent local features and edge reactivity of even or odd soot platelets. These results will contrast those previously collected for the internal positions of similar even (J. Phys. Chem. A 2005, 109, 10929.) or odd systems (J. Phys. Chem. A 2008, 112, 973.). Topologically different peripheral positions, representative of armchair and zigzag borders, exhibit different reactivity right from the beginning. Ozone attacks start off either to give primary ozonides by concerted addition or, nonconcertedly, to first produce trioxyl intermediates. Then, a variety of pathways are described, whose viability depends on both model and position. They can open the way to the possible formation of epoxide, aldehyde, and phenol groups (all entailing O(2) production) or ether (+CO(2)), lactone (+H(2)CO), and ketone functionalities. To sum up, functionalization, regardless of how achieved, can give a number of groups, most of which actually observed in PAH ozonization experimental studies. This picture can be matched up to the results on internal sites of our preceding papers, for which epoxidation was the only outcome. Most interestingly, formation of a ketone group may turn an even system into an odd one (and conversely) while involving production of HOO(center dot).

In this work, the GLOB model, an effective and reliable computational approach well suited for ab initio and QM/MM molecular dynamics simulations of complex molecular systems in solution, has been applied to study two representative open-shell systems, the cobalt(II) ion and the glycine radical in aqueous solution, with special reference to their structural and magnetic properties. The main structural features of the solvent cage around the cobalt ion and the hydrogen bonding patterns around the neutral and zwitterionic forms of the glycine radical have been investigated in some detail. The general good agreement with experiments supports the use of the present model to investigate more challenging and biological/technological relevant open-shell systems.

"Classical" MgCl2-supported Ziegler-Natta catalysts (ZNCs) continue to dominate the industrial prodn. of isotactic polypropylene. There is a growing awareness of the inherent competitive edge of these low-cost systems over single-center (primarily metallocene) catalysts and of the potential for further improvement, particularly if deeper insight into the structure of the catalytic surfaces and the mechanisms of their modification by electron donors can be achieved. In the framework of a project ultimately aiming at the implementation of ZNCs with known and controlled surface structures, we are revisiting this whole area by using a combination of advanced computational (periodic DFT) and spectroscopic (high-resoln. magic-angle-spinning 1H NMR spectroscopy) tools. In this article, we report on the neat MgCl2 matrix and on model MgCl2/electron-donor adducts. The (104) surface, with five-coordinate Mg cations, is the dominant lateral termination in well-formed large crystals, and in highly activated MgCl2 samples prepd. by ball-milling. In the latter case, a minor fraction of surface Mg sites with a higher extent of coordinative unsatn. [e.g., four-coordinate Mg cations on (110) edges and/or at crystal corners or other defective locations] also appear to be present. RMe2Si(OMe) (R=octadecyl) binds to both types of Mg sites, albeit with different strengths resulting in different mobilities. The less-electron-donating RMeSi(OMe)2, in contrast, binds to the more unsatd. Mg sites only. The approach described herein is currently being extended to MgCl2/TiCln systems, and to their adducts with internal and external donors of different natures, strengths, and steric demands.

Graphite, and particularly defective graphite, is chosen to model soot particles. Quantum mechanical calculations are first carried out on molecular polycyclic-aromatic-hydrocarbon-type systems and then extended to a periodic representation of one graphite layer. The features of the interaction of H, HO, NO, NO2, and NO3 with these model systems are examined, with the aim of defining a suitable representation of the atmospheric or combustion gas−solid interactions by which functionalization reactions can take place. The more interesting interactions with small reactive molecules regard the edge of the graphene sheet and the in-plane carbon vacancies. While these interactions can be well described by sufficiently extended molecular models, periodic models are necessary to describe accurately the equilibrium geometries because they introduce the necessary geometric constraints. The ability of a graphene sheet to easily accommodate unpaired electrons in σ or π orbitals is the basis for its interesting intera...

Maximum Probability Domain (MPD) analysis has been recently applied to pure covalent and ionic crystals. The present study is devoted to a first MPD analysis of semi ionic crystals, Silicon Oxide, Aluminum Oxide and Titanium Oxide. These crystals are involved in important catalytic and photo-catalytic processes occurring on their surfaces. For this reason the study has been performed on bulk crystal and on surface slab models. Also surface neutral oxygen vacancy, the F-0 surface defect, has been considered. The Electron Localization Function (ELF) analysis has also been performed, due to its holistic approach to electronic structures. (C) 2015 Published by Elsevier BM.

5 pags.; 2 figs.; 2 tabs., In this work we propose a general strategy to calculate accurate He-surface interaction potentials. It extends the dispersionless density functional (dlDF) approach recently developed by Pernal et al. [Phys. Rev. Lett. 109, 263201 (2009)] to adsorbate-surface interactions by including periodic boundary conditions. We also introduce a scheme to parametrize the dispersion interaction by calculating two- and three-body dispersion terms at CCSD(T) level via the method of increments [Stoll, J. Chem. Phys. 97, 8449 (1992)]. The performance of the composite approach is tested on 4He/graphene by determining the energies of the low-lying selective adsorption states, finding an excellent agreement with the best available theoretical data. Second, the capability of the approach to describe dispersionless correlation effects realistically is used to extract dispersion effects in time-dependent density functional simulations on the collision of 4He droplets with a single graphene sheet. It is found that dispersion effects play a key role in the fast spreading of the 4He nanodroplet, the evaporation-like process of helium atoms, and the formation of solid-like helium structures. These characteristics are expected to be quite general and highly relevant to explain experimental measurements with the newly developed helium droplet mediated deposition technique. © 2014 AIP Publishing LLC, This work has been supported by Grants Nos. CCG08-CSIC/ESP-3680 from CSIC-CM, FIS2011-29596-C02-01 and FIS2011-28617- C02-01 from DGI, Spain (FEDER), and 2009SGR1289 from Generatitat de Catalunya.

The present study discusses the main features of the two programs CRYSTAL and EMBED developed for the ab initio study of the electronic properties of perfect periodic structures and of crystals with local defects, respectively. After a brief historical introduction, the structure of CRYSTAL is outlined and some specific aspects are discussed in detail: the use of local basis functions, the way of dealing with the Coulomb and the exchange series, the exploitation of point symmetry, the possibility to adopt either the Hartree–Fock approach or one among a variety of Kohn–Sham Hamiltonians. The present capabilities of the program are illustrated by a survey of selected applications from existing literature. Information is provided concerning work in progress aimed at removing some of the limitations of the code and improving its performance. The characteristics of EMBED are analyzed with emphasis given to the critical aspects of the method: limits of validity of the fundamental approximation on which the embedding technique relies, problems of convergence of the self-consistent procedure; and the delicate issue of the estimate of the defect formation energy. Again, a critical survey of applications clarifies the capabilities of the code. The envisaged improvements to be introduced in a forthcoming release of EMBED are presented. The present and prospective role of the two programs in the field of computational studies of condensed matter problems is outlined. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 1032–1048, 2000

By using the CRYSTAL95 program, ab initio periodic Hartree-Fock (PHF) calculations with the full potential and 6-21G* basis set are applied to interpretation of scanning tunneling microscope (STM) images on hexagonal graphite. Our results show asymmetry similar to the experimental and previous pseudopotential calculations. The dominant feature is the deep hollow in the middle of the carbon hexagon. The three carbon atoms that do not have a neighbor in the second layer appear as small hills, the other three with neighbors in the second layer simply appear as saddle points. Our calculation has been successfully used to reproduce experimental features such as the effect of increasing the magnitude of the bias voltage and the effect of increasing the tip-to-surface separation.

Vibrational spectroscopy and Density Functional Theory (DFT) quantum chemical calculations, have been used to investigate structure and chemical properties of some fluorinated ethers. A good agreement between experimental and theoretical results has been obtained. Interesting conformational behaviour has been observed and interpreted for C-H bonds in difluoromethyl end groups. These effects help in understanding the chemical behaviour of these compounds, since C-H equilibrium bond length and bond dissociation energy (BDE) can be correlated with vibrational properties. (C) 1998 Elsevier Science S.A. All rights reserved.

The capabilities of the CRYSTAL14 program are presented, and the improvements made with respect to the previous CRYSTAL09 version discussed. CRYSTAL14 is an ab initio code that uses a Gaussian-type basis set: both pseudopotential and all-electron strategies are permitted; the latter is not much more expensive than the former up to the first-second transition metal rows of the periodic table. A variety of density functionals is available, including as an extreme case Hartree-Fock; hybrids of various nature (global, range-separated, double) can be used. In particular, a very efficient implementation of global hybrids, such as popular B3LYP and PBE0 prescriptions, allows for such calculations to be performed at relatively low computational cost. The program can treat on the same grounds zero-dimensional (molecules), one-dimensional (polymers), two-dimensional (slabs), as well as three-dimensional (3D; crystals) systems. No spurious 3D periodicity is required for low-dimensional systems as happens when plane-waves are used as a basis set. Symmetry is fully exploited at all steps of the calculation; this permits, for example, to investigate nanotubes of increasing radius at a nearly constant cost (better than linear scaling!) or to perform self-consistent-field (SCF) calculations on fullerenes as large as (10,10), with 6000 atoms, 84,000 atomic orbitals, and 20 SCF cycles, on a single core in one day. Three versions of the code exist, serial, parallel, and massive-parallel. In the second one, the most relevant matrices are duplicated, whereas in the third one the matrices in reciprocal space are distributed for diagonalization. All the relevant vectors are now dynamically allocated and deallocated after use, making CRYSTAL14 much more agile than the previous version, in which they were statically allocated. The program now fits more easily in low-memory machines (as many supercomputers nowadays are). CRYSTAL14 can be used on parallel machines up to a high number of cores (benchmarks up to 10,240 cores are documented) with good scalability, the main limitation remaining the diagonalization step. Many tensorial properties can be evaluated in a fully automated way by using a single input keyword: elastic, piezoelectric, photoelastic, dielectric, as well as first and second hyperpolarizabilies, electric field gradients, Born tensors and so forth. Many tools permit a complete analysis of the vibrational properties of crystalline compounds. The infrared and Raman intensities are now computed analytically and related spectra can be generated. Isotopic shifts are easily evaluated, frequencies of only a fragment of a large system computed and nuclear contribution to the dielectric tensor determined. New algorithms have been devised for the investigation of solid solutions and disordered systems. The topological analysis of the electron charge density, according to the Quantum Theory of Atoms in Molecules, is now incorporated in the code via the integrated merge of the TOPOND package. Electron correlation can be evaluated at the Moller-Plesset second-order level (namely MP2) and a set of double-hybrids are presently available via the integrated merge with the CRYSCOR program. (C) 2014 Wiley Periodicals, Inc.

Our purpose is to identify a computational level sufficiently dependable and affordable to assess trends in the interaction of a variety of radical or closed shell unsaturated hydro-carbons A adsorbed on soot platelet models B. These systems, of environmental interest, would unavoidably have rather large sizes, thus prompting to explore in this paper the performances of relatively low-level computational methods and compare them with higher-level reference results. To this end, the interaction of three complexes between non-polar species, vinyl radical, ethyne, or ethene (A) with benzene (B) is studied, since these species, involved themselves in growth processes of polycyclic aromatic hydrocarbons (PAHs) and soot particles, are small enough to allow high-level reference calculations of the interaction energy ΔEAB. Counterpoise-corrected interaction energies ΔEAB are used at all stages. (1) Density Functional Theory (DFT) unconstrained optimizations of the A-B complexes are carried out, using the B3LYP-D, ωB97X-D, and M06-2X functionals, with six basis sets: 6-31G(d), 6-311 (2d,p), and 6-311++G(3df,3pd); aug-cc-pVDZ and aug-cc-pVTZ; N07T. (2) Then, unconstrained optimizations by Møller-Plesset second order Perturbation Theory (MP2), with each basis set, allow subsequent single point Coupled Cluster Singles Doubles and perturbative estimate of the Triples energy computations with the same basis sets [CCSD(T)//MP2]. (3) Based on an additivity assumption of (i) the estimated MP2 energy at the complete basis set limit [EMP2/CBS] and (ii) the higher-order correlation energy effects in passing from MP2 to CCSD(T) at the aug-cc-pVTZ basis set, ΔECC-MP, a CCSD(T)/CBS estimate is obtained and taken as a computational energy reference. At DFT, variations in ΔEAB with basis set are not large for the title molecules, and the three functionals perform rather satisfactorily even with rather small basis sets [6-31G(d) and N07T], exhibiting deviation from the computational reference of less than 1 kcal mol(-1). The zero-point vibrational energy corrected estimates Δ(EAB+ZPE), obtained with the three functionals and the 6-31G(d) and N07T basis sets, are compared with experimental D0 measures, when available. In particular, this comparison is finally extended to the naphthalene and coronene dimers and to three π-π associations of different PAHs (R, made by 10, 16, or 24 C atoms) and P (80 C atoms).

Most spectroscopic parameters are influenced by nuclear dynamics and by the chemical environment. However, proper inclusion of these effects still represents a challenge in computational spectroscopy studies. In many cases, a route coupling satisfactory accuracy with reasonable computational costs consists in the integration of DFT-based methods to compute spectroscopic parameters, with ab initio molecular dynamics simulations to sample from the classical phase space of the system. Here, we discuss the application of this approach in two case studies of remarkable practical interest, namely the simulation of the absorption spectrum of 9-methyladenine, an adenine nucleoside model; and the prediction of electron spin resonance parameters for nitroxyl radicals, the prototypical spin probes. In both cases, the accuracy of the results increases significantly when the subtle interplay of intra-molecular dynamics and solvent effects is introduced.

The 50 nm-thick polystyrene (PS) film, involved in some innovative memory devices, is sandwiched between two Al electrodes and contains 8-hydroxyquinoline (8HQ) molecules and gold nanoparticles. A periodic DFT study of the molecular adsorption of two different tautomers of the 8HQ on an Al (111) surface has been performed. The changes of the metal work function under molecular adsorption have been calculated to determine the possible effect at the interface on the electron and hole injection barriers. The strong adsorption of the 8HQ(II) tautomer makes the work function change up to 1.2 eV with respect to the bare Al (111) surface. Hence, the stabilization of the 8HQ(II) at the interface could play a key role in the switching mechanism of the device.

The electronic structure of different magnetic states of alpha, beta and gamma modifications of TiCl(3) has been computed employing the density functional theory with periodic boundary conditions and localized Gaussian basis sets. The analysis of the density of the electronic states (DOS) and of the spin density makes it possible to classify these halides as Mott-Hubbard insulators, where the band gap appears a result of large on-site Coulomb interaction. For each crystalline phase, the relative stability of different magnetic states has been analyzed in terms of exchange mechanisms. The electronic population data along with the spin density maps support the assumption of a d(1) Titanium ion in a distorted octahedral crystal field, notwithstanding the not fully ionic character of TiCl(3) modifications. Dispersion forces are particularly important for this material: a classical correction (of the type f(R)/R(6)) has been added to the DFT energies and gradients, providing a good agreement with structural data.

MgCl(2) is the preferred support for the industrial Ziegler-Natta catalysts, and is believed to act as a template for the epitactic chemisorption of the active Ti species. As the first step of a thorough computational modeling of these systems, we studied the bulk and surface structure of the ordered alpha and beta phases of MgCl(2) by means of periodic DFT (B3LYP) methods using localized basis sets. The layer structure of both phases was reproduced satisfactorily with the inclusion of a (small) empirical dispersion correction ("DFT-D") as a practical method to describe the attraction between the layers. Surface models were studied on slabs with adequate thickness. It appears that various surfaces exposing 5-coordinated Mg are very similar in energy and are the lowest non-trivial surfaces. Cuts exposing 4-coordinated Mg are significantly less stable; both kinetic and equilibrium models of crystal growth indicate that they should normally not be formed to a significant extent. "Nano-ribbons" of single, flat chains of MgCl(2), sometimes proposed as components of the disordered delta phase, were also evaluated, but are predicted to be unstable to rearrangement. Implications for the role of MgCl(2) as catalyst support are discussed.

Extensive DFT calculations are reported for the carbon centered glycine radical embedded in a-glycine crystal. Combined use of periodic and embedded cluster models leads to accurate structural and magnetic properties. Environmental effects are not negligible for all properties, whereas vibrational averaging effects connected to inversion at the radical center are significant only for C-alpha and H-alpha atoms. The good agreement between computed and experimental ESR parameters confirms the nature of one of the radicals obtained upon irradiation of a-glycine crystals and the reliability of DFT computations properly including both vibrational averaging and environmental effects. (c) 2007 Elsevier B.V. All rights reserved.

PAHs made from an odd number of unsaturated carbon atoms and pi electrons (odd PAHs) have been detected in flames and flank the more familiar even PAHs, having approximately the same quantitative importance, particularly for PAHs containing more than 25 carbon atoms. Similarly, soot platelets containing an odd number of carbon atoms can be reasonably assumed to form during combustion. PAHs are intended here as small models for the investigation of some, of their local features. To this end, quantum mechanical calculations were also carried out on periodic models. The spin density patterns were found to be highly dependent on the PAH size and shape. PAHs and soot, once released in the environment, can undergo several oxidation processes. Ozone is then taken as a probe of the reactivity properties of some internal exposed portions of a platelet. A primary ozonide (PO) corresponds to an energy minimum, but the relevant concerted addition pathway does not exist, because a PO-like saddle point is second-order. The reaction begins with a nonconcerted attack that produces a trioxyl radical (TR). Subsequent O-2 loss from the TR leaves either an epoxide with a pi-delocalized electron or a pi-delocalized oxepine, by cleavage of the ring carbon-carbon bond. The initial doublet spin multiplicity thus provides a description of the reaction surface unlike that for the internal reactivity of the closed-shell even systems investigated in a previous work, even though the final functionalization is the same.

The electron density, the electrostatic potential and the electric field of the MgO (100) surface, both regular and containing an oxygen vacancy (F s center), are compared in order to understand the modifications induced in the surface-absorbate interaction by the presence of the defect, with particular attention to the metal-oxide case. The spin-density for a gold atom absorbing on the most characteristic sites of the regular and F s -defected surface is also shown. It is found that in the defected surface the electron pair in the vacancy protrudes appreciably out of the surface, thus shifting the electrostatic potential to negative values (but producing a similar electric field) and being able to chemically interact with neighboring absorbed species. These results rationalize the rotational invariance and double frustration effects previously described for the metal/F s -defected MgO (100) surface.

The potentiality of the time dependent density functional theory (TDDFT) for the description of core excitation spectra (XAS) in transition metal oxides is analyzed, considering the rutile form of TiO(2) as a test case. Cluster models are adopted to mimic the bulk, embedded within an array of point charges to simulate the Madelung potential. All of the edges, titanium and oxygen K and titanium L edges, are considered, and the TDDFT results are compared with the experimental data in order to assess the performance of the theoretical approach in dealing with this complex class of compounds. Satisfactory results have been obtained for the Ti and O K edges, while in the case of the Ti L edge some discrepancies with the experiment are still present. The configuration mixing explicitly included in the TDDFT model strongly influences the distribution of the 2p metal oscillator strength. The origin of the spectral features is investigated with the help of the partial density of the virtual states (PDOS) calculated for each core hole considered, which can be qualitatively compared with the theoretical spectra calculated in the Kohn-Sham one-electron approach.

The desorption mechanism for oxygenated functionalities on soot is investigated by quantum mechanical calculations on functionalized polycyclic aromatic hydrocarbon PAH models and compared with recently published temperature programed desorption-mass spectrometry results. Substituents on PAHs of increasing size up to 46 carbon atoms in the parent PAH are chosen to reproduce the local features of an oxidized graphenic soot platelet. Initially, the study is carried out on unimolecular fragmentation extrusion, in some cases processes producing HO, CO, or CO2, in model ketones, carboxylic acids, lactones, anhydrides, in one aldehyde, one peroxyacid, one hydroperoxide, one secondary alcohol, and one phenol. Then, a bimolecular process is considered for one of the carboxylic acids. Furthermore, some cooperative effect which can take place by involving two vicinal carboxylic groups derived from anhydride hydrolysis is investigated for other four bifunctionalized models. The comparison between the computed fragmentation desorption barriers for the assessed mechanisms and the temperature at which maxima occur in TPD spectra for HO, CO, or CO2 desorption offers a suggestion for the assignment of these maxima to specific functional groups, i.e., a key to the description of the oxidized surface. Notably, the computations suggest that 1 the desorption mode from a portion of a graphenic platelet functionalized by a carboxylic or lactone groups is significantly dependent from the chemical and geometric local environment. Consequently, we propose that 2 not all carboxylic groups go lost at the relatively low temperatures generally stated, and 3 lactone groups can be identified as producing not only CO2 but also CO. © 2006 American Institute of Physics.

The ozonization mechanism for polycyclic aromatic hydrocarbons (PAHs) and soot is investigated by quantum mechanical calculations carried out on molecular and periodic systems. PAHs, interesting per se, serve also to model the local features of the graphenic soot platelets, for which another model is provided by a periodic representation of one graphenic layer. A concerted addition leads to a primary ozonide, while a nonconcerted attack produces a trioxyl diradical (in which one of the two unpaired electrons is pi-delocalized). Easy loss of (i) (1)O(2) or (ii) (3)O(2) from either intermediate, with spin conservation, would yield stable (i) singlet or (ii) triplet pi-delocalized species which carry an epoxide group. The trioxyl diradical pathway is estimated to be preferred, in these systems. An intersystem crossing, taking place in the trioxyl diradicals, can be invoked to allow the even easier loss of a ground-state oxygen molecule with the formation of a ground-state epoxide in a more exoergic and less demanding step. We propose that soot ozonization can take place by such a process, with ultimate functionalization of the graphenic platelets by epoxide groups.

The reaction rate constants of a hydroxyl radical with some hydrofluorocarbons and hydrofluoroethers have been calculated using “transition state theory”. A new computational method has been developed for transferring along the reaction coordinate the very high accuracy corrections, obtained at G2 level of theory, from small models (CH 4 and CF 3 H) to larger molecules (halogenated ethanes and methyl ethers). The substituent effects on the CH bond strengths and on the reaction rate constants have been examined, as well as the geometrical and electronic features of the transition states.

The atomic and electronic structure of the Ag/MgO(100) interface are calculated by means of the ab initio Hartree-Fock approach combined with a supercell model. The electronic density distribution and the interface binding energy/distance are analyzed for different Ag adsorption positions, slabs of different thicknesses and varying Ag surface coverage. It is demonstrated that the adhesion energy arises mainly due to the electrostatic interaction of substrate atoms with a complicated charge redistribution in the metal layer(s), characterized by large quadrupole moments as well as electron density redistribution towards bridge and hollow positions between the nearest and next-nearest Ag atoms. (C) 1998 Elsevier Science B.V.

We present an analysis of local or semilocal density functionals for the exchange)correlation energy by decomposing them into their gradients r local Seitz s . . . radius , z relative spin polarization , and s reduced density gradient . We explain the numerical method pertaining to this kind of analysis and present results for a few atoms and ions. The atomic shell structure is prominent, and only the ranges 0 - r - 10 and s 0 - s - 3 are important. The low-density and large-gradient domains, where the approximations for the exchange)correlation energy are least trustworthy, have very little weight. Q 1997 John Wiley & Sons, Inc.

We present electronic and structural parameters for bulk rutile (TiO2) determined by means of different theoretical methods, namely, the periodic Hartree-Fock approach, subsequent post-Hartree-Fock density functional correlation correction, and a linear-combination-of-atomic-orbitals approach based entirely on density functionals. Differences between the use of effective core potentials, a frozen-core treatment, and all-electron calculations are investigated. (C) 1996 John Wiley & Sons, Inc.

Density functional theory is taken into account to correct ‘a posteriori’ the total energy calculated with ab initio Hartree-Fock theory. Two models are adopted: the first one is the local density approximation parametrized by Ceperley and Alder while the second one is non-local and corresponds to the formulation proposed by Perdew. The changes in physical properties of ionic sulfide compounds (equilibrium geometry, formation energy, elastic constants, phase transition B1 → B3 in MgS and surface relaxation) associated with the electronic correlation are evaluated.

Numerous density functional theory formulas proposed in the literature are used to compute the total electronic energy of solids, integrating a posteriori functionals of electronic density, obtained with the periodic linear combination of the atomic orbitals Hartree-Fock method. Formation energies, equilibrium geometries and bulk moduli are calculated for a set of 18 covalent, ionic and semi-ionic crystals. A set of 11 density functional models is used. The calculated observables are compared with experimental data; a statistical analysis is performed to compare the models with greater ease.

Corundum, or alpha-Al2O3, is investigated at the ab initio periodic Hartree-Fock LCAO level using the program CRYSTAL. A split valence 6-21G basis set is adopted, and the influence of d polarization functions is explored. The four geometrical parameters of the unit cell are optimized, the difference with respect to the experimental data always being less than 1%. The electric field gradient tensor at both nuclei (Al and O) is calculated; agreement with the experimental data is quite satisfactory. Band structure, density of stages and charge density maps are also given. The resulting electronic structure is highly ionic.

Periodic pseudopotential Hartree–Fock calculations have been performed on the quartz structure of silica and germanium dioxide. It is shown that the tilt angle is the most important structural parameter. This parameter emphasizes the role of the oxygen polarization function which appear to be necessary to recover the correct geometry of low quartz. The structure of β‐quartz is rationalized on the basis of the energy dependence upon the tilt angle and the absence of β‐structure in GeO2 is explained. Finally the behavior of both oxides under pressure is discussed.

The electronic structure of crystalline urea (two molecules, 16 atoms per unit cell) is investigated at an ab initio level with CRYSTAL, a Hartree–Fock linear combination of atomic orbitals (LCAO) program for periodic systems. The influence of the basis set and of the computational parameters which control the treatment of the Coulomb and exchange series and the reciprocal space integration is documented; results include total and interaction energy, Mulliken analysis data and interaction (solid minus molecules) density maps, band structure, and density of states. The crystal field modifies the electronic structure of the isolated molecule, the main effect being an increase in the ionicity of bonds. The interaction energy obtained with a 6‐21** basis set is 28 kcal/mol, (16 kcal/mol after a correction of the basis set superposition error by using the counterpoise method) to be compared with 21±0.5 kcal/mol from experiment. This preliminary application shows that accurate ab initio calculations of hydrogen...

The electronic structure of corundum (α-Al 2 O 3 ) is calculated at the ab initio Hartree-Fock level. Cohesive energy, total and projected densities of states, atomic multipoles, bond populations and electron charge density distribution maps are given. The oxygen-aluminium bond is found to be partially covalent in nature; the atomic charges are −0.73 e and +1.09 e for O and Al respectively.

The core pseudopotentials of Durand and Barthelat have been introduced in the periodic Hartree-Fock program CRYSTAL. This new option allows the treatment of crystals involving heavy atoms and therefore extends the applicability of the method. A comparison between all-electron and pseudopotential results is presented. Binding energies, optimized cell parameters, bulk moduli, band structures and electron density curves calculated with the two methods are presented for covalent (diamond and silicon), metallic (lithium) and ionic (MgO) systems. As expected the pseudopotential results for covalent compounds are in good agreement with those obtained from all-electron calculations; for MgO larger discrepancies are observed, possibly due to core polarization effects.

Infrared spectroscopy in connection with quantum mechanical calculations based on density functional theory (DFT) method have been used for the characterization of fluorinated vinyl ethers. A good correspondence between experimental and calculated spectra has been found. The peculiar behaviour in terms of frequency and intensity parameters relative to the CC stretching normal mode along a series of related compounds has been analysed and interpreted in terms of atomic polar tensor (APT) and equilibrium charge and charge fluxes (ECCFs) models. This study indicates that the IR intensity of the CC stretching mode is almost completely due to the flux terms, suggesting that the presence of fluorine atoms make the charge along CC bond very mobile. This approach, based on comparative analysis of experimental and theoretical data, allows a critical comparison between various chemically related molecules and provides information on the various effects of chemical substitution on the molecular charge distribution.