Search

Iron (Fe)-doped β-nickel oxyhydroxide (β-NiOOH) is a highly active, noble-metal-free electrocatalyst for the oxygen evolution reaction (OER), with the latter being the bottleneck in electrochemical water splitting for sustainable hydrogen production. The mechanisms underlying how the Fe dopant modulates this host material's water electro-oxidation activity are still not entirely clear. Here, we combine hybrid density functional theory (DFT) and Hubbard-corrected DFT to investigate the OER activity of the most thermodynamically favorable (and therefore, expected to be the majority) crystallographic facets of β-NiOOH, namely (0001) and (101̄0). By considering active sites involving both oxidation and reduction of the transition-metal active center during the redox cycle on these two different facets, we show that six-fold-lattice-coordinated Fe in β-NiOOH is redox inactive towards both oxidation and reduction while five-fold-lattice-coordinated Fe in β-NiOOH does exhibit redox activity. However, the determined redox activity of Fe (or lack of it) is not indicative of good (or bad) performance as a dopant on these two facets. Three of the four active sites investigated (oxo and hydroxo sites on (0001) and a hydrated site on (101̄0)) exhibit only a marginal (<0.1 V) decrease or increase in the thermodynamic overpotential upon doping with Fe. Only one of the redox-active sites investigated, the hydroxo site on (101̄0), exhibits a large attenuation in the thermodynamic overpotential upon doping (to ∼0.52 V from 0.86 V), although the doped overpotential is larger than that observed experimentally for Fe-doped NiOOH. Thus, although pure β-NiOOH facets containing four-, five-, or six-fold lattice-coordinated Ni sites have roughly equal OER activities, yielding similar OER onset potentials (shown in A. Govind Rajan, J. M. P. Martirez and E. A. Carter, J. Am. Chem. Soc. , 2020, 142 , 3600-3612), only those facets containing four-fold lattice-coordinated Fe ( e.g. , as shown in J. M. P. Martirez and E. A. Carter, J. Am. Chem. Soc. , 2019, 141 , 693-705) would be active under analogous conditions for the Fe-doped material. It follows that, while undoped β-NiOOH demonstrates a roughly facet-independent oxygen evolution activity, the activity of Fe-doped β-NiOOH strongly depends on the crystallographic facet. Our study further motivates the investigation of strategies for the selective growth of facets with low iron coordination number to enhance the water splitting activity of Fe-doped β-NiOOH.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

We benchmark calculated interlayer spacings, average topotactic voltages, thermodynamic stabilities, and band gaps in layered lithium transition-metal oxides (TMOs) and their de-lithiated counterparts, which are used in lithium-ion batteries as positive electrode materials, against available experimental data. Specifically, we examine the accuracy of properties calculated within density functional theory (DFT) using eight different treatments of electron exchange-correlation: the strongly constrained and appropriately normed (SCAN) and Perdew-Burke-Ernzerhof (PBE) density functionals, Hubbard- U -corrected SCAN and PBE ( i.e. , SCAN+ U and PBE+ U ), and SCAN(+ U ) and PBE(+ U ) with added long-range dispersion (D) interactions ( i.e. , DFT(+ U )+D). van der Waals interactions are included respectively via the revised Vydrov-Van Voorhis (rVV10) for SCAN(+ U ) and the DFT-D3 for PBE(+ U ). We find that SCAN-based functionals predict larger voltages due to an underestimation of stability of the MO 2 systems, while also predicting smaller interlayer spacings compared to their PBE-based counterparts. Furthermore, adding dispersion corrections to PBE has a greater effect on voltage predictions and interlayer spacings than with SCAN, indicating that DFT-SCAN - despite being a ground-state theory - fortuitously captures some short and medium-range dispersion interactions better than PBE. While SCAN-based and PBE-based functionals yield qualitatively similar band gap predictions, there is no significant quantitative improvement of SCAN-based functionals over the corresponding PBE-based versions. Finally, we expect SCAN-based functionals to yield more accurate property predictions than the respective PBE-based functionals for most TMOs, given SCAN's stronger theoretical underpinning and better predictions of systematic trends in interlayer spacings, intercalation voltages, and band gaps obtained in this work.

