Toward an efficient f-in-core/f-in-valence switchable description for DFTB calculations of Ce 4f states in ceria.

A computational protocol is developed for efficient studies of partially reduced redox-active oxides using the self-consistent charge density functional tight-binding method. The protocol is demonstrated for ceria, which is a prototypical reducible oxide material. The underlying idea is to achieve a consistent (and harmonized) set of Slater-Koster (SK) tables with connected repulsive potentials that enable switching on and off the in-valence description of the Ce 4f states without serious loss of accuracy in structure and energetics. The implicit treatment of the Ce 4f states, with the use of f-in-core SK-tables, is found to lead to a significant decrease in computational time. More importantly, it allows for explicit control of the oxidation states of individual Ce atoms. This makes it possible to "freeze" the electronic configuration, thereby allowing the exploration of the energetics for various meta-stable configurations. We anticipate that the outlined strategy can help to shed light on the interplay between the size, shape, and redox activity for nanoceria and other related materials.

Water on ceria{111}: Comparison between 23 experimental vibrational studies in the literature and new modeling.

Theoretical and experimental vibrational signatures of H2O and OH- (dissociated water) adsorbed on stoichiometric ceria{111} surfaces are compared. The experimental ones were collected from low-coverage experiments in the literature, and the theoretical anharmonic frequencies were generated using density functional theory calculations employing the optPBE-vdW functional for coverages from 0.5 to a few monolayers. It is found that (i) the experiments and our calculations overall agree well, lending credibility to both; (ii) the calculations manage to resolve the large class of H-bonded motifs into frequency classes that can guide experimental assignments; (iii) it is possible to find a geometrical H-bond definition that also captures the OH vibrational frequency downshifts well: R(H⋯O) ≤2.5 Å and the O-H⋯O angle θ ≥ 100°; and (iv) the frequency vs electric field relations for water and hydroxides (i.e., dissociated water) follow different and well-separated curves.

H2O2(s) and H2O2·2H2O(s) crystals compared with ices: DFT functional assessment and D3 analysis.

The H2O and H2O2 molecules resemble each other in a multitude of ways as has been noted in the literature. Here, we present density functional theory (DFT) calculations for the H2O2(s) and H2O2·2H2O(s) crystals and make selected comparisons with ice polymorphs. The performance of a number of dispersion-corrected density functionals-both self-consistent and a posteriori ones-are assessed, and we give special attention to the D3 correction and its effects. The D3 correction to the lattice energies is large: for H2O2(s) the D3 correction constitutes about 25% of the lattice energy using PBE, much more for RPBE, much less for SCAN, and it primarily arises from non-H-bonded interactions out to about 5 Å.The large D3 corrections to the lattice energies are likely a consequence of several effects: correction for missing dispersion interaction, the ability of D3 to capture and correct various other kinds of limitations built into the underlying DFT functionals, and finally some degree of cell-contraction-induced polarization enhancement. We find that the overall best-performing functionals of the twelve examined are optPBEvdW and RPBE-D3. Comparisons with DFT assessments for ices in the literature show that where the same methods have been used, the assessments largely agree.

Phonon-Assisted Hot Carrier Generation in Plasmonic Semiconductor Systems.

Plasmonic materials have optical cross sections that exceed by 10-fold their geometric sizes, making them uniquely suitable to convert light into electrical charges. Harvesting plasmon-generated hot carriers is of interest for the broad fields of photovoltaics and photocatalysis; however, their direct utilization is limited by their ultrafast thermalization in metals. To prolong the lifetime of hot carriers, one can place acceptor materials, such as semiconductors, in direct contact with the plasmonic system. Herein, we report the effect of operating temperature on hot electron generation and transfer to a suitable semiconductor. We found that an increase in the operation temperature improves hot electron harvesting in a plasmonic semiconductor hybrid system, contrasting what is observed on photodriven processes in nonplasmonic systems. The effect appears to be related to an enhancement in hot carrier generation due to phonon coupling. This discovery provides a new strategy for optimization of photodriven energy production and chemical synthesis.

Photoinduced Adsorption and Oxidation of SO2 on Anatase TiO2(101).

