Challenges in the use of density functional theory to examine catalysis by M‐doped ceria surfaces

For CeO₂ or M‐doped CeO₂ catalysts, reliable energetics associated with surface reactivity requires accurate representation of oxidized and reduced metal states. Density functional theory (DFT) is used extensively for metals and metal oxides; however, for strongly correlated electron materials, conventional DFT fails to predict both qualitative and quantitative properties. This is the result of a localized electron self‐interaction error that is inherit to DFT. DFT+U has shown promise in correcting energetic errors due to the self‐interaction error, however, its transferability across processes relevant to surface catalysis remains unclear. Hybrid functionals, such as HSE06, can also be used to correct this self‐interaction error. These hybrid functionals are computationally intensive, and especially demanding for periodic surface slab models. This perspective details the challenges in representing the energetics of M‐doped ceria catalyzed processes and examines using DFT extensions to model the localized electronic properties. © 2013 Wiley Periodicals, Inc.     

Systematic study of the effect of HSE functional internal parameters on the electronic structure and band gap of a representative set of metal oxides

The effect of the amount of Hartree–Fock mixing parameter (α ) and of the screening parameter (w ) defining the range separated HSE type hybrid functional is systematically studied for a series of seven metal oxides: TiO₂, ZrO₂, CuO₂, ZnO, MgO, SnO₂, and SrTiO₃. First, reliable band gap values were determined by comparing the optimal α reproducing the experiment with the inverse of the experimental dielectric constant. Then, the effect of the w in the HSE functional on the calculated band gap was explored in detail. Results evidence the existence of a virtually infinite number of combinations of the two parameters which are able to reproduce the experimental band gap, without a unique pair able to describe the full studied set of materials. Nevertheless, the results point out the possibility of describing the electronic structure of these materials through a functional including a screened HF exchange and an appropriate correlation contribution. © 2017 Wiley Periodicals, Inc.     

Accurate prediction of vertical electronic transitions of Ni(II) coordination compounds via time dependent density functional theory

Time dependent density functional theory calculations are completed for five Ni(II) complexes formed by polydentate peptides to predict the electronic absorption spectrum. The ligands examined were glycyl‐glycyl‐glycine (GGG), glycyl‐glycyl‐glycyl‐glycine (GGGG), glycyl‐glycyl‐histidine (GGH), glycyl‐glycyl‐cysteine (GGC), and triethylenetetramine (trien). Fifteen functionals and two basis sets were tested. On the basis of the mean absolute percent deviation (MAPD), the ranking among the functionals is: HSE06 ∼ MPW1PW91 ∼ PBE0 > ω‐B97x‐D ∼ B3P86 ∼ B3LYP ∼ CAM‐B3LYP > PBE ∼ BLYP ∼ BP86 > TPSS > TPSSh > BHandHLYP > M06 ≫ M06‐2X. Concerning the basis sets, the triple‐ζ def2‐TZVP performs better than the double‐ζ LANL2DZ. With the functional HSE06 and basis set def2‐TZVP the MAPD with respect to the experimental λ_max is 1.65% with a standard deviation of 1.26%. The absorption electronic spectra were interpreted in terms of vertical excitations between occupied and virtual MOs based on Ni‐d atomic orbitals. The electronic structure of the Ni(II) species is also discussed.     

Evaluation of range-separated hybrid density functionals for the prediction of vibrational frequencies, infrared intensities, and Raman activities

We present an assessment of different density functionals, with emphasis on range-separated hybrids, for the prediction of fundamental and harmonic vibrational frequencies, infrared intensities, and Raman activities. Additionally, we discuss the basis set convergence of vibrational properties of H2O with long-range corrected hybrids. Our results show that B3LYP is the best functional for predicting vibrational frequencies (both fundamental and harmonic); the screened-PBE hybrid (HSE) density functional works best for infrared intensities, and the long-range corrected PBE (LC-omegaPBE), M06-HF, and M06-L density functionals are almost as good as MP2 for predicting Raman activities. We show the predicted Raman spectrum of adenine as an example of a medium-size molecule where a DFT/Sadlej pVTZ calculation is affordable and compare our results against the experimental spectrum.     

