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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Multiscale Modeling of Agglomerated Ceria Nanoparticles: Interface Stability and Oxygen Vacancy Formation.

The interface formation and its effect on redox processes in agglomerated ceria nanoparticles (NPs) have been investigated using a multiscale simulation approach with standard density functional theory (DFT), the self-consistent-charge density functional tight binding (SCC-DFTB) method, and a DFT-parameterized reactive force-field (ReaxFF). In particular, we have modeled Ce40O80 NP pairs, using SCC-DFTB and DFT, and longer chains and networks formed by Ce40O80 or Ce132O264 NPs, using ReaxFF molecular dynamics simulations. We find that the most stable {111}/{111} interface structure is coherent whereas the stable {100}/{100} structures can be either coherent or incoherent. The formation of {111}/{111} interfaces is found to have only a very small effect on the oxygen vacancy formation energy, E vac. The opposite holds true for {100}/{100} interfaces, which exhibit significantly lower E vac values than the bare surfaces, despite the fact that the interface formation eliminates reactive {100} facets. Our results pave the way for an increased understanding of ceria NP agglomeration.

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.

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.

Dedoping of Lead Halide Perovskites Incorporating Monovalent Cations.

We report significant improvements in the optoelectronic properties of lead halide perovskites with the addition of monovalent ions with ionic radii close to Pb2+. We investigate the chemical distribution and electronic structure of solution processed CH3NH3PbI3 perovskite structures containing Na+, Cu+, and Ag+, which are lower valence metal ions than Pb2+ but have similar ionic radii. Synchrotron X-ray diffraction reveals a pronounced shift in the main perovskite peaks for the monovalent cation-based films, suggesting incorporation of these cations into the perovskite lattice as well as a preferential crystal growth in Ag+ containing perovskite structures. Furthermore, the synchrotron X-ray photoelectron measurements show a significant change in the valence band position for Cu- and Ag-doped films, although the perovskite bandgap remains the same, indicating a shift in the Fermi level position toward the middle of the bandgap. Such a shift infers that incorporation of these monovalent cations dedope the n-type perovskite films when formed without added cations. This dedoping effect leads to cleaner bandgaps as reflected by the lower energetic disorder in the monovalent cation-doped perovskite thin films as compared to pristine films. We also find that in contrast to Ag+ and Cu+, Na+ locates mainly at the grain boundaries and surfaces. Our theoretical calculations confirm the observed shifts in X-ray diffraction peaks and Fermi level as well as absence of intrabandgap states upon energetically favorable doping of perovskite lattice by the monovalent cations. We also model a significant change in the local structure, chemical bonding of metal-halide, and the electronic structure in the doped perovskites. In summary, our work highlights the local chemistry and influence of monovalent cation dopants on crystallization and the electronic structure in the doped perovskite thin films.

Photoinduced Stark Effects and Mechanism of Ion Displacement in Perovskite Solar Cell Materials.

Organometallic halide perovskites (OMHPs) have recently emerged as a promising class of materials in photovoltaic technology. Here, we present an in-depth investigation of the physics in these systems by measuring the photoinduced absorption (PIA) in OMHPs as a function of materials composition, excitation wavelength, and modulation frequency. We report a photoinduced Stark effect that depends on the excitation wavelength and on the dipole strength of the monovalent cations in the A position of the ABX3 perovskite. The results presented are corroborated by density functional theory calculations and provide fundamental information about the photoinduced local electric field change under blue and red excitation as well as insights into the mechanism of light-induced ion displacement in OMHPs. For optimized perovskite solar cell devices beyond 19% efficiency, we show that excess thermalization energy of blue photons plays a role in overcoming the activation energy for ion diffusion.

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.