Adsorption of molecules is a fundamental step in all heterogeneous catalytic reactions. Nevertheless, the basic mechanism by which photon-mediated adsorption processes occur on solid surfaces is poorly understood, mainly because they involve excited catalyst states that complicate the analysis. Here we demonstrate a method by which density functional theory (DFT) can be used to quantify photoinduced adsorption processes on transition metal oxides and reveal the fundamental nature of these reactions. Specifically, the photoadsorption of SO2 on TiO2(101) has been investigated by using a combination of DFT and in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). The combined experimental and theoretical approach gives a detailed description of the photocatalytic desulfurization process on TiO2, in which sulfate forms as a stable surface product that is known to poison the catalytic surface. This work identifies surface-SO42- as the sulfate species responsible for the surface poisoning and shows how this product can be obtained from a stepwise oxidation of SO2 on TiO2(101). Initially, the molecule binds to a lattice O2- ion through a photomediated adsorption process and forms surface sulfite, which is subsequently oxidized into surface-SO42- by transfer of a neutral oxygen from an adsorbed O2 molecule. The work further explains how the infrared spectra associated with this oxidation product change during interactions with water and surface hydroxyl groups, which can be used as fingerprints for the surface reactions. The approach outlined here can be generalized to other photo- and electrocatalytic transition metal oxide systems.

Fast, Low-Cost Synthesis of ZnO:Eu Nanosponges and the Nature of Ln Doping in ZnO.

A low-cost template-free solution chemical route to highly porous nanocrystalline sponges of ZnO-EuO1.5 with 0-5 mol % Eu is presented. The process uses Zn- and Eu-acetate-nitrate and triethanolamine as precursors in methanol. After evaporation of the solvent and heating at 200 °C for 3 min, crystalline ZnO:Eu sponges with minor amounts of organic residues were obtained. Heating to 400 °C replaced the organics with carbonate, which in its turn was decomposed at temperatures below 600 °C, forming ZnO:Eu sponges. Samples heated to 200-1000 °C for 3 min were studied with XRD, SEM, TEM, TG, XPS, and IR spectroscopy. The ZnO:Eu crystallite sizes could be tuned from below 10 nm for sponges prepared at 200-500 °C, to over 100 nm range at 900 °C, without sintering of the overall microstructure. XRD showed the presence of hexagonal ZnO:Eu (or at 700-1000 °C, ZnO:Eu and cubic Eu2O3) as the only phases present. The ZnO:Eu had slightly larger unit cell dimensions than the literature value of ZnO for samples obtained at 200-600 °C, while the unit cells of samples obtained at higher temperatures were quite close to the value of undoped ZnO. XPS showed that Eu was mainly in its 3+ state and well-distributed within the sponges but segregated at the ZnO sponge surface upon heating at 700-1000 °C, in accordance with XRD studies showing Eu2O3 formation.

The water/ceria(111) interface: Computational overview and new structures.

Thin film structures of water on the CeO2(111) surface for coverages between 0.5 and 2.0 water monolayers have been optimized and analyzed using density functional theory (optPBE-vdW functional). We present a new 1.0 ML structure that is both the lowest in energy published and features a hydrogen-bond network extending the surface in one-dimension, contrary to what has been found in the literature, and contrary to what has been expected due to the large bulk ceria cell dimension. The adsorption energies for the monolayer and multilayered water structures agree well with experimental temperature programmed desorption results from the literature, and we discuss the stability window of CeO2(111) surfaces covered with 0.5-2.0 ML of water.

Anion-mediated electronic effects in reducible oxides: Tuning the valence band of ceria via fluorine doping.

Combining experimental spectroscopy and hybrid density functional theory calculations, we show that the incorporation of fluoride ions into a prototypical reducible oxide surface, namely, ceria(111), can induce a variety of nontrivial changes to the local electronic structure, beyond the expected increase in the number of Ce3+ ions. Our resonant photoemission spectroscopy results reveal new states above, within, and below the valence band, which are unique to the presence of fluoride ions at the surface. With the help of hybrid density functional calculations, we show that the different states arise from fluoride ions in different atomic layers in the near surface region. In particular, we identify a structure in which a fluoride ion substitutes for an oxygen ion at the surface, with a second fluoride ion on top of a surface Ce4+ ion giving rise to F 2p states which overlap the top of the O 2p band. The nature of this adsorbate F--Ce4+ resonant enhancement feature suggests that this bond is at least partially covalent. Our results demonstrate the versatility of anion doping as a potential means of tuning the valence band electronic structure of ceria.