Applicability of density functional theory in reproducing accurate vibrational spectra of surface bound species

The structural equilibrium parameters, the adsorption energies, and the vibrational frequencies of the nitrogen molecule and the hydrogen atom adsorbed on the (111) surface of rhodium have been investigated using different generalized‐gradient approximation (GGA), nonlocal correlation, meta‐GGA, and hybrid functionals, namely, Perdew, Burke, and Ernzerhof (PBE), Revised‐RPBE, vdW‐DF, Tao, Perdew, Staroverov, and Scuseria functional (TPSS), and Heyd, Scuseria, and Ernzerhof (HSE06) functional in the plane wave formalism. Among the five tested functionals, nonlocal vdW‐DF and meta‐GGA TPSS functionals are most successful in describing energetics of dinitrogen physisorption to the Rh(111) surface, while the PBE functional provides the correct chemisorption energy for the hydrogen atom. It was also found that TPSS functional produces the best vibrational spectra of the nitrogen molecule and the hydrogen atom on rhodium within the harmonic formalism with the error of ?2.62 and ?1.1% for the N? N stretching and Rh? H stretching frequency. Thus, TPSS functional was proposed as a method of choice for obtaining vibrational spectra of low weight adsorbates on metallic surfaces within the harmonic approximation. At the anharmonic level, by decoupling the Rh? H and N? N stretching modes from the bulk phonons and by solving one‐ and two‐dimensional Schrödinger equation associated with the Rh? H, Rh? N, and N? N potential energy we calculated the anharmonic correction for N? N and Rh? H stretching modes as ?31 cm^?1 and ?77 cm^?1 at PBE level. Anharmonic vibrational frequencies calculated with the use of the hybrid HSE06 function are in best agreement with available experiments. © 2014 Wiley Periodicals, Inc.     

A First‐Principles Study on the Electronic,Vibrational, and Thermodynamic Properties of Jadeite and its Tentative Low‐Density Polymorph

The high‐pressure clinopyroxene mineral jadeite (chemical composition NaAlSi₂O₆) was studied by density‐functional theory with respect to its electronic, vibrational, and thermodynamic properties, correctly reproducing the available experimental data. At a larger‐than‐normal volume, however, a low‐density alumosilicate phase with tetrahedral instead of octahedral Al–O coordination was identified. This low‐density phase was investigated theoretically, too, and the results were compared to jadeite and experimentally observed properties of what has been dubbed jadeite glass. It turned out that the theoretically obtained properties of this hypothetical polymorph, such as the bulk modulus, the molar volume and the vibrational frequencies, agree well with the corresponding properties of the glass phase. Hence, from an inverse structure‐property relationship we propose to model jadeite glass with the aforementioned low‐density alumosilicate phase and tetrahedral Al–O motifs, as suggested from first principles and corroborated from experiment. The possibilities as well as the limitations of electronic, phonon and thermodynamic properties calculations applied to such a polymorph are also discussed.     

Diamondoids approach to electronic,structural, and vibrational properties of GeSi superlattice nanocrystals: a first-principles study

