Search

In this article, the results of a recently implemented DFT a posteriori and Kohn-Sham (Ks) linear combination of atomic orbital computational scheme for solids are presented. The equilibrium lattice parameters, bulk moduli, and lattice energies are calculated for eight crystallized systems. Local density approximation (LDA) and generalized gradient approximation (GGA) functionals and potentials are used. The maps of the Hartree-Fock (HF) and Ks electronic densities and band structures are depicted. The Ks results confirm the trend of the a posteriori scheme. Very good agreement between calculated and experimental lattice energies has been found for GGA potentials. (C) 1995 John Wiley & Sons, Inc.

The amplification of microscopic molecular motions so as to produce a controlled macroscopic body effect is the main challenge in the development of molecular mechanical devices. That amplification requires the coherent and ordered movement of each molecule of a whole macroscopic set, such as that taking place in a single-crystal-to-single-crystal transition. Actually, single-crystal-to-single-crystal transitions in molecular crystals can produce a variety of mechanical effects potentially useful in the development of smart materials. A challenging issue in these dynamic crystals, propedeutic to many possible applications in devices, is the gaining of a strict control over the mechanical effects associated with the transition. Here we report an example in which the control of the mechanical effects was successfully obtained. The compound studied undergoes a reversible single-crystal-to-single-crystal transition at 71 °C, from a planar stacked to a herringbone type packing. To this transition, a reversible macroscopic self-reshaping of the crystal is associated. Depending on the morphology, the crystal specimen undergoes a reversible longitudinal expansion of about 20% or a reversible transverse expansion of 20%, the other two dimensions of the crystal specimen being substantially unchanged. The amount of the macroscopic reshaping effect (20%) fully matches the relative variation of the sole unit cell parameter that changes during the transition (from 8.139 to 9.666 Å) in a sort of scale-invariant process. This represents striking evidence of controlled translation of sub-nanometer molecular motions up to the macroscopic scale of body motion.

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 present study investigates the coal char combustion by a combination of thermochemical and X-ray photoemission spectroscopy (XPS) analyses. Thermoanalytical methods (differential thermogravimetry, differential scanning calorimetry, and temperature-programmed desorption) are used to identify the key reactive steps that occur upon oxidation and heating of coal char (chemisorption, structural rearrangement and switchover of surface oxides, and desorption) and their energetics. XPS is used to reveal the chemical nature of the surface oxides that populate the char surface and to monitor their evolution throughout thermochemical processing. XPS spectra show the presence on the carbon surface of three main components. It is shown that the most abundant oxygen functionality in the raw char is epoxy. It decreases with preoxidation at 300 degrees C and even more at 500 degrees C, where carboxyl and ether oxygen functionalities prevail. The rearrangement of epoxy during preoxidation goes together with activation of the more stable and less reactive carbon sites. Results are in good agreement with semi-lumped kinetic models of carbon oxidation, which include (1) formation of "metastable" surface oxides, (2) complex switchover, and (3) desorption into CO and CO2.

The sole crystal structure generally observed for a chemical compound (actual crystal structure) corresponds to only one of the many minima present in the crystal energy landscape of that compound. A challenging problem of crystallography is to explain why the many other possible crystal structures (in terms of energy and density) present in the landscape (virtual structures) are not observed as polymorphs. In the present paper, with specific reference to semicarbazides of substituted benzoic acids, the concept of actual and virtual structures is extended and it is proven that optimization of a molecule in a different crystal packing of another molecule (interchangeability of molecules and packings) can lead to low minima of the lattice energy and to a high lattice density also in the case of molecules significantly different in shape (para/meta substituted compounds). This suggests that a crystal structure can be viewed as a sort of elastic system that can be deformed to fit molecules of different electronic character and shape. Surprisingly, in many cases, this deformation is accomplished at low cost of energy and density. Our results also suggest that the small number of polymorphs cannot be explained only by inspection of the crystal energy landscape, which would otherwise suggest many, but also, and perhaps mostly, by the analysis of the initial steps of crystal formation even before nucleation, when many different molecular clusters corresponding to different packings are in competition with each other.

Imines obtained by condensation of 4-hydroxybenzohydrazide with aliphatic ketones are a rare example of a class of compounds showing a remarkable tendency to crystallize in acentric polar space groups (Pna21 or Cc). In fact, all of the (seven) compounds studied up to now show at least one polar polymorph. In some cases, polymorphism was detected, and a nonpolar centrosymmetric phase was also identified (P21/c or P21/n space group). With the aim to disclose the conditions that can favor the formation of acentric structures in molecular crystals, we report, in this paper, a theoretical analysis (ab initio density functional theory with periodic boundary) of the lattice energy and density of all the packing modes observed in the whole set of imines. The computational analysis has been performed by optimizing each compound in its own experimental packings (actual crystal structures) and also in the packings of the other compounds of the class (virtual structures). The experimental crystallographic data and the theoretical analysis suggest that two conformers, basically differing for the orientation of the phenolic H atom in the plane of the phenyl ring, compete, in solution, for the formation of polar or centrosymmetric packings. The transitions between polar and centrosymmetric polymorphs are of diffusive type, and single crystals are not preserved, while the transitions between different polar polymorphs can be of single-crystal-to-single-crystal type.