Maximally resolved anharmonic OH vibrational spectrum of the water/ZnO(101¯0) interface from a high-dimensional neural network potential.

Unraveling the atomistic details of solid/liquid interfaces, e.g., by means of vibrational spectroscopy, is of vital importance in numerous applications, from electrochemistry to heterogeneous catalysis. Water-oxide interfaces represent a formidable challenge because a large variety of molecular and dissociated water species are present at the surface. Here, we present a comprehensive theoretical analysis of the anharmonic OH stretching vibrations at the water/ZnO(101¯0) interface as a prototypical case. Molecular dynamics simulations employing a reactive high-dimensional neural network potential based on density functional theory calculations have been used to sample the interfacial structures. In the second step, one-dimensional potential energy curves have been generated for a large number of configurations to solve the nuclear Schrödinger equation. We find that (i) the ZnO surface gives rise to OH frequency shifts up to a distance of about 4 Å from the surface; (ii) the spectrum contains a number of overlapping signals arising from different chemical species, with the frequencies decreasing in the order ν(adsorbed hydroxide) > ν(non-adsorbed water) > ν(surface hydroxide) > ν(adsorbed water); (iii) stretching frequencies are strongly influenced by the hydrogen bond pattern of these interfacial species. Finally, we have been able to identify substantial correlations between the stretching frequencies and hydrogen bond lengths for all species.

Fluorine impurities at CeO2(111): Effects on oxygen vacancy formation, molecular adsorption, and surface re-oxidation.

We investigate the effects of anion doping with fluorine impurities on the chemistry of the CeO2 (111) facet, using the results of DFT + U calculations. We consider three prototypical processes: the formation of oxygen vacancies, the adsorption of O2 and H2O molecules, and the re-oxidation of the surface with fragments of the two molecules. We find that the first two of these processes are not strongly affected, but that the presence of F lowers the energy gained in the re-oxidation of the surface in comparison to the healing of an oxygen vacancy, by 1.47 eV in the case of O2 (provided that the F is part of a cluster) and by 0.92 eV in the case of H2O. Based on these results, we suggest that F could enhance the redox chemistry of ceria by toggling between being in the surface and on the surface, effectively facilitating the release of lattice O by acting as a "place holder" for it. Finally, we find that the desorption of F as either 12F2 or HF is energetically unfavourable, suggesting that F doped ceria should be stable in the presence of O2 and H2O.

Many competing ceria (110) oxygen vacancy structures: from small to large supercells.

We present periodic "DFT+U" studies of single oxygen vacancies on the CeO(2)(110) surface using a number of different supercells, finding a range of different local minimum structures for the vacancy and its two accompanying Ce(III) ions. We find three different geometrical structures in combination with a variety of different Ce(III) localization patterns, several of which have not been studied before. The desired trapping of electrons was achieved in a two-stage optimization procedure. We find that the surface oxygen nearest to the vacancy either moves within the plane towards the vacancy, or rises out of the surface into either a symmetric or an unsymmetric bridge structure. Results are shown in seven slab geometry supercells, p(2 × 1), p(2 × 2), p(2 × 3), p(3 × 2), p(2 × 4), p(4 × 2), and p(3 × 3), and indicate that the choice of supercell can affect the results qualitatively and quantitatively. An unsymmetric bridge structure with one nearest and one next-nearest neighbour Ce(III) ion (a combination of localizations not previously found) is the ground state in all (but one) of the supercells studied here, and the relative stability of other structures depends strongly on supercell size. Within any one supercell the formation energies of the different vacancy structures differ by up to 0.5 eV, but the same structure can vary by up to ∼1 eV between supercells. Furthermore, finite size scaling suggests that the remaining errors (compared to still larger supercells) can also be ∼1 eV for some vacancy structures.