Germanium silicide diamondoids are used to determine electronic, structural, and vibrational properties of GeSi superlattice nanocrystals and bulk as their building block limit. Density functional theory at the generalized gradient approximation level of Perdew, Burke, and Ernzerhof (PBE) with 6-31G(d) basis including polarization functions is used to investigate the electronic structure of these diamondoids. The investigated molecules and diamondoids range from GeSiH₆ to Ge₆₃Si₆₃H₉₂. The variation of the energy gap is shown from nearly 7 eV toward bulk value which is slightly higher than the average of Si and Ge energy gaps. Variations of bond lengths, tetrahedral, and dihedral angles as the number of atoms increases are shown taking into account the effect of shape fluctuations. Localized and delocalized electronic charge distribution and bonds for these molecules are discussed. Vibrational radial breathing mode (RBM) converges from its initial molecular value at 332 cm^?1 to its bulk limit at 0 cm^?1 (blue shift). Longitudinal optical-highest reduced mass mode (HRMM) converges from its initial molecular value 332 cm^?1 to experimental bulk limit at 420.7 cm^?1 (red shift). Hydrogen vibrational modes are nearly constant in their frequencies as the size of diamondoids increases in contrast with lower frequency Ge–Si vibrational modes. GeSi diamondoids can be identified from surface hydrogen vibrational modes fingerprint, while the size of these diamondoids can be identified from Ge–Si vibrational modes.     

Temperature dependence on phase evolution in the BaTiO3 polytypes studied using ab initio calculations

Identifying the forces that drive a temperature-induced phase transition is always challenging in the prospect of the first-principles methods. Herein, we perform a first-principles study of the temperature effects on structural, energetic, electronic, and vibrational properties of four BaTiO₃ polymorphs using quasi-harmonic approximations. Study of the stability between these four phases, which we break into contributions arising from the vibration of the lattice, electronic structure, and volume expansion/contraction, is helpful to confirm the sequence of phase transitions as cubic → tetragonal → orthorhombic → rhombohedral, as well as its transition temperatures. A general mechanism was proposed based on the combination between structural distortions at [TiO₆] clusters, vibrational characteristics, and electronic structure. These findings confirm the power of quasi-harmonic approximations to disclose the main fingerprints associated with both thermic and mechanical phase transitions, serving as a guide for further theoretical studies.     

Theoretical studies on the structural,spectroscopic, thermodynamic,and electronic properties of zoledronic acid

The structure, spectroscopic, thermodynamic, and electronic properties of zoledronic acid (ZL, 1-hydroxy- 2-(1H-imidazol-1-yl)ethane-1,1-diyldiphosphonic acid), typical third-generation nitrogen-containing bisphosphonates (N-BPs), have been investigated systematically. Six conformations are taken into account, including three unprotonated and three protonated structures. They are optimized by four different density functional theory (DFT) methods combined with four different basis sets to evaluate their performance in predicting the structural and spectral features of ZL. Thermodynamic properties are calculated based on the harmonic vibrational analysis, including the standard heat capacity (C _p,m ⁰ ), entropy (S _m ⁰ ), and enthalpy (S _m ⁰ ). The ¹H and ¹³C NMR chemical shifts are calculated using the GIAO method and compared with the experimental data. Molecular electrostatic potential (MEP) and frontier molecular orbital (FMO) analyses are also performed to study the electronic characteristics of the title compound.     

New Insights into the Band‐Gap Narrowing of (N,P)‐Codoped TiO2 from Hybrid Density Functional Theory Calculations

The electronic properties of anatase‐TiO₂ codoped by N and P at different concentrations have been investigated via generalized Kohn–Sham theory with the Heyd–Scuseria–Ernzerhof (HSE06) hybrid functional for exchange‐correlation in the context of density functional theory. At high doping concentrations, we find that the high photocatalytic activity of (N, P)‐codoped anatase TiO₂ vis‐à‐vis the N‐monodoped case can be rationalized by a double‐hole‐mediated coupling mechanism [Yin et al., Phys. Rev. Lett. 2011, 106, 066801] via the formation of an effective N? P bond. On the other hand, Ti³⁺ and Ti⁴⁺ ions’ spin double‐exchange results in more substantial gap narrowing for larger separations between N and P atoms. At low doping concentrations, double‐hole‐coupling is dominant, regardless of the N? P distance.     

Correlation functional in screened‐exchange density functional theory procedures