Series of ZnII complexes with coordinated pyridine and a O,N,O tridentate core are reported. Four species of the general formula ZnLiPy2 were synthesized, where Li is a substituted ((aroylhydrazono)methyl)-3-hydroxyphenyl-4-(octyloxy)benzoate. The structure of the four mononuclear complexes and of dinuclear Zn2L24Py2 was solved by single crystals X-ray analysis. By TGA-DSC analysis and optical observation a nematogenic phase behavior was detected for ligands only and not for complexes. Photoluminescence spectra were recorded both in solution and in solid state Only the complexes showed a noteworthy photoluminescence in the solid state with tunable wavelength by varying the aroyl group. Medium to very high PL quantum yields were recorded on the crystalline complexes obtained in their mononuclear coordination environment. In the case of 4-nitrobenzoyl substituent a dinuclear complex with formula Zn2L24Py2 showed an unprecedent exceptionally high PLQY value. Both periodic and molecular DFT calculations were performed to rationalize the absorbance of these complexes at “one particle” level.

Reaction of N-rich conjugated bis(3,4-diamino-1,2,4-triazole)s with dilute hydrochloric acid affords bis(3,4-diamino-1,2,4-triazol-2-ium)chlorides. These compounds form layered structures in which planar layers of molecules are parallelly stacked. The π-stacking distance of the layers is relatively short, as compared with all-carbon-containing aromatic compounds, ranging between 3.00 and 3.22 Å at 173 K, and several contacts shorter than the sum of van der Waals radii are observed in the crystal structures. The features of the crystal structures are discussed in terms of the high nitrogen content of the compounds and of the H-bonding patterns. Periodic ab initio theoretical calculations of the crystal structures have allowed decomposing the lattice energy into various contributions in order to put up the relevance of van der Waals interactions for the π-stacking. In particular, it is found that van der Waals interactions account for about 10% of the total lattice energy and about 50% of the stacking energy (interlayer energy).

The ability of thermosalient solids, organic analogues of inorganic martensites, to move by rapid mechanical reconfiguration or ballistic event remains visually appealing and potentially useful, yet mechanistically elusive phenomenon. Here, with a material that undergoes both thermosalient and non-thermosalient phase transitions, we demonstrate that the thermosalient effect is preceded by anomalous thermal expansion of the unit cell. The crystal explosion occurs as sudden release of the latent strain accumulated during the anisotropic, exceedingly strong expansion of the unit cell with αa = 225.9 × 10−6 K−1, αb = 238.8 × 10−6 K−1 and αc = −290.0 × 10−6 K−1, the latter being the largest negative thermal expansivity observed for an organic compound thus far. The results point out to the occurence of the thermosalient effect in phase transitions as means to identify new molecular materials with strong positive and/or negative thermal expansion which prior to this work could only be discovered serendipitously.

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...

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

The role of surface oxygen vacancies in the optical properties of tin dioxide nanobelts is investigated in this paper. Using a first-principles approach, based on the density functional theory combined to a very accurate exchange correlation functional, we characterize SnO(2) (101), that is the nanobelt largest surface. We show that the presence of surface oxygen vacancies leads to the appearance of (i) occupied states located at about 1 eV above the valence band and (ii) unoccupied states lying in resonance with the conduction band. Photoluminescence characterization performed on samples of SnO(2) nanobelts at low temperature shows that the basic spectral features of luminescence are in excellent agreement with theoretical predictions.

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 wave functions of boron, aluminium and gallium phosphides have been calculated at both the Hartree - Fock and local density approximation levels using the ab initio linear combination of atomic orbitals method implemented in the program CRYSTAL. The Mulliken populations, bandstructure, electron charge densities, and structure factors are evaluated and compared. The chemical bonds are mainly covalent, as revealed by the difference electron charge-density maps obtained by subtracting from the total densities the superposition of the spherical atomic densities. A charge-density accumulation between the nearest neighbours slightly shifted towards the phosphorus appears in the AlP and GaP compounds, whereas an opposite charge transfer from phosphorus to boron is observed in boron phosphide from the Mulliken population analysis.

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).

