Search

Monolayer WSe 2 has aroused strong interest as a semiconducting 2D material. The tuning of its electronic properties is important for expanding the range of applications. In this work, by hybrid density functional theory calculations, we investigate the effect of Na adsorption on the electronic and magnetic properties of stoichiometric and defective WSe 2 monolayers and their WSe 2x S 2(1-x) alloys. Our results suggest that Na adsorption on defective WSe 2x S 2(1-x) monolayers would split and partially fill the empty defect states in the gap region inducing a local magnetic moment of 1 μ B per alkaline atom. The efficient tuning of the electronic and magnetic properties of WSe 2 monolayer by Na adsorption combined with the Se vacancy and S alloying would benefit its wide applications., (© 2020 IOP Publishing Ltd.)

In this work, we present a quantum mechanical investigation, based on the self-consistent charge density functional tight-binding (SCC-DFTB) method, of the functionalization with silane-type ligands (TETT) of a spherical TiO 2 nanoparticle of realistic size (2.2 nm containing 700 atoms) to create an efficient nanosystem for simultaneous photodynamic therapy and drug transport. We determine the mechanism of the TETT ligand anchoring and its stability under thermal treatment, through molecular dynamics simulations at 300 K. Then, we build a medium and a full coverage model (22 and 40 TETTs, respectively) and analyze the interaction among TETT ligands and between the ligands and the surface. Finally, on the fully covered nanoparticle, we succeed in localizing two minimum energy structures for an attached doxorubicin anticancer molecule (DOX) and provide the atomistic details for both the covalent and the non-covalent (electrostatic) types of interaction. A future development of this work will be the investigation of the loading capacity of this drug delivery system and of the pH effect of the surrounding aqueous environment., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

