Assessment and validation of a screened Coulomb hybrid density functional

This paper presents a revised and improved version of the Heyd-Scuseria-Ernzerhof screened Coulomb hybrid functional. The performance of this functional is assessed on a variety of molecules for the prediction of enthalpies of formation, geometries, and vibrational frequencies, yielding results as good as or better than the successful PBE0 hybrid functional. Results for ionization potentials and electron affinities are of slightly lower quality but are still acceptable. The comprehensive test results presented here validate our assumption that the screened, short-range Hartree-Fock (HF) exchange exhibits all physically relevant properties of the full HF exchange. Thus, hybrids can be constructed which neglect the computationally demanding long-range part of HF exchange while still retaining the superior accuracy of hybrid functionals, compared to pure density functionals.     

A screened hybrid density functional study on energetic complexes: Metal carbohydrazide nitrates

The molecular geometry, electronic structure and thermochemistry of a series of metal carbohydrazide nitrates were investigated using the Heyd–Scuseria–Ernzerhof (HSE) screened hybrid density functional. The results show that Ca, Sr, and Ba complexes have additional coordinated oxygen atoms from the nitrate ion, which differ obviously from Cu, Ni, Co, and Mg complexes in terms of the geometric structure. Detailed NBO analyses clearly indicate that the metal–ligand interactions in Cu, Ni, and Co complexes are covalent, whereas those of Mg, Ca, Sr, and Ba complexes are ionic in nature. Furthermore, the donor–acceptor interactions result in a reduction of occupancies of σ_{C? O} and σ_{N? H} orbitals. Consequently, the bond lengths increase and the bond orders decrease. Finally, the calculated heats of formation predict that the ionic alkaline‐earth metal carbohydrazide nitrates are more stable than the covalent transition metal carbohydrazide nitrates. It agrees well with the available experimental thermal stabilities, indicating that the metal–ligand bonding character plays an important role in the stabilities of these energetic complexes. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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.     

A screened hybrid density functional study on energetic complexes: Alkaline-earth metal carbohydrazide perchlorates

The molecular geometries, electronic structures and stabilities of a series of alkaline-earth metal carbohydrazide perchlorates were investigated using the Heyd–Scuseria–Ernzerhof (HSE) screened hybrid density functional. The results show that Be and Mg complexes have six-coordinated octahedron features, as previously reported for the transition metal complexes. However, Ca, Sr and Ba complexes have additional coordinated oxygen atoms from the perchlorate ion. Detailed NBO analyses indicate that the metal–ligand interactions are essentially ionic and play an important role in the stabilities of these energetic complexes. The donor–acceptor interactions result in a reduction of occupancies of σ_C=O and σ_N–H bond orbitals, and also their subsequent impact on bond length and bond order.     

COS与O2反应的密度泛函理论研究

用量子化学密度泛函理论（DFT）方法，对COS与O2的反应进行了理论研究．在UB3LYP／6—31G^＊，UB3LYP／6—311＋＋G^＊＊水平上，优化了反应势能面上各驻点（反应物、产物、中间体和过渡态）的几何构型，在UB3LYP／6—31G^＊水平上通过内禀反应坐标（IRC）计算和振动分析，对过渡态进行了确认．在CCSD（T）／6—311＋＋G（2d，2p）水平上进行了单点能量计算，并确定了反应机理．研究结果表明，反应主要产物为CO2和SO.     

Adsorption and migration growth of the Ce,O, and F adatoms on the CeO2 (111) surface: A density functional theory study

In order to investigate the microscopic behavior of the crystal surface growth of the fluorinated cerium dioxide polishing powder, the adsorption and migration of the Ce, O, and F atoms on the CeO₂ (111) surface were studied by using density functional theory with Hubbard correction +U. The adsorption energies of three single atoms at five high-symmetry sites and the migration activation energies along the migration pathway on the CeO₂ (111) surface were calculated. Results show that the most stable adsorption sites of the Ce, O, and F atoms were the O_h, Ce_bri, and Ce_t sites, respectively. The Ce atom migrated from the O_h to the O_t site. The O atom migrated from the Ce_bri to the O_bri site. The F atom migrated from the Ce_t to the O_h site. The migration activation energies of the Ce, O, and F atoms along the migration pathways were 1.526, 0.597, and 0.263 eV, respectively. The F adatom does not change the spatial configuration of the Ce and the O atoms. When the O vacancy occurs on the CeO₂ (111) surface, the F adatom can make up for the O vacancy defect.     

Theoretical study of the benzene excimer using time-dependent density functional theory

