Search

The accurate description of iron oxides/water interfaces requires reliable force field parameters that can be developed through comparison with sophisticated quantum mechanical calculations. Here, a set of CLASS2 force field parameters is optimized to describe the Fe–Owater cross-interaction through comparison with hybrid density functional theory (HSE06) calculations of the potential energy function for a single water molecule adsorbed on the Fe3O4 (001) surface and with density functional tight binding (DFTB+U) molecular dynamics simulations for a water trilayer on the same surface. The performance of the new parameters is assessed through the analysis of the number density profile of a water bulk (12 nm) sandwiched between two magnetite slabs of large surface area. Their transferability is tested for water adsorption on the curved surface of a spherical Fe3O4 nanoparticle of realistic size (2.5 nm). [ABSTRACT FROM AUTHOR]

In this work, we investigate the Fe3O4 (001) surface/water interface by combining several theoretical approaches, ranging from a hybrid functional method (HSE06) to density-functional tight-binding (DFTB) to molecular mechanics (MM). First, we assess the accuracy of the DFTB method to correctly reproduce HSE06 results on structural details and energetics and available experimental data for adsorption of isolated water, dimers, and trimers up to a water monolayer. Second, we build two possible configurations of a second and a third overlayer and perform molecular dynamics simulations with DFTB, monitoring the water orientation, the H-bond network, and the ordered water structure formation. To make our models more realistic, we then build a 12 nm-thick water multilayer on top of the Fe3O4 (001) surface slab model, which we investigate through MM-molecular dynamics (MD). The water layer structuring, revealed by the analysis of the atomic positions from a long MM-MD run for this large MM model, extends up to about 6–7 Å and nicely compares with that observed for a water trilayer model. However, MM and DFTB MD simulations show some discrepancy due to the poor description of the Fe⋯OH2 distance in MM that calls for further work in the parameterization of the model. [ABSTRACT FROM AUTHOR]

Magnetite has attracted increasing attention in recent years due to its promising and diverse applications in biomedicine. Theoretical modelling can play an important role in understanding magnetite-based nanomaterials at the atomic scale for a deeper insight into the experimental observations. However, calculations based on density functional theory (DFT) are too costly for realistically large models of magnetite nanoparticles. Classical force field methods are very fast but lack of precision and of the description of electronic effects. Therefore, a cheap and efficient quantum mechanical simulation method with comparable accuracy to DFT is highly desired. Here, a less computationally demanding DFT-based method, i.e., self-consistent charge density functional tight-binding (SCC-DFTB), is adopted to investigate magnetite bulk and low-index (001) surfaces with newly proposed parameters for Fe–O interactions. We report that SCC-DFTB with on-site Coulomb correction provides results in quantitatively comparable agreement with those obtained by DFT + U and hybrid functional methods. Therefore, SCC-DFTB is valued as an efficient and reliable method for the description magnetite. This assessment will promote SCC-DFTB computational studies on magnetite-based nanostructures that attract increasing attention for medical applications. [ABSTRACT FROM AUTHOR]

TiO2 nanoparticles (NPs) are nowadays considered fundamental building blocks for many technological applications. Morphology is found to play a key role with spherical NPs presenting higher binding properties and chemical activity. From the experimental point of view, the characterization of these nano-objects is extremely complex, opening a large room for computational investigations. In this work, TiO2 spherical NPs of different sizes (from 300 to 4000 atoms) have been studied with a two-scale computational approach. Global optimization to obtain stable and equilibrated nanospheres was performed with a self-consistent charge density functional tight-binding (SCC-DFTB) simulated annealing process, causing a considerable atomic rearrangement within the nanospheres. Those SCCDFTB relaxed structures have been then optimized at the DFT(B3LYP) level of theory. We present a systematic and comparative SCC-DFTB vs DFT(B3LYP) study of the structural properties, with particular emphasis on the surface-to-bulk sites ratio, coordination distribution of surface sites, and surface energy. From the electronic point of view, we compare HOMO-LUMO and Kohn-Sham gaps, total and projected density of states. Overall, the comparisons between DFTB and hybrid density functional theory show that DFTB provides a rather accurate geometrical and electronic description of these nanospheres of realistic size (up to a diameter of 4.4 nm) at an extremely reduced computational cost. This opens for new challenges in simulations of very large systems and more extended molecular dynamics. [ABSTRACT FROM AUTHOR]