Predicting Frequency from the External Chemical Environment: OH Vibrations on Hydrated and Hydroxylated Surfaces.

Robust correlation curves are essential to decipher structural information from IR-vibrational spectra. However, for surface-adsorbed water and hydroxides, few such correlations have been presented in the literature. In this paper, OH vibrational frequencies are correlated against 12 structural descriptors representing the quantum mechanical or geometrical environment, focusing on those external to the vibrating molecule. A nonbiased fitting procedure based on Gaussian process regression (GPR) was used alongside simple analytical functional forms. The training data consist of 217 structurally unique OH groups from 38 water/metal oxide interface systems for MgO, CaO and CeO2, all optimized at the DFT level, and the fully anharmonic and uncoupled OH vibrational signatures were calculated. Among our results, we find the following: (i) The intermolecular R(H···O) hydrogen bond distance is particularly strong, indicating the primary cause of the frequency shift. (ii) Similarly, the electric field along the H-bond vector is also a good descriptor. (iii) Highly detailed machine learning descriptors (ACSF, SOAP) are less intuitive but were found to be more capable descriptors. (iv) Combinations of geometric and QM descriptors give the best predictions, supplying complementary information.

SOx on ceria from adsorbed SO2.

Results from first-principles calculations present a rather clear picture of the interaction of SO(2) with unreduced and partially reduced (111) and (110) surfaces of ceria. The Ce(3+)∕Ce(4+) redox couple, together with many oxidation states of S, give rise to a multitude of SO(x) species, with oxidation states from +III to +VI. SO(2) adsorbs either as a molecule or attaches via its S-atom to one or two surface oxygens to form sulfite (SO(3)(2-)) and sulfate (SO(4)(2-)) species, forming new S-O bonds but never any S-Ce bonds. Molecular adsorption is found on the (111) surface. SO(3)(2-) structures are found on both the (111) and (110) surfaces of both stoichiometric and partially reduced ceria. SO(4)(2-) structures are observed on the (110) surface together with the formation of two reduced Ce(3+) surface cations. SO(2) can also partially heal the ceria oxygen vacancies by weakening a S-O bond, when significant electron transfer from the surface (Ce4f) into the lowest unoccupied molecular orbital of the SO(2) adsorbate takes place and oxidizes the surface Ce(3+) cations. Furthermore, we propose a mechanism that could lead to monodentate sulfate formation at the (111) surface.

Using DFTB to Model Photocatalytic Anatase-Rutile TiO2 Nanocrystalline Interfaces and Their Band Alignment.

Band alignment effects of anatase and rutile nanocrystals in TiO2 powders lead to electron-hole separation, increasing the photocatalytic efficiency of these powders. While size effects and types of possible alignments have been extensively studied, the effect of interface geometries of bonded nanocrystal structures on the alignment is poorly understood. To allow conclusive studies of a vast variety of bonded systems in different orientations, we have developed a new density functional tight-binding parameter set to properly describe quantum confinement in nanocrystals. By applying this set, we found a quantitative influence of the interface structure on the band alignment.

Curvature Constrained Splines for DFTB Repulsive Potential Parametrization.

The Curvature Constrained Splines (CCS) methodology has been used for fitting repulsive potentials to be used in SCC-DFTB calculations. The benefit of using CCS is that the actual fitting of the repulsive potential is performed through quadratic programming on a convex objective function. This guarantees a unique (for strictly convex) and optimum two-body repulsive potential in a single shot, thereby making the parametrization process robust, and with minimal human effort. Furthermore, the constraints in CCS give the user control to tune the shape of the repulsive potential based on prior knowledge about the system in question. Herein, we developed the method further with new constraints and the capability to handle sparse data. We used the method to generate accurate repulsive potentials for bulk Si polymorphs and demonstrate that for a given Slater-Koster table, which reproduces the experimental band structure for bulk Si in its ground state, we are unable to find one single two-body repulsive potential that can accurately describe the various bulk polymorphs of silicon in our training set. We further demonstrate that to increase transferability, the repulsive potential needs to be adjusted to account for changes in the chemical environment, here expressed in the form of a coordination number. By training a near-sighted Atomistic Neural Network potential, which includes many-body effects but still essentially within the first-neighbor shell, we can obtain full transferability for SCC-DFTB in terms of describing the energetics of different Si polymorphs.