A theoretical characterization of the potential energy surfaces of the singlet benzene excimer states derived from the B2u monomer excited state has been performed using time-dependent density functional theory. The excited-state potential energy surfaces were initially characterized by computations along the parallel and perpendicular intermolecular translational coordinates. These calculations predict that the lowest excited state for parallel translation is bound with a minimum at 3.15 angstroms and with a binding energy of 0.46 eV, while the perpendicular translational coordinate was essentially found to be a repulsive state. At the calculated minimum distance, the effect of in-plane rotation, out-of-plane rotation, and slipped-parallel translation were examined. The rotational calculations predict that deviations from the D6h geometry lead to a destabilization of the excimer state; however, small angular variations in the range of 0 degrees -10 degrees are predicted to be energetically feasible. The slipped-parallel translational calculations also predict a destabilizing effect on the excimer state and were found to possess barriers to this type of dissociation in the range of 0.50-0.61 eV. When compared to experimentally determined values for the benzene excimer energetics, the calculated values were found to be in semiquantitative agreement. Overall, this study suggests that the time-dependent density functional theory method can be used to characterize the potential energy surfaces and the energetics of aromatic excimers with reasonable accuracy.     

Dispersion in the Mott insulator UO2: A comparison of photoemission spectroscopy and screened hybrid density functional theory

We present a comparison between the screened hybrid density functional theory of Heyd, Scuseria, and Enzerhof (HSE06) and high-resolution photoemission (PES) measurement on a single crystal of UO(2). Angle-resolved photoemission data show a slight dispersion in the f-orbital derived bands in good agreement with the HSE band structure. The effect of spin-orbit coupling on the HSE band gap has also been calculated and found to be negligible.     

Structures and vibrations of Nb3O and Nb3O-: a density functional study

Density functional calculations of neutral and anionic niobium trimer monoxides are presented. The calculations were performed employing scalar quasirelativistic effective core potentials. In order to test the accuracy of the used effective core potentials in the framework of density functional theory the pulsed field ionization-zero electron kinetic energy photoelectron spectrum of Nb(3)O was simulated and compared to experiment. Different isomers of Nb(3)O and Nb(3)O(-) were studied in order to determine the ground state structures. For both neutral and anionic systems a planar C(2v) structure with an edge-bound oxygen atom was found as a ground state. Equilibrium structure parameters, harmonic frequencies, and adiabatic electron affinity are reported. The calculated electron affinity and frequencies are in good agreement with the available experimental data obtained recently from vibrationally resolved negative ion photoelectron spectroscopy.     

Adsorption and protonation of CO2 on partially hydroxylated gamma-Al2O3 surfaces: a density functional theory study

Adsorption and protonation of CO2 on the (110) and (100) surfaces of gamma-Al2O3 have been studied using density functional theory slab calculations. On the dry (110) and (100) surfaces, the O-Al bridge sites were found to be energetically favorable for CO2 adsorption. The adsorbed CO2 was bound in a bidentate configuration across the O-Al bridge sites, forming a carbonate species. The strongest binding with an adsorption energy of 0.80 eV occurs at the O3c-Al5c bridge site of the (100) surface. Dissociation of water across the O-Al bridge sites resulted in partially hydroxylated surfaces, and the dissociation is energetically favorable on both surfaces. Water dissociation on the (110) surface has a barrier of 0.42 eV, but the same process on the (100) surface has no barrier with respect to the isolated water molecule. On the partially hydroxylated gamma-Al2O3 surfaces, a bicarbonate species was formed by protonating the carbonate species with the protons from neighboring hydroxyl groups. The energy difference between the bicarbonate species and the coadsorbed bidentate carbonate species and hydroxyls is only 0.04 eV on the (110) surface, but the difference reaches 0.97 eV on the (100) surface. The activation barrier for forming the bicarbonate species on the (100) surface, 0.42 eV, is also lower than that on the (110) surface (0.53 eV).     

Microstructural evolution in a CeO2-Gd2O3 system

Microstructural evolution in a CeO2-Gd2O3 system at atomic and nanoscale levels with increasing Gd concentration has been comprehensively investigated by transmission electron microscopy. When the Gd concentration was increased from 10 to 80 at.%, the phase transformation from ceria with fluorite structure to solid solution with C-type structure was not a sudden change but an evolution in the sequence of clusters, domains, and precipitates with C-type structure in the fluorite-structured matrix. Moreover, the ordering of aggregated Gd cations and oxygen vacancies in these microstructural inhomogeneities developed continuously with increasing Gd concentration. This microstructural evolution can be further described based on the development of defect clusters containing Gd cations and oxygen vacancies.     

Theoretical study of the syn and anti thiophene-2-aldehyde conformers using density functional theory and normal coordinate analysis