In the present study, we have explored several prospects for the further development of screened‐exchange density functional theory (SX‐DFT) procedures. Using the performance of HSE06 as our measure, we find that the use of alternative correlation functionals (as oppose to PBEc in HSE06) also yields adequate results for a diverse set of thermochemical properties. We have further examined the performance of new SX‐DFT procedures (termed HSEB‐type methods) that comprise the HSEx exchange and a (near‐optimal) reparametrized B97c (c _OS,0 = c _SS,0 = 1, c _OS,1 = −1.5, c _OS,2 = −0.644, c _SS,1 = −0.5, and c _SS,2 = 1.10) correlation functionals. The different variants of HSEB all perform comparably to or slightly better than the original HSE‐type procedures. These results, together with our fundamental analysis of correlation functionals, point toward various directions for advancing SX‐DFT methods. © 2017 Wiley Periodicals, Inc.     

Band gap anomalies of the ZnM2(III)O4 (M(III)=Co, Rh, Ir) spinels

Development of high figure-of-merit p-type transparent conducting oxides has become a global research goal. ZnM(2)(III)O(4) (M(III) = Co, Rh, Ir) spinels have been identified as potential p-type materials, with ZnIr(2)O(4) reported to be a transparent conducting oxide. In this article the geometry and electronic structure of ZnM(2)(III)O(4) are studied using the Perdew-Purke-Ernzerhof generalized gradient approximation (PBE-GGA) to density functional theory and a hybrid density functional, HSE06. The valence band features of all the spinels indicate that they are not conducive to high p-type ability, as there is insufficient dispersion at the valence band maxima. The trend of increasing band-gap as the atomic number of the M(III) cation increases, as postulated from ligand field theory, is not reproduced by either level of theory, and indeed is not seen experimentally in the literature. GGA underestimates the band-gaps of these materials, while HSE06 severely overestimates the band-gaps. The underestimation (overestimation) of the band-gaps by GGA (HSE06) and the reported transparency of ZnIr(2)O(4) is discussed.     

Efficient hybrid density functional calculations in solids: assessment of the Heyd-Scuseria-Ernzerhof screened Coulomb hybrid functional

The present work introduces an efficient screening technique to take advantage of the fast spatial decay of the short range Hartree-Fock (HF) exchange used in the Heyd-Scuseria-Ernzerhof (HSE) screened Coulomb hybrid density functional. The screened HF exchange decay properties and screening efficiency are compared with traditional hybrid functional calculations on solids. The HSE functional is then assessed using 21 metallic, semiconducting, and insulating solids. The examined properties include lattice constants, bulk moduli, and band gaps. The results obtained with HSE exhibit significantly smaller errors than pure density functional theory (DFT) calculations. For structural properties, the errors produced by HSE are up to 50% smaller than the errors of the local density approximation, PBE, and TPSS functionals used for comparison. When predicting band gaps of semiconductors, we found smaller errors with HSE, resulting in a mean absolute error of 0.2 eV (1.3 eV error for all pure DFT functionals). In addition, we present timing results which show the computational time requirements of HSE to be only a factor of 2-4 higher than pure DFT functionals. These results make HSE an attractive choice for calculations of all types of solids.     

An S‐Oxygenated [NiFe] Complex Modelling Sulfenate Intermediates of an O2‐Tolerant Hydrogenase

To understand the molecular details of O₂‐tolerant hydrogen cycling by a soluble NAD⁺‐reducing [NiFe] hydrogenase, we herein present the first bioinspired heterobimetallic S‐oxygenated [NiFe] complex as a structural and vibrational spectroscopic model for the oxygen‐inhibited [NiFe] active site. This compound and its non‐S‐oxygenated congener were fully characterized, and their electronic structures were elucidated in a combined experimental and theoretical study with emphasis on the bridging sulfenato moiety. Based on the vibrational spectroscopic properties of these complexes, we also propose novel strategies for exploring S‐oxygenated intermediates in hydrogenases and similar enzymes.     