Iron-doped nickel oxyhydroxide has been identified as one of the most active alkaline oxygen evolution reaction (OER) catalysts, exhibiting an overpotential lower than values observed for state-of-the-art precious metal catalysts. Several computational investigations have found widely varying effects of doping on the theoretical overpotential of the OER on NiOx. Comparisons of these results are made difficult by the numerous differences in the structural and computational parameters used in these studies. In this work, within a consistent framework, we calculate the theoretical overpotentials for reactions occurring on the most stable, basal plane of undoped and doped β-NiOOH. We compare the activities of Fe(iii), Co(iii), and Mn(iii) doping using density functional theory with Hubbard-like U corrections on the transition-metal d orbitals. We compare the effect of surface and subsurface doping in order to establish whether the dopants act as new active sites for the reaction or whether they induce more widespread changes in the material. The results of our study find only a small reduction in the overpotential (∼0.1 and ≤0.05 V when doped in the surface and subsurface layers, respectively) for the three dopants, if doped in the dominant basal plane. This is much less than the reductions of 0.3 V experimentally observed for the most active Fe-doped systems. Furthermore, the magnitudes of reductions in overpotentials for the three dopants are similar. This work therefore disqualifies the possibility of enhancing the activity of the dominant exposed basal plane of β-NiOOH through substitutional doping.

Adsorbed protons that develop hydride character have been proposed to play a role in the mechanism of CO 2 reduction catalyzed by pyridine on GaP photoelectrodes. Investigating their stability represents an important step towards vetting this mechanism. In this contribution, the relative stability of the adsorbed protons is determined using cluster models with dispersion-corrected density functional theory and continuum solvation. Proton acidity constants computed under typical experimental conditions are compared to the acidity constants of other relevant species. The adsorbed protons are predicted to be very stable, suggesting that they will be present on the surface and available to be reduced to surface hydrides that could possibly react with adsorbed pyridine to form adsorbed dihydropyridine, a previously proposed co-catalyst. However, the high stability of such protons also suggests that the surface does not represent a significant proton source; as a consequence, protons required in the proposed mechanism must be provided by a different source such as the acidified aqueous solution in contact with the electrode surface.

Fourier transform infrared (FTIR) microspectroscopy and confocal imaging have been used to demonstrate that the neutral rhenium(i) tricarbonyl 1,10-phenanthroline complex bound to 4-cyanophenyltetrazolate as the ancillary ligand is able to localise in regions with high concentrations of polar lipids such as phosphatidylethanolamine (PE), sphingomyelin, sphingosphine and lysophosphatidic acid (LPA) in mammalian adipocytes.

Reducing operating temperatures is a key step in making solid oxide fuel cell (SOFC) technology viable. A promising strategy for accomplishing this goal is employing mixed ion-electron conducting (MIEC) cathodes. La1-xSrxCo1-yFeyO3-δ (LSCF) is the most widely employed MIEC cathode material; however, rational optimization of the composition of LSCF requires fundamental insight linking its electronic structure to its defect chemistry. To provide the necessary insight, density functional theory plus U (DFT+U) calculations are used to investigate the electronic structure of LSCF (xSr = 0.50, yCo = 0.25). The DFT+U calculations show that LSCF has a significantly different electronic structure than La1-xSrxFeO3 because of the addition of cobalt, but that minimal electronic structure differences exist between La0.5Sr0.5Co0.25Fe0.75O3 and La0.5Sr0.5Co0.5Fe0.5O3. The oxygen vacancy formation energy (ΔEf,vac) is calculated for residing in different local environments within La0.5Sr0.5Co0.25Fe0.75O3. These results show that configurations have the highest ΔEf,vac, while have the lowest ΔEf,vac and may act as traps for . We conclude that compositions with more Fe than Co are preferred because the additional sites would lead to higher overall ΔEf,vac (and lower concentrations), while the trapping strength of the sites is relatively weak (∼0.3 eV).