An extensive computational study of thiophene-2-aldehyde conformers syn and anti has been carried out using density functional (DFT). From these calculations, B3LYP/6-31G(d) has been chosen as it produces results remarkably close in comparison with experimental ones, with less demanding computational time. Data obtained from DFT computation were used to perform a normal coordinate analysis to complement and give insight in the experimental vibrational assignment. Calculated dipole moments and relative stabilities of isomers coherently support experimental statements given in the literature.     

Syngas production from methane over CeO2-Fe2O3 mixed oxides using a chemical-looping method

A series of mixed oxides Ce_{1 ? x} Fe _x O₂ was prepared by a hydrothermal method. XRD and Raman spectra were measured to study the structure of the prepared materials. The temperature-programmed reduction was undertaken to estimate reducibility of the oxides. Syngas generation from methane using these materials as oxygen carriers/catalysts via a chemical-looping procedure was investigated in detail. This procedure includes catalytic oxidation and decomposition of methane to produce H₂-rich gas at the first step followed by the production of the CO-rich gas by oxidizing the carbon deposited on deactivated catalysts. The results showed that all iron ions were incorporated into the ceria lattice with the formation of oxygen vacancies in the Ce_0.9Fe_0.1O₂ sample, while isolated Fe₂O₃ particles were distributed on the surface of the Ce_0.8Fe_0.2O₂ sample. TPR measurements and the analysis of the two-step chemical-looping reactions indicated a strong interaction between the Ce and Fe species which accounts for an increased activity of the mixed oxides in the syngas generation compared to that of individual oxides. Among the several samples, the Ce_0.8Fe_0.2O₂ catalyst showed the highest activity for methane partial oxidation due to the synergetic effects caused by the interaction of surface iron entities and Ce-Fe solid solution. In addition, selective oxidation of carbon by oxygen to CO can also be found over this material since gaseous products are formed at the carbon oxidation step with the selectivity to CO reaching 91.2%. Evidence is presented that syngas can be feasibly produced from methane with high selectivity via the chemical-looping procedure over the CeO₂-Fe₂O₃ mixed oxides.     

Copper impurities in bulk ZnO: a hybrid density functional study

Transition metal doping of ZnO is considered as a promising way to obtain a diluted magnetic semiconducting oxide. In this work we investigate copper doping of ZnO by means of density functional theory, using a hybrid exchange-correlation functional and a periodic approach with localized atomic basis functions. Isolated copper species, such as copper substitutional to zinc, Cu(s), and Cu interstitial, Cu(i), are analyzed in terms of transition energy levels and hyperfine coupling constants with reference to available spectroscopic data. We also examine the potential magnetic interaction between copper species, their interaction with oxygen vacancies, and the possibility of copper clustering. The relative stability of the various copper impurities considered in this study is finally compared on the basis of their formation energy at different oxygen chemical potentials and Fermi level values.     

水在Au1/CeO2单原子催化剂上吸附与解离的理论研究

自Haruta和Hutchings发现负载的纳米金催化剂的催化活性后,负载型金催化剂一直是非均相催化的研究重点之一.近年来,单原子催化剂因其优异的活性、选择性,超高的原子利用效率,引起了科学家们的广泛关注.越来越多的单原子金催化剂被成功制备,并被证实具有很好的催化活性.水,作为环境中最常见的物质,在实际的催化体系中往往难以避免,即使在超高真空环境中也会有痕量的水气存在.水的解离不仅是水煤气反应的重要步骤之一,而且对别的反应也有一定的促进作用.尽管水和纳米团簇催化剂之间的研究已经颇有成效,但水和单原子金催化剂之间的作用还不是非常清晰.因此,我们采用密度泛函理论从原子尺度研究了水和Au1/CeO2单原子催化剂的相互作用.我们首先研究了水在完美CeO2表面和含有一个氧空位的CeO2–x表面上的解离过程,研究发现分子态的水和解离态的水在完美CeO2表面可以共存,而一旦在表面形成氧空位后,由于较低的能垒和极大的放热,解离态的水将占据绝对优势.接下来探索了水在完美Au1/CeO2表面和含有一个氧空位的Au1/CeO2–x表面上的解离过程,发现结论恰好和CeO2表面相反.水的解离过程在完美的Au1/CeO2表面几乎是一个无能垒的过程,并且解离会放出大量的热量.而一旦在表面形成氧空位后,单原子Au的轨道处于满占状态,无法提供水的吸附位点.水的解离过程在Ce位点进行,分子吸附能与解离吸附能相当,分子态与解离态共存.为了进一步理解单原子金在水的解离过程中起到的作用,我们分析了水和Au1/CeO2之间的电子相互作用.研究结果表明,单原子金不仅为水的吸附提供了位点,金的5d轨道和水的2p轨道之间的相互作用还有效减弱了水中氧氢键的强度,使水的解离更容易进行.由此可见,在涉及到水解离的反应中,以Au1/CeO2为代表的单原子催化剂有望带来新的突破.最后,我们还测试了范德华力对研究体系的影响.研究发现尽管范德华力会使吸附能的绝对值增加,但是并不影响我们得到的结论.     