The hydrate crystal lithium hydroxide monohydrate LiOH.H2O has been studied by ab initio periodic Hartree-Fock calculations. The influence of the crystalline environment on the local molecular properties (molecular geometry, atomic charges, electron density, molecular vibrations and deuterium quadrupole coupling constants) of the water molecule, the lithium and hydroxide ions has been calculated. A number of crystalline bulk properties are also presented, optimized crystalline structure, lattice energy and electronic band structure. The optimized cell parameters from calculations with a large basis set of triple-zeta quality differ by only 1-3% from the experimental neutron-determined cell, whereas the STO-3g basis set performs poorly (differences of 5-10%). With the triple-zeta basis also the atomic positions and intermolecular distances agree very well with the experiment. The lattice energy differs by approximately 8% from the experimental value, and by at most 3% when a density-functional electron correlation correction is applied. Large electron-density rearrangements occur in the water molecule and in the hydrogen bond and are in qualitative and quantitative agreement with experimental X-ray diffraction results. The quadrupole-coupling constants of the water and hydroxide deuterium atoms are found to be very sensitive to the O-H bond length and are in good agreement with experimental values when the calculation is based on the experimental structure. The anharmonic O-H stretching vibrations in the crystal are presented and found to be very close to results from calculations on molecular clusters. The electronic band and density-of-states spectra are discussed. Model calculations on a hydrogen fluoride chain were used to rationalize the results.

Calculations are reported of the enthalpies of mixing of MnO/MgO and NiO/MgO based on temperature-dependent lattice statics and dynamics and ab initio periodic Hartree-Fock methodology with a posteriori corrections for electron correlation derived from density-functional theory. For MnO/MgO both methodologies predict enthalpies of mixing that are close to those reported by Hahn and Muan (J. Phys. Chem. Solids, 1961, 19, 338). For NiO/MgO, on the other hand, classical simulations predict positive deviations from ideality in contrast to Hartree-Fock calculations which predict small negative deviations in agreement with the experimental data of Davies and Navrotsky (J. Solid State Chem., 1981, 38, 264).

Preliminary results of two annul simulations of the air quality in a region of the Northern Italy are presented. The simulations, which refer to two different year, were performed using a modeling system based on a meteorological model, an interface code and a transport-chemistry model. Comparisons with measured data are presented and discussed.

The effect of the Madelung field corresponding to formal point charges ± e arranged according to the fcc lattice of LiF upon CO, N2 and NO molecules is studied in order to isolate the purely electrostatic interaction from the other forces acting in the related physisorption processes. The polarization which is produced in those molecules by the Madelung field is analyzed in terms of dipole and quadrupole moments, as well as in terms of electron charge distributions. The influence of the electron correlation is taken into account by using the QCISD technique, which requires the semi-infinite lattice being simulated with a suitable cluster of point charges.

The static structural properties and B1-B2 phase transition of SrO were calculated at an ab initio level within the Hartree-Fock (HF) approximation using the effective core potential. The HF data were corrected a posteriori by integrating the HF charge density according to the correlation---only gradient-density-functional formulas proposed by Colle and Salvetti [Theor. Chim. Acta 37, 329 (1975)] and by Perdew [J. P. Perdew, Electronic Structure of Solids, edited by P. Ziesche and H. Eschring (Akademie Verlag, Berlin, in press)]. The HF binding energy, structural parameters, lattice parameter, bulk modulus, and elastic constants are in reasonable agreement with experiment. The correlation correction brings the binding energy into better agreement with the experiment. The transition pressure between B1 and B2 phases is in good agreement with the experimental value at the HF level, whereas it is corrected poorly by the addition of an a posteriori correlation correction.

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.