The kinetics of hydrogen abstraction by five radicals (H, CH3, O((3)P), OH, and HO2) from a biodiesel surrogate, methyl propanoate (MP), is theoretically investigated. We employ high-level ab initio quantum chemistry methods, coupled-cluster singles and doubles with perturbative triples correction (CCSD(T)) and multi-reference singles and doubles configuration interaction (MRSDCI) with the Davidson-Silver (DS) correction, and obtain chemically accurate reaction energetics. Overall, MRSDCI + DS predicts comparable energetics to CCSD(T) for MP + H/CH3/O/OH. The rate constants are computed using transition state theory (TST-Rice-Ramsperger-Kassel-Marcus theory) in conjunction with the separable-hindered-rotor approximation, variable reaction coordinate TST, and the multi-structure all-structure (MS-AS) approach. A simplified method, semi-multi-structure, is also employed for MP + OH/HO2, and the rate coefficients with this less expensive method are in good agreement with the results obtained with the MS-AS method. The fitted modified Arrhenius expressions are provided over a temperature range of 250 to 2000 K. The predicted rate coefficients for MP + OH agree remarkably well with experimental data over a wide temperature range. Branching ratio analysis of all the studied reactions shows that abstractions of the secondary H atoms within MP are expected to dominate the consumption of fuel at low temperatures, and the contributions of abstractions from the two methyl groups increase with temperature for all abstracting radicals.

The unimolecular dissociation and isomerization kinetics of the three methyl propanoate (MP) radicals, CH3CH2C(=O)OĊH2 (MP-m), CH3ĊHC(=O)OCH3 (MP-α), and ĊH2CH2C(=O)OCH3 (MP-β), are theoretically investigated using high-level ab initio methods and the Rice-Ramsperger-Kassel-Marcus (RRKM)/master equation (ME) theory. Stationary-point energies are obtained using the coupled cluster singles and doubles with perturbative triples correction (CCSD(T)), multi-reference singles and doubles configuration interaction (MRSDCI) with the Davidson-Silver (DS) correction, and multi-reference averaged coupled pair functional (MRACPF2) theories. The isomerization barriers between the three radicals are predicted to be generally lower than the corresponding bond dissociation channels, leading to a strongly coupled reaction system in subsequent kinetics studies. The phenomenological temperature- and pressure-dependent rate coefficients are computed using the RRKM/ME theory over a temperature range of 500 to 2000 K and at a pressure range of 0.01 atm to the high-pressure limit, which are then fitted to modified Arrhenius expressions. The β-scission rate coefficients of MP-α to CH3CHC(=O) and CH3O[combining dot above] are predicted to be the smallest because of its highest activation energy among all studied unimolecular reactions channels. Analysis of branching fractions shows that both MP-m and MP-α radicals mainly decompose to the bimolecular products CH3CH2Ċ(=O) and H2CO, whereas the MP-β radical primarily decomposes via cleavage of a C-C bond to form C2H4 and CH3OĊ(=O). The isomerization channels dominate at low temperatures, the branching fractions of which decrease with increasing temperature and become very minor at about 2000 K. Our accurate rate coefficients and branching fractions help to illuminate the unique combustion properties of MP.

Alloying nickel oxide (NiO) with lithium oxide (Li2O) at high Li concentrations may reduce NiO's band gap and expand its use as a light absorber in photocatalytic and tandem dye-sensitized solar cell technologies. In this work, we evaluate the viability of this alloy as a p-type hole transport material. We use embedded cluster models, along with unrestricted Hartree-Fock and complete active space self-consistent field theories, to study the impact of alloying on polaronic transport of holes. Our calculated energy barrier for hole transfer in undoped NiO is in excellent agreement with the experimental value of ∼0.1 eV. We predict that hole transport in NiO is anisotropic and mostly confined parallel to the (111) ferromagnetic planes. Applying the same model to Li-doped NiO indicates that isolated Li ions do not introduce free holes into NiO samples. However, free holes can be created in the homogeneous Li0.125Ni0.875O alloy, in which the Li concentration is very high. Our kinetic Monte Carlo calculations show that hole mobility in this alloy is lower than in undoped NiO. However, the additional free holes and the predicted lower band gap of Li0.125Ni0.875O should increase hole conductivity compared to NiO upon alloy formation. Therefore, Li0.125Ni0.875O alloys have potential for use as a hole transporter, as well as a sunlight absorber, in a variety of solar energy applications.