Au(II)化合物[Au(CH2)2PH2]2X2 (X＝F, Cl, Br, I)的量子化学理论研究

采用从头计算HF, MP2方法和密度泛函理论, 对Au(II)系列化合物[Au(CH₂)₂PH₂]₂X₂ (X＝F, Cl, Br, I)的几何结构、电子结构和振动频率进行了研究. 研究表明Au的5d和6s电子参与Au—Au以及Au—X之间的成键. Au—Au, Au—X键强烈的电子相关作用使HF方法不适于该体系的研究, BP86和B3LYP两种泛函给出较大的Au—Au和Au—X键长, 而MP2方法和局域的密度泛函方法则给出了合理的结构参数. 局域密度泛函方法计算得到的Au—Au键和 Au—X键振动频率也与实验数据符合较好. 还运用含时密度泛函理论计算了[Au(CH₂)₂PH₂]₂X₂的电子激发能, 对分子在紫外-可见光谱范围内的电子跃迁进行了分析, 考察了卤素配体对激发能的影响, 并结合分子轨道能级的变化对此给予了解释.     

First-principles study of structural,elastic, electronic and optical properties of perovskites XCaH3 (X = Cs and Rb) under pressure

The stability, structural parameters, elastic constants, electronic and optical properties of perovskites CsCaH₃ and RbCaH₃ were investigated by the density functional theory. The calculated lattice parameters are in agreement with previous calculation and experimental data. The energy band structures, density of states, born-effective-charge and Mulliken charge population were obtained. The perovskites CsCaH₃ and RbCaH₃ present a direct band gap of 3.15 eV and 3.27 eV at equilibrium. The top of the valence bands reflects the s electronic character for both structures. Furthermore, the absorption spectrum, refractive index, extinction coefficient, reflectivity, energy-loss spectrum, and dielectric function were calculated. The origin of the spectral peaks was interpreted based on the electronic structures. The static dielectric constant and refractive index are indeed, inverse proportional to the direct band gap.     

On the issue of structure-chemical features of silicon oxides and thermal emission of singlet oxygen

Special features of the chemical structure of silicon oxides were analyzed by the methods of vibrational and electronic spectroscopy and it was found that cationic (Si-Si) and anionic (O₂ and O₃) sublattices contain isolated individual Si-CO oscillators, including electron-excited oscillators, and homoatomic bonds. By means of chromatographic and sorption-calorimetric analyses inverse dependence of singlet oxygen specific thermal emission on specific surface area of samples was found that corresponds to volume generation of excited states.     

Theoretical study of CeO2 and Ce2O3 using a screened hybrid density functional

The predicted structures and electronic properties of CeO(2) and Ce(2)O(3) have been studied using conventional and hybrid density functional theory. The lattice constant and bulk modulus for CeO(2) from local (LSDA) functionals are in good agreement with experiment, while the lattice parameter from a generalized gradient approximation (GGA) is too long. This situation is reversed for Ce(2)O(3), where the LSDA lattice constant is much too short, while the GGA result is in reasonable agreement with experiment. Significantly, the screened hybrid HSE functional gives excellent agreement with experimental lattice constants for both CeO(2) and Ce(2)O(3). All methods give insulating ground states for CeO(2) with gaps for the 4f band lying between 1.7 eV (LSDA) and 3.3 eV (HSE) and 6-8 eV for the conduction band. For Ce(2)O(3) the local and GGA functionals predict a semimetallic ground state with small (0-0.3 eV) band gap but weak ferromagnetic coupling between the Ce(+3) centers. By contrast, the HSE functional gives an insulating ground state with a band gap of 3.2 eV and antiferromagnetic coupling. Overall, the hybrid HSE functional gives a consistent picture of both the structural and electronic properties of CeO(2) and Ce(2)O(3) while treating the 4f band consistently in both oxides.