The adsorption of TiCl 4 on the surfaces of MgCl 2 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 MgCl 2-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 TiCl 4 on MgCl 2 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 TiCl 4 on well-formed MgCl 2 crystals under conditions relevant for catalysis can only occur on MgCl 2(1 1 0) or equivalent lateral faces, whereas the interaction with MgCl 2(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

In magnetic resonance imaging (MRI), the degree of contrast can be improved by using suitable contrast agents. The diethylenetriaminepentaacetic acid (DTPA) complexes of paramagnetic ions have been proposed for this purpose. This paper deals with extensive simulations of the distribution of species in blood plasma when solutions of manganese(II), copper(II), iron(III) and gadolinium(III) ions (as their soluble salts and as DTPA complexes) are injected. The various interaction equilibria, for which formation constants are known, are considered in order to assess toxic side-effects associated with their use in MRI. The data obtained support from the thermodynamic point of view, the use of GdDTPA, and suggest that the administration of a slight excess of ligand would guarantee complete coordination of the toxic gadolinium ion, and only minor interaction with the metal ions naturally present in blood plasma.

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 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 generalization of Becke’s numerical integration scheme to periodic functions is presented, which allows the LCAO-Kohn-Sham equations for crystals to be solved efficiently. The computational implementation of the scheme and its calibration are briefly discussed.

The preliminary results of a simulation of the air quality in a region of Northern Italy are presented. The simulation was performed using a modelling system based on a meteorological model, an interface code and a transport-chemistry model. Comparisons with measured data are presented and discussed.

The layered mineral brucite Mg(OH${)}_{2}$ is investigated theoretically using an ab initio all-electron linear combination of atomic orbitals Hartree-Fock (HF) approximation. At the HF level, the interlayer interaction is weak and the interlayer distance is larger than the experimental one. Bonding is discussed on the basis of density of states and charge-density maps. No hydrogen bond is characterized. A posteriori correction of the energy for the correlation error is performed by use of the functional approach. The three semilocal functional formulas used yield similar results. This brings in extra interlayer bonding interaction, and yields a calculated geometry in agreement with experiments. The analysis of the interlayer bondings shows that it is mainly of dispersion type, and that the used functionals account for dispersion, in particular at short interatomic distances.

The regular adsorption of CO at various concentrations on an infinite three-layer slab of LiF has been studied by using the CRYSTAL ab initio program. The CO molecules are adsorbed vertically above the surface cations by purely electrostatic attraction counterbalanced by the Pauli repulsion. The induced polarization has been analysed through the changes occurring in the molecular orbitals of the adsorbed species. The adsorption via the oxygen atom appears to be preferred, which contradicts IR experimental data. Such contradiction arises from the correlation error in the Hartree-Fock SCF method which is applied in the CRYSTAL program. Consideration of simple correlated clusters enables a semi-quantitative evaluation which additively corrects the above contradictory results.

Three different crystal phases of magnesium oxide [B1 (NaCl), B2 (CsCl), and B${8}_{1}$ (inverse NiAs)] are investigated theoretically. An ab initio all-electron linear combination of atomic orbitals Hartree-Fock approximation is adopted; extended basis sets are used which have been variationally optimized in each case. In ordinary conditions, the B${8}_{1}$ and B2 structures result less stable with respect to the rocksalt structure by 0.44 and 1.77 eV, respectively. With increasing pressure, the transition B1\ensuremath{\rightarrow}B2 is estimated to occur around 2.2 Mbar, while it is not excluded that a pressure interval exists around 2 Mbar where the B${8}_{1}$ phase is the most stable of the three structures. Data on the electronic properties of the three phases are provided and discussed: Mulliken populations, charge-density distribution, electron momentum distribution and anisotropy, magnesium core deformation, and band structure.

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 expressions for the corrections to the frequency of the3P2−3P1 and3P1−3P0 π-transitions of Mg and Ca are derived, taking into account the anomalous magnetic moment of the electron as well as the diamagnetic term of the Hamiltonian. Numerical results for24Mg are given.

The application of the maximum irreversible entropy variation principle to the calculation of the lattice constant is suggested as a new methodological approach.

A secondary pollution modelling system for simulating airborne dispersion and chemical reactions is applied over a regional scale domain located in the North-West of Italy, where urban and industrial areas are present. It was found an overestimation of NO(2) in the urban areas, probably due to an underestimation of the vertical diffusivity; the analysis of O(x) confirms that the discrepancies in O(3) and NO(2) are mainly due to local scale effects; the model shows a general underestimation of the observed PM10 concentrations due to the uncertainties in the emission inventories, spatial resolution, and the adopted aerosol modelling approach.

IR and Raman spectra of morpholine, a structure-directing agent, entrapped within the chabasite-related structures of both ALPO-34 and SAPO-34, are discussed. It is revealed that some vibrational features of the morpholine, which is present in the protonic form, are very sensitive to both host–guest (framework–morpholine) and guest–guest (morpholine–morpholine and morpholine–water) interactions occurring within the zeolites. FTIR spectroscopy shows that: (a) the NH 2 group of the protonated morpholine is H-bonded to the zeolitic frameworks (more strongly in SAPO-34); (b) the CH 2 vibrations in the 1500–1300 cm −1 region are sensitive to the surrounding being different for the ALPO-34 compared with the SAPO-34. Raman spectroscopy reveals that the ring-breathing mode of the morpholine has a small red-shift (from 822 cm −1 for the protonated molecule in solution, to 815 cm −1 in SAPO-34), and is deeply decreased in intensity when the molecule is in the constrained space of the zeolitic cages. In this case, two new modes at 472 and 245 cm −1 , tentatively assigned to morpholine–framework motions, are also present. IR and Raman bands are better resolved in the case of morpholine in ALPO-34 where the template–framework interactions are weaker.

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.