We use ab initio density functional theory + U calculations to characterize the oxide ion diffusion process in bulk Sr2Fe1.5Mo0.5O(6-δ) (SFMO) by analyzing the formation and migration of oxygen vacancies. We show that SFMO's remarkable ionic conductivity arises from its intrinsic content of oxygen vacancies and a predicted very low migration barrier of such vacancies. Theoretical analysis of the electronic structure reveals a crucial role played by strongly hybridized Fe 3d/O 2p states to achieve the attendant mixed ion-electron conductor character so important for intermediate temperature fuel cell operation. We predict a next-nearest-neighbor-type migration pathway for the O(2-) ion should dominate. The low energy barrier of this pathway is mainly related to electrostatic interactions with homogeneously distributed Mo in the SFMO sublattice. We identify the reasons why Fe-rich perovskites, with the key addition of a certain concentration of Mo, produce excellent electronic and ionic transport properties so crucial for efficient operation of intermediate temperature solid oxide fuel cells.

Multi-reference Configuration Interaction (MRCI) is often used to predict the electronic structures and reaction energetics of small molecules with very high accuracy. Unfortunately, MRCI is inapplicable to large or even medium-sized molecules for two reasons: its computational cost scales poorly with molecule size and MRCI methods are not size extensive, leading to large energy errors. We have developed a new local (L) and approximately size extensive MRCI method that addresses both shortcomings. Truncating long-range electron correlation in a local orbital basis as well as efficient processing of two-electron integrals via Cholesky decomposition (CD) and integral screening reduce the computational cost to O(N(3)) with a small prefactor. A priori and a posteriori size extensivity corrections are both considered. The Multi-reference Averaged Coupled-Pair Functional (MRACPF) provides approximate size extensivity by modifying the energy functional. The very inexpensive Davidson-Silver and Pople a posteriori schemes also produce quite accurate corrections over a large range of molecular size. Employing CD-LMRACPF is slightly more expensive than using a Davidson-type correction, but the former gives superior results. Molecules with up to 50 heavy atoms can be treated with our CD-LMRACPF method thus far.

The positions of electronic band edges are one important metric for determining a material's capability to function in a solar energy conversion device that produces fuels from sunlight. In particular, the position of the valence band maximum (conduction band minimum) must lie lower (higher) in energy than the oxidation (reduction) reaction free energy in order for these reactions to be thermodynamically favorable. We present first principles quantum mechanics calculations of the band edge positions in five transition metal oxides and discuss the feasibility of using these materials in photoelectrochemical cells that produce fuels, including hydrogen, methane, methanol, and formic acid. The band gap center is determined within the framework of DFT+U theory. The valence band maximum (conduction band minimum) is found by subtracting (adding) half of the quasiparticle gap obtained from a non-self-consistent GW calculation. The calculations are validated against experimental data where possible; results for several materials including manganese(ii) oxide, iron(ii) oxide, iron(iii) oxide, copper(i) oxide and nickel(ii) oxide are presented.

Quantitative characterization of low-lying excited electronic states in materials is critical for the development of solar energy conversion materials. The many-body Green's function method known as the GW approximation (GWA) directly probes states corresponding to photoemission and inverse photoemission experiments, thereby determining the associated band structure. Several versions of the GW approximation with different levels of self-consistency exist in the field. While the GWA based on density functional theory (DFT) works well for conventional semiconductors, less is known about its reliability for strongly correlated semiconducting materials. Here we present a systematic study of the GWA using hematite (α-Fe(2)O(3)) as the benchmark material. We analyze its performance in terms of the calculated photoemission/inverse photoemission band gaps, densities of states, and dielectric functions. Overall, a non-self-consistent G(0)W(0) using input from DFT+U theory produces physical observables in best agreement with experiments., (This journal is © the Owner Societies 2011)