We investigate the behavior of oxygen vacancies in three different metal-oxide semiconductors (rutile and anatase TiO2, monoclinic WO3, and tetragonal ZrO2) using a recently proposed hybrid density-functional method in which the fraction of exact exchange is material-dependent but obtained ab initio in a self-consistent scheme. In particular, we calculate charge-transition levels relative to the oxygen-vacancy defect and compare computed optical and thermal excitation/emission energies with the available experimental results, shedding light on the underlying excitation mechanisms and related materials properties. We find that this novel approach is able to reproduce not only ground-state properties and band structures of perfect bulk oxide materials but also provides results consistent with the optical and electrical behavior observed in the corresponding substoichiometric defective systems. [ABSTRACT FROM AUTHOR]

We investigate the long-standing problem of hole localization at the Al impurity in quartz SiO2, using a relatively recent DFT hybrid-functional method in which the exchange fraction is obtained ab initio, based on an analogy with the static many-body COHSEX approximation to the electron self-energy. As the amount of the admixed exact exchange in hybrid functionals has been shown to be determinant for properly capturing the hole localization, this problem constitutes a prototypical benchmark for the accuracy of the method, allowing one to assess to what extent self-interaction effects are avoided. We obtain good results in terms of description of the charge localization and structural distortion around the Al center, improving with respect to the more popular B3LYP hybrid-functional approach. We also discuss the accuracy of computed hyperfine parameters, by comparison with previous calculations based on other self-interaction-free methods, as well as experimental values. We discuss and rationalize the limitations of our approach in computing defect-related excitation energies in lowdielectric- constant insulators. [ABSTRACT FROM AUTHOR]

Transition metal doping of ZnO is considered as a promising way to obtain a diluted magnetic semiconducting oxide. In this work we investigate copper doping of ZnO by means of density functional theory, using a hybrid exchange-correlation functional and a periodic approach with localized atomic basis functions. Isolated copper species, such as copper substitutional to zinc, Cus, and Cu interstitial, Cui, are analyzed in terms of transition energy levels and hyperfine coupling constants with reference to available spectroscopic data. We also examine the potential magnetic interaction between copper species, their interaction with oxygen vacancies, and the possibility of copper clustering. The relative stability of the various copper impurities considered in this study is finally compared on the basis of their formation energy at different oxygen chemical potentials and Fermi level values. [ABSTRACT FROM AUTHOR]

A hybrid density functional study based on a periodic approach with localized atomic orbital basis functions has been performed in order to compute the optical and thermodynamic transition levels between different charge states of defect impurities in bulk ZnO. The theoretical approach presented allows the accurate computation of transition levels starting from single particle Kohn-Sham eigenvalues. The results are compared to previous theoretical findings and with available experimental data for a variety of defects ranging from oxygen vacancies, zinc interstitials, and hydrogen and nitrogen impurities. We find that H and Zn impurities give rise to shallow levels; the oxygen vacancy is stable only in the neutral VO and doubly charged VO2+ variants, while N-dopants act as deep acceptor levels. [ABSTRACT FROM AUTHOR]

The removal of lattice O atoms, as well as the addition of interstitial H atoms, in TiO2 is known to cause the reduction in the material and the formation of “Ti3+” ions. By means of electronic structure calculations we have studied the nature of such oxygen vacancy and hydrogen impurity states in the bulk of the anatase polymorph of TiO2. The spin polarized nature of these centers, the localized or delocalized character of the extra electrons, the presence of defect-induced states in the gap, and the polaronic distortion around the defect have been investigated with different theoretical methods: standard density functional theory (DFT) in the generalized-gradient approximation (GGA), GGA+U methods as a function of the U parameter, and two hybrid functionals with different admixtures of Hartree–Fock exchange. The results are found to be strongly dependent on the method used. Only GGA+U or hybrid functionals are able to reproduce the presence of states at about 1 eV below the conduction band, which are experimentally observed in reduced titania. The corresponding electronic states are localized on Ti 3d levels, but partly delocalized solutions are very close in energy. These findings show the limited predictive power of these theoretical methods to describe the electronic structure of reduced titania in the absence of accurate experimental data. [ABSTRACT FROM AUTHOR]