Theoretical study on infrared vibrational spectra of boric-acid in gas-phase using density functional methods

Density functional theory (DFT) methods with various exchange-correlation functionals such as SVWN, BVWN, BVWN5, BLYP, B1LYP, B3LYP, B3PW91, and BH and H are employed in a theoretical study of molecular boric-acid in gas-phase. In the calculations, the split valence 6-311++G** and 6-31G* basis sets were used. The geometry, zero-point vibrational energies (ZPVEs), and harmonic infrared vibrational (IR) frequencies are predicted. The calculated C_3h-symmetry geometrical parameters are compared with Hartree–Fock (HF) calculation results and experimental data. IR frequencies predicted by the BLYP, B3LYP, and B3PW91 calculations are in good agreement with experimental data. The frequency calculations presented here also suggest that the C_3h-symmetrical structure corresponds to a minimum in the potential energy surface, but neither is D_3h- or C₃-symmetrical structure.     

Reactivity of Al3O3- cluster toward H2O studied by density functional theory

Density functional theory calculations (Becke's three parameter hybrid functional) have been done on a wide range of possible structures for the complexes formed in the reaction between Al(3)O(3) (-) and one or two water molecules. Both energetically competitive structural isomers of Al(3)O(3) (-) (kitelike and distorted rectangle) were considered. The structures of neutral complexes accessed from detachment of the stable anion structures were also optimized. The calculations predict that hydroxide complexes are energetically favored over Lewis acid-base and charge-dipole complexes. For Al(3)O(3) (-)/H(2)O complexes, the kite-based hydroxide and rectangle-based hydroxide are predicted to be nearly isoenergetic, while for Al(3)O(3) (-)/(H(2)O)(2), the rectangle-based dihydroxide emerges as being 0.5 eV more stable than the lowest energy kite-based dihydroxide. The structures of these and their neutrals are used to analyze anion PE spectra of Al(3)O(4)H(2) (-) and Al(3)O(5)H(4) (-) obtained previously [F. A. Akin and C. C. Jarrold, J. Chem. Phys. 118, 5841 (2003)].     

Y2O3,CeO2掺杂的Mg-PSZ陶瓷材料研究

采用传统的陶瓷工艺制备少量掺杂Ｙ２Ｏ３,ＣｅＯ２的Ｍｇ－ＰＳＺ陶瓷,材料在较低的温度下烧结致密并实现了微晶化,探讨了Ｙ２Ｏ３,ＣｅＯ２的复合稳定作用和１１００℃热处理过程对材料相组成、显微结构和力学性能的影响。经１１００℃适当时间热处理,Ｙ２Ｏ３,ＣｅＯ２的复合稳定作用有效报制亚共析分解反应发生,优化调整ｃ－ＺｒＯ２晶粒中的ｔ－ＺｒＯ２析出体成核长大过程。断裂特征为穿晶、沿晶兼有,相变增韧和微裂纹     

A dual-level ab initio and hybrid density functional theory dynamics study on the unimolecular decomposition reaction C2H5O-->CH2O + CH3

We present a direct ab initio and hybrid density functional theory dynamics study of the thermal rate constants of the unimolecular decomposition reaction of C2H5O-->CH2O + CH3 at a high-pressure limit. MPW1K/6-31+G(d,p), MP2/6-31+G(d,p), and MP2(full)/6-31G(d) methods were employed to optimize the geometries of all stationary points and to calculate the minimum energy path (MEP). The energies of all the stationary points were refined at a series of multicoefficient and multilevel methods. Among all methods, the QCISD(T)/aug-cc-pVTZ energies are in good agreement with the available experimental data. The rate constants were evaluated based on the energetics from the QCISD(T)/aug-cc-pVTZ//MPW1K/6-31+G(d,p) level of theory using both microcanonical variational transition state theory (microVT) and RRKM theory with the Eckart tunneling correction in the temperature range of 300-2500 K. The calculated rate constants at the QCISD(T)/aug-cc-pVTZ/MPW1K/6-31+G(d,p) level of theory are in good consistent with experimental data. The fitted three-parameter Arrhenius expression from the microVT/Eckart rate constants in the temperature range 200-2500 K is k = 2.52 x 10(12)T(0.41)e(-8894.0/T) s(-1). The falloff curves of pressure-dependent rate constants are performed using master-equation method within the temperature range of 391-471 K. The calculated results are in good agreement with the available experimental data.