B3LYP calculations of cerium oxides.

In this paper we evaluate the performance of density functional theory with the B3LYP functional for calculations on ceria (CeO(2)) and cerium sesquioxide (Ce(2)O(3)). We demonstrate that B3LYP is able to describe CeO(2) and Ce(2)O(3) reasonably well. When compared to other functionals, B3LYP performs slightly better than the hybrid functional PBE0 for the electronic properties but slightly worse for the structural properties, although neither performs as well as LDA+U(U=6 eV) or PBE+U(U=5 eV). We also make an extensive comparison of atomic basis sets suitable for periodic calculations of these cerium oxides. Here we conclude that there is currently only one type of cerium basis set available in the literature that is able to give a reasonable description of the electronic structure of both CeO(2) and Ce(2)O(3). These basis sets are based on a 28 electron effective core potential (ECP) and 30 electrons are attributed to the valence space of cerium. Basis sets based on 46 electron ECPs fail for these materials.

Density Functional Tight Binding Theory Approach for the CO2 Reduction Reaction Paths on Anatase TiO2 Surfaces.

Herein, we have investigated the CO2 reduction paths on the (101) anatase TiO2 surface using an approach based on the density functional tight binding (DFTB) theory. We analyzed the reaction paths for the conversion of carbon dioxide to methane by performing a large number of calculations with intermediates placed in various orientations and locations at the surface. Our results show that the least stable intermediate is CO2H and therefore a key bottleneck is the reduction of CO2 to formic acid. Hydrogen adsorption is also weak and would also be a limiting factor, unless very high pressures of hydrogen are used. The results from our DFTB approach are in good agreement with the hybrid functional based density functional theory calculations presented in the literature.

Oligomer Electrolytes for Light-Emitting Electrochemical Cells: Influence of the End Groups on Ion Coordination, Ion Binding, and Turn-on Kinetics.

The electrolyte is an essential constituent of the light-emitting electrochemical cell (LEC), since its operating mechanism is dependent on the redistribution of mobile ions in the active layer. Recent developments of new ion transporters have yielded high-performance devices, but knowledge about the interactions between the ionic species and the ion transporters and the influence of these interactions on the LEC performance is lacking. We therefore present a combined computational and experimental effort that demonstrates that the selection of the end group in a star-branched oligomeric ion transporter based on trimethylolpropane ethoxylate has a paramount influence on the ionic interactions in the electrolyte and thereby also on the performance of the corresponding LECs. With hydroxyl end groups, the cation from the salt is strongly coordinated to the ion transporter, which leads to suppression of ion pairing, but the penalty is a hindered ion release and a slow turn-on for the LEC devices. With methoxy end groups, an intermediate coordination strength is seen together with the formation of contact ion pairs, but the LEC performance is very good with fast turn-on. Using a series of ion transporters with alkyl carbonate end groups, the ion transporter:cation coordination strength is lowered further, but the turn-on kinetics are slower than what is seen for devices comprising the methoxy end-capped ion transporter.

STM Images of Anionic Defects at CeO2(111)-A Theoretical Perspective.

We present a theoretically oriented analysis of the appearance and properties of plausible candidates for the anionic defects observed in scanning tunneling microscopy (STM) experiments on CeO2(111). The simulations are based on density functional theory (DFT) and cover oxygen vacancies, fluorine impurities and hydroxyl groups in the surface and sub-surface layers. In the surface layer, all three appear as missing spots in the oxygen sublattice in filled state simulated STM images, but they are distinguishable in empty state images, where surface oxygen vacancies and hydroxyls appear as, respectively, diffuse and sharp bright features at oxygen sites, while fluorine defects appear as triangles of darkened Ce ions. In the sub-surface layer, all three defects present more complex patterns, with different combinations of brightened oxygen ion triangles and/or darkened Ce ion triangles, so we provide image maps to support experimental identification. We also discuss other properties that could be used to distinguish the defects, namely their diffusion rates and distributions.