In this work, we investigate the Fe 3 O 4 (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 Fe 3 O 4 (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⋯OH 2 distance in MM that calls for further work in the parameterization of the model.

TiO 2 nanoparticles (NPs) are attracting materials for biomedical applications, provided that they are coated with polymers to improve solubility, dispersion and biocompatibility. Conformation, coverage density and solvent effects largely influence their functionality and stability. In this work, we use atomistic molecular dynamics simulations to study polyethylene glycol (PEG) grafting to highly curved TiO 2 NPs (2-3 nm) in different solvents. We compare the coating polymer conformations on NPs with those on (1 0 1) flat surfaces. In water, the transition from mushroom to brush conformation starts only at high density (σ = 2.25 chains/nm 2 ). In dichloromethane (DCM), at low-medium coverage (σ < 1.35 chains/nm 2 ), several interactions between the PEG chains backbone and undercoordinated Ti atoms are established, whereas at σ = 2.25 chains/nm 2 the conformation clearly becomes brush-like. Finally, we demonstrate that these spherical brushes, when immersed in water, but not in DCM, follow the Daoud-Cotton (DC) classical scaling model for the polymer volume fraction dependence with the distance from the center of star-shaped systems., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Iron oxide magnetic nanoparticles (NPs) are stimuli-responsive materials at the forefront of nanomedicine. Their realistic finite temperature simulations are a formidable challenge for first-principles methods. Here, we use density functional tight binding to open up the required time and length scales and obtain global minimum structures of Fe&#95;{3}O&#95;{4} NPs of realistic size (1400 atoms, 2.5 nm) and of different shapes, which we then refine with hybrid density functional theory methods to accomplish proper electronic and magnetic properties, which have never been accurately described in simulations. On this basis, we develop a general empirical formula and prove its predictive power for the evaluation of the total magnetic moment of Fe&#95;{3}O&#95;{4} NPs. By converting the total magnetic moment into the macroscopic saturation magnetization, we rationalize the experimentally observed dependence with shape. We also reveal interesting reconstruction mechanisms and unexpected patterns of charge ordering.

PEGylation of metal oxide nanoparticles is the common approach to improve their biocompatibility and in vivo circulation time. In this work, we present a combined experimental and theoretical study to determine the operating condition that guarantee very high grafting densities, which are desirable in any biomedical application. Moreover, we present an insightful conformational analysis spanning different coverage regimes and increasing polymer chain lengths. Based on 13 C NMR measurements and molecular dynamics simulations, we show that classical and popular models of polymer conformation on surfaces fail in determining the mushroom-to-brush transition point and prove that it actually takes place only at rather high grafting density values., Competing Interests: The authors declare no competing financial interest., (Copyright © 2019 American Chemical Society.)

Drug delivery systems are based on reversible interactions between carriers and drugs. Spacers are often introduced to tailor the type of interaction and to keep drugs intact. Here, we model a drug delivery system based on a functionalized curved TiO 2 nanoparticle of realistic size (700 atoms - 2.2 nm) by the neurotransmitter dopamine to carry the anticancer chemotherapeutic agent doxorubicin (DOX). The multiscale quantum chemical study aims at unraveling the nature and mechanism of the interactions between the components and the electronic properties of the composite system. We simulate the temperature effect through molecular dynamics runs of thermal annealing. Dopamine binds preferentially to low coordinated Ti sites on the nanoparticle through dissociated bidentate and chelate modes involving the diol groups. DOX is tethered by H-bonds, π-π stacking, dipole-dipole interactions and dispersion forces. Comparing different coverage densities of the spacer on the nanoparticle surface, we assess the best conditions for an effective drug transport and release: only at full coverage, DOX does not slip among the dopamine molecules to reach the nanoparticle surface, which is crucial to avoid the formation of stable coordinative bonds with under-coordinated Ti atoms. Finally, given the strong absorption properties and fluorescence of DOX and of the TiO 2 photocatalyst, we model the effect of light irradiation through excited state calculations to localize excitons and to follow the charge carrier's life path. This fundamental study on the nature and mechanism of drug/carrier interaction provides a solid ground for the rational design of new experimental protocols for a more efficient drug transport and release and its combination with photodynamic therapy.

Metal-doped and metal-supported two-dimensional materials are attracting a lot of interest as potentially active electrocatalysts for reduction and oxidation processes. Previously, when a non-regular 2 D h-BN layer was grown on a Cu(111) surface, metal adatoms were found to spontaneously emerge from the bulk to fill the atomic holes in the structure and become available for surface catalysis. Herein, computational electrochemistry is used to investigate and compare the performance of Cu-doped and Cu-supported pristine and defective h-BN systems for the electrocatalytic water oxidation reaction. For the various model systems, the intermediate species of this multistep oxidation process are identified and the free-energy variations for each step of reaction are computed, even for those steps that do not involve an electron or a proton transfer. Both associative and O 2 direct evolution mechanisms are considered. On this thermodynamic basis, the potential-determining step, the thermodynamic-determining step, and consequently the theoretical overpotential are determined for comparison with experiments. Small Cu clusters (tetramers) trapped in the h-BN defective lattice on a Cu(111) support are found to be very active for the water oxidation reaction since such systems are characterized by a low overpotential and by a small energy cost for O 2 release from the catalyst, which is often observed to be a major limit for other potential electrocatalysts., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Despite the many successful syntheses and applications of dopamine-functionalized TiO 2 nanohybrids, there has not yet been an atomistic understanding of the interaction of this 1,2-dihydroxybenzene derivative ligand with the titanium dioxide surfaces. In this work, on the basis of a wide set of dispersion-corrected hybrid density functional theory (DFT) calculations and density functional tight binding (DFTB) molecular dynamics simulations, we present a detailed study of the adsorption modes, patterns of growth, and configurations of dopamine on the anatase (101) TiO 2 surface, with reference to the archetype of 1,2-dihydroxybenzene ligands, i.e., catechol. At low coverage, the isolated dopamine molecule prefers to bend toward the surface, coordinating the NH 2 group to a Ti 5c ion. At high coverage, the packed molecules succeed in bending toward the surface only in some monolayer configurations. When they do, we observe a proton transfer from the surface to the ethyl-amino group, forming terminal NH 3 + species, which highly interact with the O atoms of a neighboring dopamine molecule. This strong Coulombic interaction largely stabilizes the self-assembled monolayer. On the basis of these results, we predict that improving the probability of dopamine molecules being free to bend toward the surface through thermodynamic versus kinetic growth conditions will lead to a monolayer of fully protonated dopamine molecules., Competing Interests: The authors declare no competing financial interest.

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.

Metal oxide nanoparticles (NPs) are regarded as good candidates for many technological applications, where their functional environment is often an aqueous solution. The correct description of metal oxide electronic structure is still a challenge for local and semilocal density functionals, whereas hybrid functional methods provide an improved description, and local atomic function-based codes such as CRYSTAL17 outperform plane wave codes when it comes to hybrid functional calculations. However, the computational cost of hybrids are still prohibitive for systems of real sizes, in a real environment. Therefore, we here present and critically assess the accuracy of our electrostatic embedding quantum mechanical/molecular mechanical (QM/MM) coupling between CRYSTAL17 and AMBER16, and demonstrate some of its capabilities via the case study of TiO₂ NPs in water. First, we produced new Lennard⁻Jones (LJ) parameters that improve the accuracy of water⁻water interactions in the B3LYP/TIP3P coupling. We found that optimizing LJ parameters based on water tri- to deca-mer clusters provides a less overstructured QM/MM liquid water description than when fitting LJ parameters only based on the water dimer. Then, we applied our QM/MM coupling methodology to describe the interaction of a 1 nm wide multilayer of water surrounding a spherical TiO₂ nanoparticle (NP). Optimizing the QM/MM water⁻water parameters was found to have little to no effect on the local NP properties, which provide insights into the range of influence that can be attributed to the LJ term in the QM/MM coupling. The effect of adding additional water in an MM fashion on the geometry optimized nanoparticle structure is small, but more evident effects are seen in its electronic properties. We also show that there is good transferability of existing QM/MM LJ parameters for organic molecules⁻water interactions to our QM/MM implementation, even though these parameters were obtained with a different QM code and QM/MM implementation, but with the same functional.

The structure and electronic properties of carbon-based nanostructures obtained by metal surface assisted synthesis is highly dependent on the nature of the precursor molecule. Here, we report on a combined scanning tunneling microscopy, soft X-ray spectroscopy and density functional theory investigation on the surface assisted polymerization of Br-corannulene at Ag(110) and on the possibility of building a mesh of π-conjugated polymers starting from buckyball shaped molecules. Indeed, the corannulene units form one-molecule-wide ribbons in which the natural concavity of the precursor molecule is maintained. These C-based nanostructures are corrugated and merge into a covalent network on the surface.

In most technological applications, nanoparticles are immersed in a liquid environment. Understanding nanoparticles/liquid interfacial effects is extremely relevant. This work provides a clear and detailed picture of the type of chemistry and physics taking place at the prototypical TiO 2 nanoparticles/water interface, which is crucial in photocatalysis and photoelectrochemistry. We present a multistep and multiscale investigation based on hybrid density functional theory (DFT), density functional tight-binding, and quantum mechanics/molecular mechanics calculations. We consider increasing water partial pressure conditions from ultra-high vacuum up to the bulk water environment. We first investigate single water molecule adsorption modes on various types of undercoordinated sites present on a realistic curved nanoparticle (2-3 nm) and then, by decorating all the adsorption sites, we study a full water monolayer to identify the degree of water dissociation, the Brønsted-Lowry basicity/acidity of the nanoparticle in water, the interface effect on crystallinity, surface energy, and electronic properties, such as the band gap and work function. Furthermore, we increase the water coverage by adding water multilayers up to a thickness of 1 nm and perform molecular dynamics simulations, which evidence layer structuring and molecular orientation around the curved nanoparticle. Finally, we clarify whether these effects arise as a consequence of the tension at the water drop surface around the nanosphere by simulating a bulk water up to a distance of 3 nm from the oxide surface. We prove that the nanoparticle/water interfacial effects go rather long range since the dipole orientation of water molecules is observed up to a distance of 5 Å, whereas water structuring extends at least up to a distance of 8 Å from the surface.

2D tungsten diselenide (2D-WSe₂) is one of the most successful bidimensional materials for optoelectronic and photonic applications, thanks to its strong photoluminescence properties and to a characteristic large exciton binding energy. Although these optical properties are widely recognized by the scientific community, there is no general understanding of the atomistic details of the excitonic species giving rise to them. In this work, we present a density functional theory investigation of excitons in 2D-WSe₂, where we compare results obtained by standard generalized gradient approximation (GGA) methods (including spin-orbit coupling) with those by hybrid density functionals. Our study provides information on the size of the self-trapped exciton, the number and type of atoms involved, the structural reorganization, the self-trapping energy, and the photoluminescence energy, whose computed value is in good agreement with experimental measurements in the literature. Moreover, based on the comparative analysis of the self-trapping energy for the exciton with that for isolated charge carriers (unbound electrons and holes), we also suggest a simplified approach for the theoretical estimation of the excitonic binding energy, which can be compared with previous estimates from different approaches or from experimental data.

Surfaces and their interaction with water play an important role in most of materials' applications. Magnetite has attracted continued interest in the fields of catalysis, spintronic devices, magnetic resonance imaging (MRI) and drug delivery. In this work, water adsorption and its effect on the stability diagram and on the electronic structure of the Fe3O4(001) surface are investigated by hybrid density functional theory calculations combined with an ab initio atomistic thermodynamic approach. We span a wide range of gaseous O2 and vapor H2O partial pressures. At low water pressure, a reconstructed SCV surface model is confirmed to be the most stable model at common working O2 partial pressures. However, at high water coverage, an unexpected stability inversion is observed that makes the hydrated bulk-terminated DBT surface the most favored. These results open up new horizons in Fe3O4 surface chemistry when working in an aqueous environment and are of key importance to develop rational strategies to surface engineering for high performance Fe3O4 nanomaterials.

Beyond two-dimensional (2D) materials, interfaces between 2D materials and underlying supports or 2D-coated metal or metal oxide nanoparticles exhibit excellent properties and promising applications. The hybrid interface between graphene and anatase TiO 2 shows great importance in photocatalytic, catalytic, and nanomedical applications due to the excellent and complementary properties of the two materials. Water, as a ubiquitous and essential element in practical conditions and in the human body, plays a significant role in the applications of graphene/TiO 2 composites for both electronic devices and nanomedicine. Carbon vacancies, as common defects in chemically prepared graphene, also need to be considered for the application of graphene-based materials. Therefore, the behavior of water on top and at the interface of defective graphene on anatase TiO 2 surface was systematically investigated by dispersion-corrected hybrid density functional calculations. The presence of the substrate only slightly enhances the on-top adsorption and reduces the on-top dissociation of water on defective graphene. However, at the interface, dissociated water is largely preferred compared with undissociated water on bare TiO 2 surface, showing a prominent cover effect. Reduced TiO 2 may further induce oxygen diffusion into the bulk. Our results are helpful to understand how the presence of water in the surrounding environment affects structural and electronic properties of the graphene/TiO 2 interface and thus its application in photocatalysis, electronic devices, and nanomedicine.

Understanding the electronic structure of metal oxide semiconductors is crucial to their numerous technological applications, such as photoelectrochemical water splitting and solar cells. The needed experimental and theoretical knowledge goes beyond that of pristine bulk crystals, and must include the effects of surfaces and interfaces, as well as those due to the presence of intrinsic defects (e.g. oxygen vacancies), or dopants for band engineering. In this review, we present an account of the recent efforts in predicting and understanding the optoelectronic properties of oxides using ab initio theoretical methods. In particular, we discuss the performance of recently developed dielectric-dependent hybrid functionals, providing a comparison against the results of many-body GW calculations, including G 0 W 0 as well as more refined approaches, such as quasiparticle self-consistent GW. We summarize results in the recent literature for the band gap, the band level alignment at surfaces, and optical transition energies in defective oxides, including wide gap oxide semiconductors and transition metal oxides. Correlated transition metal oxides are also discussed. For each method, we describe successes and drawbacks, emphasizing the challenges faced by the development of improved theoretical approaches. The theoretical section is preceded by a critical overview of the main experimental techniques needed to characterize the optoelectronic properties of semiconductors, including absorption and reflection spectroscopy, photoemission, and scanning tunneling spectroscopy (STS).

Heterogeneous photocatalysis is vital in solving energy and environmental issues that this society is confronted with. Although photocatalysts are often operated in the presence of water, it has not been yet clarified how the interaction with water itself affects charge dynamics in photocatalysts. Using water-coverage-controlled steady and transient infrared absorption spectroscopy and large-model (∼800 atoms) ab initio calculations, we clarify that water enhances hole trapping at the surface of TiO 2 nanospheres but not of well-faceted nanoparticles. This water-assisted effect unique to the nanospheres originates from water adsorption as a ligand at a low-coordinated Ti-OH site or through robust hydrogen bonding directly to the terminal OH at the highly curved nanosphere surface. Thus, the interaction with water at the surface of nanospheres can promote photocatalytic reactions of both oxidation and reduction by elongating photogenerated carrier lifetimes. This morphology-dependent water-assisted effect provides a novel and rational basis for designing and engineering nanophotocatalyst morphology to improve photocatalytic performances.

Magnetite exhibits a famous phase transition, called Verwey transition, at the critical temperature T V of about 120 K. Although numerous efforts have been devoted to the understanding of this interesting transition, up to now, it is still under debate whether a charge ordering and a band gap exist in magnetite above T V . Here, we systematically investigate the charge ordering and the electronic properties of magnetite in its cubic phase using different methods based on density functional theory: DFT+U and hybrid functionals. Our results show that, upon releasing the symmetry constraint on the density but not on the geometry, charge disproportionation (Fe 2+ /Fe 3+ ) is observed, resulting in a band gap of around 0.2 eV at the Fermi level. This implies that the Verwey transition is probably a semiconductor-to-semiconductor transition and that the conductivity mechanism above T V is small polaron hopping.

TiO 2 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, TiO 2 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 SCC-DFTB 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.

The photoactivity of methanol adsorbed on the anatase TiO 2 (101) surface was studied by a combination of scanning tunneling microscopy (STM), temperature-programmed desorption (TPD), X-ray photoemission spectroscopy (XPS), and density functional theory (DFT) calculations. Isolated methanol molecules adsorbed at the anatase (101) surface show a negligible photoactivity. Two ways of methanol activation were found. First, methoxy groups formed by reaction of methanol with coadsorbed O 2 molecules or terminal OH groups are photoactive, and they turn into formaldehyde upon UV illumination. The methoxy species show an unusual C 1s core-level shift of 1.4 eV compared to methanol; their chemical assignment was verified by DFT calculations with inclusion of final-state effects. The second way of methanol activation opens at methanol coverages above 0.5 monolayer (ML), and methyl formate is produced in this reaction pathway. The adsorption of methanol in the coverage regime from 0 to 2 ML is described in detail; it is key for understanding the photocatalytic behavior at high coverages. There, a hydrogen-bonding network is established in the adsorbed methanol layer, and consequently, methanol dissociation becomes energetically more favorable. DFT calculations show that dissociation of the methanol molecule is always the key requirement for hole transfer from the substrate to the adsorbed methanol. We show that the hydrogen-bonding network established in the methanol layer dramatically changes the kinetics of proton transfer during the photoreaction.

The presence of defects in the graphenic layers deposited on metal surfaces modifies the nature of the interaction. Unsaturated carbon atoms, due to vacancies in the lattice, form strong organometallic bonds with surface metal atoms that highly enhance the binding energy between the two materials. We investigate by means of a wide set of dispersion-corrected density functional theory calculations how such strong chemical bonds affect both the electronic properties of these hybrid interfaces and the chemical reactivity with water, which is commonly present in the working conditions. We compare different metal substrates (Cu vs Pt) that present a different type of interaction with graphene and with defective graphene. This comparative analysis allows us to unravel the controlling factors of water reactivity, the role played by the carbon vacancies and by the confinement or "graphene cover effect". Water is capable of breaking the C-Cu bond by dissociating at the undercoordinated carbon atom of the vacancy, restoring the weak van der Waals type of interaction between the two materials that allows for an easy detachment of graphene from the metal, but the same is not true in the case of Pt, where C-Pt bonds are much stronger. These conclusions can be used to rationalize water reactivity at other defective graphene/metal interfaces.

A water/(101) anatase TiO 2 interface has been investigated with the DFT-based self-consistent-charge density functional tight-binding theory (SCC-DFTB). By comparison of the computed structural, energetic, and dynamical properties with standard DFT-GGA and experimental data, we assess the accuracy of SCC-DFTB for this prototypical solid-liquid interface. We tested different available SCC-DFTB parameters for Ti-containing compounds and, accordingly, combined them to improve the reliability of the method. To better describe water energetics, we have also introduced a modified hydrogen-bond-damping function (HBD). With this correction, equilibrium structures and adsorption energies of water on (101) anatase both for low (0.25 ML) and full (1 ML) coverages are in excellent agreement with those obtained with a higher level of theory (DFT-GGA). Furthermore, Born-Oppenheimer molecular dynamics (MD) simulations for mono-, bi-, and trilayers of water on the surface, as computed with SCC-DFTB, evidence similar ordering and energetics as DFT-GGA Car-Parrinello MD results. Finally, we have evaluated the energy barrier for the dissociation of a water molecule on the anatase (101) surface. Overall, the combined set of parameters with the HBD correction (SCC-DFTB+HBD) is shown to provide a description of the water/water/titania interface, which is very close to that obtained by standard DFT-GGA, with a remarkably reduced computational cost. Hence, this study opens the way to the future investigations on much more extended and realistic TiO 2 /liquid water systems, which are extremely relevant for many modern technological applications.

The nature of the interaction of water with the WO 3 surface is of crucial importance for the use of this semiconductor oxide in photocatalysis. In this work, we investigate water adsorption and dissociation on both clean and O-deficient (001) WO 3 surfaces by means of an accurate DFT approach. The O vacancy formation energy (computed with respect to O 2 ) has been evaluated for all possible surface configurations, and the removal of the terminal O atom along the c axis is found to be preferred, costing about half the corresponding energy in the bulk. The presence of oxygen vacancies leads to a semiconductor to metal transition, confirming the experimental evidence of n-type conductivity in defective WO 3 films. H 2 O preferably adsorbs on WO 3 in a molecular undissociated form, due to the presence of W ions at the surface that act as Lewis acid sites. This interaction, about -1 eV per H 2 O molecule, is not very strong. Contrary to what is usually expected, the presence of oxygen vacancies does not significantly affect H 2 O adsorption. Finally, we investigated the H 2 O desorption from a hydroxylated surface. This suggests that the exposure of WO 3 to H 2 directly results in a hydroxylated surface and the corresponding H 2 O desorption turns out to be a very efficient mechanism to generate a reduced oxide surface, with important consequences on the electronic structure of this oxide.

By a combination of scanning tunneling microscopy, X-ray spectroscopic techniques and density functional theory calculations, we prove the formation of extended patterns of parallel, graphene nanoribbons with alternate zig-zag and armchair edges and selected width by surface-assisted Ullmann coupling polymerization and dehydrogenation of 1,6-dibromopyrene (C 16 H 8 Br 2 ). Besides the relevance of these nanostructures for their possible application in nanodevices, we demonstrate the peculiarity of halogenated pyrene derivatives for the formation of nanoribbons, in particular on Ag(110). These results open the possibility of tuning the shape and dimension of nanoribbons (and hence the correlated electronic properties) by choosing suitably tailored or on-purpose designed molecular precursors.

N-dopants in bulk monoclinic ZrO2 and their magnetic interactions have been investigated by DFT calculations, using the B3LYP hybrid functional. The electronic and magnetic properties of the paramagnetic N species, substitutionals and interstitials, are discussed. Their thermodynamic stability has been estimated as a function of the oxygen partial pressure. At 300 K, N prefers interstitial sites at any range of oxygen pressure, while at higher temperatures (700-1000 K), oxygen poor-conditions facilitate substitutional dopants. We have considered the interaction of two N defects in various positions in order to investigate the possible occurrence of ferromagnetic ordering. A very small magnetic coupling constant has been calculated for several 2N-ZrO2 configurations, thus demonstrating that magnetic ordering can be achieved only at very low temperatures, well below liquid nitrogen. Furthermore, when N atoms replace O at different sites, resulting in slightly different positions of the corresponding N 2p levels, a direct charge transfer can occur between the two dopants with consequent quenching of the magnetic moment. Another mechanism that contributes to the quenching of the N magnetic moments is the interplay with oxygen vacancies. These effects contribute to reduce the concentration of magnetic impurities, thus limiting the possibility to establish magnetic ordering.

The "catalysis under cover" involves chemical processes which take place in the confined zone between a 2D material, such as graphene, h-BN, or MoS2, and the surface of an underlying support, such as a metal or a semiconducting oxide. The hybrid interface between graphene and anatase TiO2 is extremely important for photocatalytic and catalytic applications because of the excellent and complementary properties of the two materials. We investigate and discuss the reactivity of O2 and H2O on top and at the interface of this hybrid system by means of a wide set of dispersion-corrected hybrid density functional calculations. Both pure and boron- or nitrogen-doped graphene are interfaced with the most stable (101) anatase surface of TiO2 in order to improve the chemical activity of the C-layer. Especially in the case of boron, an enhanced reactivity toward O2 dissociation is observed as a result of both the contribution of the dopant and of the confinement effect in the bidimensional area between the two surfaces. Extremely stable dissociation products are observed where the boron atom bridges the two systems by forming very stable B-O covalent bonds. Interestingly, the B defect in graphene could also act as the transfer channel of oxygen atoms from the top side across the C atomic layer into the G/TiO2 interface. On the contrary, the same conditions are not found to favor water dissociation, proving that the "catalysis under cover" is not a general effect, but rather highly depends on the interfacing material properties, on the presence of defects and impurities and on the specific reaction involved.

Graphene inertness towards chemical reactivity can be considered as an accepted postulate by the research community. This limit has been recently overcome by chemically and physically modifying graphene through non-metal doping or interfacing with acceptor/donor materials (metals or semiconductors). As a result, outstanding performances as catalytic, electrocatalytic, and photocatalytic material have been observed. In this critical Review we report computational work performed, by our group, on the reactivity of free-standing, metal- and semiconductor-supported B-doped graphene towards oxygen, which is at the basis of extremely important energy-related chemical processes, such as the oxygen reduction reaction. It appears that a combination of doping and interfacing approaches for the activation of graphene can open unconventional and unprecedented reaction paths, thus boosting the potential of modified graphene in many chemical applications., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The mechanism of water photooxidation on TiO2 surfaces is still controversial. Here we report a first-principles density functional study based on a hybrid functional method in which an adsorbed water molecule is found to directly interact with a self-trapped hole at a bridging oxygen site and to transform into an OH(•) radical species through a concerted proton/hole transfer. This study analyzes both the thermodynamics and kinetics of this step of the reaction, which is generally considered to be the rate determining one. The fate of the OH(•) radical is then investigated in terms of its reactivity with different surface species, with a second OH(•) radical, or with a second water molecule coming from the environment. We find that OH(•) radicals can either acquire a hydrogen from surrounding water molecules or, if they meet, couple to form hydrogen peroxide with highly associated energy gain.

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.

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 low-dielectric-constant insulators.

Dibromotetracene molecules are deposited on the Cu(110) surface at room temperature. The complex evolution of this system has been monitored at different temperatures (i.e., 298, 523, 673, and 723 K) by means of a variety of complementary techniques that range from STM and temperature-programmed desorption (TPD) to high-resolution X-ray spectroscopy (XPS) and near-edge X-ray absorption fine structure spectroscopy (NEXAFS). State-of-the-art density-functional calculations were used to determine the chemical processes that take place on the surface. After deposition at room temperature, the organic molecules are transformed into organometallic monomers through debromination and carbon-radical binding to copper adatoms. Organometallic dimers, trimers, or small oligomers, which present copper-bridged molecules, are formed by increasing the temperature. Surprisingly, further heating to 673 K causes the formation of elongated chains along the Cu(110) close-packed rows as a consequence of radical-site migration to the thermodynamically more stable molecule heads. Finally, massive dehydrogenation occurs at the highest temperature followed by ring condensation to nanographenic patches. This study is a paradigmatic example of how intermolecular coupling can be modulated by the stepwise control of a simple parameter, such as temperature, through a sequence of domino reactions., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

CONSPECTUS: Very rarely do researchers use metal oxides in their pure and fully stoichiometric form. In most of the countless applications of these compounds, ranging from catalysis to electronic devices, metal oxides are either doped or defective because the most interesting chemical, electronic, optical, and magnetic properties arise when foreign components or defects are introduced in the lattice. Similarly, many metal oxides are diamagnetic materials and do not show a response to specific spectroscopies such as electron paramagnetic resonance (EPR) spectroscopy. However, doped or defective oxides may exhibit an interesting and informative paramagnetic behavior. Doped and defective metal oxides offer an expanding range of applications in contemporary condensed matter science; therefore researchers have devoted enormous effort to the understanding their physical and chemical properties. The interplay between experiment and computation is particularly useful in this field, and contemporary simulation techniques have achieved high accuracies with these materials. In this Account, we show how the direct comparison between spectroscopic experimental and computational data for some selected and relevant materials provides ways to understand and control these complex systems. We focus on the EPR properties and electronic transitions that arise from the presence of dopants and defects in bulk metal oxide materials. We analyze and compare the effect of nitrogen doping in TiO2 and ZnO (two semiconducting oxides) and MgO (a wide gap insulator) and examine the effect of oxygen deficiency in the semiconducting properties of TiO2-x, ZnO1-x, and WO3-x materials. We chose these systems because of their relevance in applications including photocatalysis, touch screens, electrodes in magnetic random access memories, and smart glasses. Density functional theory (DFT) provides the general computational framework used to illustrate the electronic structure of these systems. However, for a more accurate description of the oxide band gap and of the electron localization of the impurity states associated with dopants or defects, we resorted to the use of hybrid functionals (B3LYP), where a portion of exact exchange in the exchange-correlation functional partly corrects for the self-interaction error inherent in DFT. In many cases, the self-interaction correction is very important, and these results can lead to a completely different physical picture than that obtained using local or semilocal functionals. We analyzed the electronic transitions in terms of their transition energy levels, which provided a more accurate comparison with experimental spectroscopic data than Kohn-Sham eigenvalues. The effects of N-doping were similar among the three oxides that we considered. The nature of the impurity state is always localized at the dopant site, which may limit their application in photocatalytic processes. Photocatalytic systems require highly delocalized photoexcited carriers within the material to effectively trigger redox processes at the surface. The nature of the electronic states associated with the oxygen deficiency differed widely in the three investigated oxides. In ZnO1-x and WO3-x the electronic states resemble the typical F-centers in insulating oxides or halides, with the excess electron density localized at the vacancy site. However, TiO2 acts as a reducible oxide, and the removal of neutral oxygen atoms reduced Ti(4+) to Ti(3+).

The dispersion of small amounts of Ce(4+) ions in the bulk of ZrO2 leads to a photoactive material sensitive to visible light. This is shown by monitoring with EPR the formation and the reactivity of photogenerated (λ > 420 nm) charge carriers. The effect, as confirmed by DFT calculations, is due to the presence in the solid of empty 4f Ce states at the mid gap, which act as intermediate levels in a double excitation mechanism. This solid can be considered an example of a third-generation photoactive material.

An important step in oxide photochemistry, the loading of electrons into shallow trap states, was studied using infrared (IR) spectroscopy on both, rutile TiO2 powders and single-crystal, r-TiO2(110) oriented samples. After UV-irradiation or n-doping by exposure to H-atoms broad IR absorption lines are observed for the powders at around 940 cm(-1). For the single crystal substrates, the IR absorption bands arising from an excitation of the trapped electrons into higher-lying final states show additional features not observed in previous work. On the basis of our new, high-resolution data and theoretical studies on the polaron binding energy in rutile we propose that the trap states correspond to polarons and are thus intrinsic in nature. We assign the final states probed by the IR-experiments to hydrogenic states within the polaron potential. Implications of these observations for photochemistry on oxides will be briefly discussed.

Understanding the interfacial hole transfer between photoirradiated anatase TiO2 nanoparticles and organic adsorbates at the surface is a crucial step toward the full comprehension and control of the photooxidation processes, which are fundamental in photocatalytic applications for energy production, environmental care, and medical cancer phototherapy. In this density functional study, we use the hybrid functional B3LYP in order to investigate how the hole transfer process takes place at the anatase (101) surface when a series of organic alcohols and acids are adsorbed. The results allow us to propose a mechanism of photoinduced proton dissociation and a scale of scavenging power of the different organic adsorbates (glycerol > tert-butanol > iso-propanol > methanol > formic acid). We observe that the surface dipole originated by the molecular adsorption causes a reduction in the cost to form a hole at a surface oxygen. This can then be trapped by the organic adsorbate only if proton dissociation takes place. In the case of glycerol, the hole transfer triggers further photoinduced chemical reactivity than the simple proton dissociation. Only in the case of cathecol, which is an excellent hole scavenger, the proton dissociation is not a required step.

(WO(3))(3) clusters deposited on the (110) rutile TiO(2) surface are excellent catalysts for the formaldehyde (CH(2)O) polymerization reaction (J. Phys. Chem. C 2010, 114, 17017). The present B3LYP study unravels the possible paths of this catalyzed reaction. According to the stoichiometry of the r-TiO(2) surface, the (WO(3))(3) clusters can be neutral, singly charged, or doubly charged. We find that only neutral (WO(3))(3) and anionic (WO(3))(3)(-) clusters are reactive toward CH(2)O molecules. In both cases it is possible to determine more than one mechanism on the basis of a nucleophilic attack of the formaldehyde O atom to the W ions of the cluster. The reaction proceeds through successive attacks of other CH(2)O molecules and the formation of acetal and polyacetal intermediates, which inhibits the chain propagation. Only in the case of the anionic (WO(3))(3)(-) catalyst is a totally different reaction path possible at low temperatures. This path involves the formation of radical species where the unpaired electron is localized on the organic moiety bound to the cluster. The polymer chain propagation follows a radical mechanism with low activation barriers. Thus, a cluster's electron charging speeds up the formaldehyde polymerization at low temperatures. On the basis of these unexpected results, we conclude that electron-rich supports and low working temperatures are the keys to kinetic control of the reaction favoring a fast radical chain propagation mechanism.

It has been experimentally observed that Ti doping of bulk ZrO(2) induces a large red-shift of the optical absorption edge of the material from 5.3 to 4.0 eV [Livraghi et al., J. Phys. Chem. C, 2010, 114, 18553-18558]. In this work, density functional calculations based on the hybrid functional B3LYP show that Ti dopants in the substitutional position to Zr in the tetragonal lattice cause the formation of an empty Ti 3d band about 0.5 eV below the bottom of the conduction band. The optical transition level ε(opt)(0/-1) from the topmost valence state to the lowest empty Ti impurity state is found at 4.9 eV in a direct band gap of 5.7 eV. The calculated shift is consistent with the experimental observation. The presence of Ti(3+) species in Ti-doped ZrO(2), probed by means of electron paramagnetic resonance (EPR), is rationalized as the result of electron transfers from intrinsic defect states, such as oxygen vacancies, to substitutional Ti(4+) centers., (This journal is © the Owner Societies 2011)

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, Cu(s), and Cu interstitial, Cu(i), 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.

The (17)O and (1)H hyperfine interactions of water ligands in the Ti(III) aquo complex in a frozen solution were determined using Hyperfine Sublevel Correlation (HYSCORE) and Pulse Electron Nuclear Double Resonance (ENDOR) spectroscopies at 9.5 GHz. The isotropic hyperfine interaction (hfi) constant of the water ligand (17)O was found to be about 7.5 MHz. (1)H Single Matched Resonance Transfer (SMART) HYSCORE spectra allowed resolution of the hfi interactions of the two inequivalent water ligand protons and the relative orientations of their hfi tensors. The magnetic and geometrical parameters extracted from the experiments were compared with the results of DFT computations for different geometrical arrangements of the water ligands around the cation. The theoretical observable properties (g tensor (1)H and (17)O hfi tensors and their orientations) of the [Ti(H(2)O)(6)](3+) complex are in quantitative agreement with the experiments for two slightly different geometrical arrangements associated with D(3d) and C(i) symmetries.

A combined experimental and first principles study of Cu defects in bulk ZnO is presented. Cu particles are epitaxially deposited on the polar O-ZnO(0001) surface at room temperature. Upon heating, a broadening of the quasielastic peak in high resolution electron energy loss spectra is observed, corresponding to an electronic doping effect of Cu atoms in bulk ZnO with an ionization energy of 88 meV. Cu impurities in ZnO, although commonly acting as acceptors, are presently observed to induce shallow donor states. We assign these to interstitial Cu species on the basis of a hybrid density functional study.

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 V(O) and doubly charged V(O) (2+) variants, while N-dopants act as deep acceptor levels.

The formation of O(2) (-) radical anions by contact of O(2) molecules with a Na pre-covered MgO surface is studied by a combined EPR and quantum chemical approach. Na atoms deposited on polycrystalline MgO samples are brought into contact with O(2). The typical EPR signal of isolated Na atoms disappears when the reaction with O(2) takes place and new paramagnetic species are observed, which are attributed to different surface-stabilised O(2) (-) radicals. Hyperfine sublevel correlation (HYSCORE) spectroscopy allows the superhyperfine interaction tensor of O(2) (-)Na(+) species to be determined, demonstrating the direct coordination of the O(2) (-) adsorbate to surface Na(+) cations. DFT calculations enable the structural details of the formed species to be determined. Matrix-isolated alkali superoxides are used as a standard to enable comparison of the formed species, revealing important and unexpected contributions of the MgO matrix in determining the electronic structure of the surface-stabilised Na(+)-O(2) (-) complexes.

The removal of lattice O atoms, as well as the addition of interstitial H atoms, in TiO(2) is known to cause the reduction in the material and the formation of "Ti(3+)" 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 TiO(2). 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.

The adsorption of small amounts of alkali metal atoms (Li, Na, K, Rb, and Cs) on the surface of MgO powders and thin films has been studied by means of EPR spectroscopy and DFT calculations. From a comparison of the measured and computed g values and hyperfine coupling constants (hfccs), a tentative assignment of the preferred adsorption sites is proposed. All atoms bind preferentially to surface oxide anions, but the location of these anions differs as a function of the deposition temperature and alkali metal. Lithium forms relatively strong bonds with MgO and can be stabilized at low temperatures on terrace sites. Potassium interacts very weakly with MgO and is stabilized only at specific sites, such as at reverse corners where it can interact simultaneously with three surface oxygen atoms (rubidium and cesium presumably behave in the same way). Sodium forms bonds of intermediate strength and could, in principle, populate more than a single site when deposited at room temperature. In all cases, large deviations of the hfccs from the gas-phase values are observed. These reductions in the hfccs are due to polarization effects and are not connected to ionization of the alkali metal, which would lead to the formation of an adsorbed cation and a trapped electron. In this respect, hydrogen atoms behave completely differently. Under similar conditions, they form (H(+))(e(-)) pairs. The reasons for this different behavior are discussed.

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.