Understanding biochemical mechanisms and changes associated with disease conditions and, therefore, development of improved clinical treatments, is relying increasingly on various biochemical mapping and imaging techniques on tissue sections. However, it is essential to be able to ascertain whether the sampling used provides the full biochemical information relevant to the disease and is free from artefacts. A multi-modal micro-spectroscopic approach, including FTIR imaging and PIXE elemental mapping, has been used to study the molecular and elemental profile within cryofixed and formalin-fixed murine brain tissue sections. The results provide strong evidence that amino acids, carbohydrates, lipids, phosphates, proteins and ions, such as Cl(-) and K(+), leach from tissue sections into the aqueous fixative medium during formalin fixation of the sections. Large changes in the concentrations and distributions of most of these components are also observed by washing in PBS even for short periods. The most likely source of the chemical species lost during fixation is the extra-cellular and intra-cellular fluid of tissues. The results highlight that, at best, analysis of formalin-fixed tissues gives only part of the complete biochemical "picture" of a tissue sample. Further, this investigation has highlighted that significant lipid peroxidation/oxidation may occur during formalin fixation and that the use of standard histological fixation reagents can result in significant and differential metal contamination of different regions of tissue sections. While a consistent and reproducible fixation method may be suitable for diagnostic purposes, the findings of this study strongly question the use of formalin fixation prior to spectroscopic studies of the molecular and elemental composition of biological samples, if the primary purpose is mechanistic studies of disease pathogenesis.

Herein is described a general sampling protocol that includes culture, differentiation and fixing of cells in their preferred morphology on the one sample substrate (Si(3)N(4)) to enable subsequent diverse modern microspectroscopic analyses. The protocol enables unprecedented correlated and complementary information on the intracellular biochemistry of metabolic processes, diseases and their treatment, which offers the opportunity to revolutionize our understanding of cell and tissue biology at a molecular level. The culture of adherent cells onto inexpensive Si(3)N(4) membranes allows microspectroscopic analyses across the electromagnetic spectrum, from hard X-ray fluorescence (both XRF and XANES), through to visible and fluorescence light microscopies, and infrared microspectroscopy without substrate interference. Adherent mammalian cell lines (3T3-L1 adipocytes and H9c2 cardiac myocytes) illustrate the in vitro application of these protocols. The cells adhered strongly to Si(3)N(4) membranes and visually displayed normal proliferative and phenotypic growth; more importantly, rapid alcohol fixation of cells did not affect their structural integrity for subsequent analyses.

One obstacle in orbital-free density functional theory (OF-DFT) is the lack of accurate and transferable local pseudopotentials (LPSs). In this work, we build high quality LPSs by inverting Kohn-Sham (KS) equations on bulk valence electron densities to obtain an atom-centered local pseudopotential. With this approach, we build LPSs for Mg, Al, and Si, and then test them in KS DFT calculations of static bulk properties for several Mg, Al, and Si bulk structures as well as beta''-Al3Mg. Our Mg, Al, and Si LPSs produce correct ground state properties and phase orderings. These LPSs are then tested in KS-DFT calculations of surface energies for several low-index Mg and Al surfaces, point defect properties in hexagonal-close-packed (hcp) Mg, face-centered cubic (fcc) Al, and diamond Si, and stacking fault energies in fcc Al. All of these LPS results agree quantitatively with the results from nonlocal pseudopotentials with errors less than or equal to 40 meV per atom. Finally, we perform OF-DFT calculations for various Mg and Al structures, employing the Wang-Govind-Carter (WGC) nonlocal kinetic energy density functional (KEDF). The OF-DFT results generally agree well with the corresponding KS-DFT results. With our new Mg and Al LPSs and the WGC KEDF, OF-DFT now provides a practical method for accurate, large-scale first principles simulations of main group metals and their alloys.

A possible mechanism for shock-induced failure in aluminium involves atomic vacancies diffusing through the crystal lattice and agglomerating to form voids, which continue to grow, ultimately resulting in ductile fracture. We employ orbital-free density functional theory, a linear-scaling first-principles quantum mechanics method, to study vacancy formation, diffusion, and aggregation in aluminium under shock loading conditions of compression and tension. We calculate vacancy formation and migration energies, and find that while nearest-neighbor vacancy pairs are unstable, next-nearest-neighbor vacancy pairs are stable. As the number of nearby vacancies increases, we predict that vacancy clusters preferentially grow through next-nearest-neighbor vacancies. The energetics are found to be greatly affected by expansion and compression, leading to insight as to how vacancies behave under shock conditions.