In this work we present a simplified procedure to use hybrid functionals and localized atomic basis sets to simulate scanning tunneling microscopy (STM) images of stoichiometric, reduced and hydroxylated rutile (110) TiO2 surface. For the two defective systems it is necessary to introduce some exact Hartree-Fock exchange in the exchange functional in order to correctly describe the details of the electronic structure. Results are compared to the standard density functional theory and planewave basis set approach. Both methods have advantages and drawbacks that are analyzed in detail. In particular, for the localized basis set approach, it is necessary to introduce a number of Gaussian function in the vacuum region above the surface in order to correctly describe the exponential decay of the integrated local density of states from the surface. In the planewave periodic approach, a thick vacuum region is required to achieve correct results. Simulated STM images are obtained for both the reduced and hydroxylated surface which nicely compare with experimental findings. A direct comparison of the two defects as displayed in the simulated STM images indicates that the OH groups should appear brighter than oxygen vacancies in perfect agreement with the experimental STM data. [ABSTRACT FROM AUTHOR]

The bonding, spin density related properties, and vibrational frequency of CO bound to single Au atom in the gas-phase or supported on MgO surfaces have been investigated with a variety of computational methods and models: periodic plane waves calculations have been compared with molecular approaches based on atomic orbital basis sets; pseudopotential methods with all electron fully relativistic calculations; various density functional theory (DFT) exchange-correlation functionals with the unrestricted coupled-cluster singles and doubles with perturbative connected triples [CCSD(T)]. AuCO is a bent molecule but the potential for bending is very soft, and small changes in the bond angle result in large changes in the CO gas-phase vibrational frequency. At the equilibrium geometry the DFT calculated vibrational shift of CO with respect to the free molecule is about -150 cm-1, whereas smaller values -60/-70 cm-1 are predicted by the more accurate CCSD(T) method. These relatively large differences are due to the weak and nonclassic bonding in this complex. Upon adsorption on MgO, the CO vibrational shift becomes much larger, about -290 cm-1, due to charge transfer from the basic surface oxide anion to AuCO. This large redshift is predicted by all methods, and is fully consistent with that measured for MgO/AuCO complexes. The strong influence of the support on the AuCO bonding is equally well described by all different approaches. [ABSTRACT FROM AUTHOR]

We report g tensors of surface color centers, so-called Fs+ centers, of MgO calculated with two density-functional approaches using accurately embedded cluster models. In line with recent UHV measurements on single-crystalline MgO film, we determined only small g-tensor anisotropies and negative shifts Δg≡g-ge for all Fs+ sites considered, namely, (001)-terrace, step, edge, and corner sites. The g values are very sensitive to the local structure of the defect: relaxation reverses the sign of Δg. However, accounting for the spin-orbit interaction either self-consistently or perturbatively yields very similar results. In addition to the values of the tensor components, their direction with respect to the surface was determined. In contrast to edges, significant deviations from ideal C2v symmetry were found for Fs+ centers at steps. Recent data on single-crystalline thin films are reevaluated in the light of these results. [ABSTRACT FROM AUTHOR]

The NO/NiO(100) system represents an excellent test case for the theory of surface chemical bond since accurate information about geometry, adsorption strength, and spin properties is available from experiments performed on NiO and Ni-doped MgO powders, single crystals, and thin films. We used cluster models to describe the NO/NiO interaction in combination with density functional theory (DFT) and wave function-based methods. We have identified four major aspects of the interaction: (1) the bonding cannot be described by a single determinant; (2) a spin-polarized DF-B3LYP approach gives reasonable adsorption properties at the price of a physically incorrect spin distribution; (3) a key ingredient of the interaction is the Coulomb repulsion within the Ni 3d shell; since this term is described very differently depending on the exchange-correlation functional it can result in overbound generalized gradient approach or Becke, Lee, Yang, and Parr or in strongly unbound (HFLYP) systems depending on the DFT approach; (4) the proper inclusion of the dynamical correlation is essential to treat the on-site Coulomb repulsion within the Ni 3d shell and to provide an accurate bond strength. In fact, the explicitly correlated complete-active-space second-order perturbation theory method gives results in overall agreement with the experiment. This shows the importance of treating on the same footing spin and electron correlation as well as the multiconfiguration character of the wave function. [ABSTRACT FROM AUTHOR]