Quasiparticle calculations of the electronic properties of ZrO2 and HfO2 polymorphs and their interface with Si

Quasiparticle calculations are performed to investigate the electronic band structures of various polymorphs of Hf and Zr oxides. The corrections with respect to density-functional-theory results are found to depend only weakly on the crystal structure. Based on these bulk calculations as well as those for bulk Si, the effect of quasiparticle corrections is also investigated for the band offsets at the interface between these oxides and Si assuming that the lineup of the potential at the interface is reproduced correctly within density-functional theory. On the one hand, the valence-band offsets are practically unchanged with a correction of a few tenths of electron volts. On the other hand, conduction-band offsets are raised by 1.3-1.5 eV. When applied to existing calculations for the offsets at the density-functional-theory level, our quasiparticle corrections provide results in good agreement with the experiment.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

Mechanistic Study of 1,3-Butadiene Formation in Acetylene Hydrogenation over the Pd-Based Catalysts Using Density Functional Calculations

Green oil, which leads to the deactivation of the catalysts used for the selective hydrogenation of acetylene, has long been observed but its formation mechanism is not fully understood. In this work, the formation of 1,3-butadiene, known to be the precursor of green oil, on both Pd(111) and Pd(211) surfaces is examined using density functional theory calculations. The pathways containing C-2 + C-2 coupling reactions as well as the corresponding hydrogenation reactions are studied in detail. Three pathways for 1,3-butadiene production, namely coupling plus hydrogenation and further hydrogenation, hydrogenation plus coupling plus hydrogenation, and a two step hydrogenation followed by coupling, are determined. By comparing the effective barriers, we identify the favored pathway on both surfaces. A general understanding toward the deactivation process of the industrial catalysts is also provided. In addition, the effects of the formation of subsurface carbon atoms as well as the Ag alloying on the 1,3-butadiene formation on Pd-based catalysts are also investigated and compared with experimental results.

Monte Carlo calculations of quantum yield in inhomogeneous PtSi/p-Si Schottky barriers

Monte Carlo calculations of quantum yield in PtSi/p-Si infrared detectors are carried out taking into account the presence of a spatially distributed barrier potential. In the 1-4 mu m wavelength range it is found that the spatial inhomogeneity of the barrier has no significant effect on the overall device photoresponse. However, above lambda = 4.0 mu m and particularly as the cut-off wavelength (lambda approximate to 5.5 mu m) is approached, these calculations reveal a difference between the homogeneous and inhomogeneous barrier photoresponse which becomes increasingly significant and exceeds 50% at lambda = 5.3 mu m. It is, in fact, the inhomogeneous barrier which displays an increased photoyield, a feature that is confirmed by approximate analytical calculations assuming a symmetric Gaussian spatial distribution of the barrier. Furthermore, the importance of the silicide layer thickness in optimizing device efficiency is underlined as a trade-off between maximizing light absorption in the silicide layer and optimizing the internal yield. The results presented here address important features which determine the photoyield of PtSi/Si Schottky diodes at energies below the Si absorption edge and just above the Schottky barrier height in particular.

Identifying the distinct features of geometric structures for hole trapping to generate radicals on rutile TiO2(110) in photooxidation using density functional theory calculations with hybrid functional

Using density functional theory calculations with HSE 06 functional, we obtained the structures of spin-polarized radicals on rutile TiO2(110), which is crucial to understand the photooxidation at the atomic level, and further calculate the thermodynamic stabilities of these radicals. By analyzing the results, we identify the structural features for hole trapping in the system, and reveal the mutual effects among the geometric structures, the energy levels of trapped hole states and their hole trapping capacities. Furthermore, the results from HSE 06 functional are compared to those from DFT + U and the stability trend of radicals against the number of slabs is tested. The effect of trapped holes on two important steps of the oxygen evolution reaction, i.e. water dissociation and the oxygen removal, is investigated and discussed.

Identifying the trend of reactivity for sp2 materials: an electron delocalization model from first principles calculations

The reactivity of sp² carbon materials is studied using the adsorption and dissociation of O₂ on graphene and graphene oxide as model systems. The reactions on the basal plane, zigzag and armchair edges of graphene and graphene oxide with different oxygen-containing groups are calculated using first principles calculations. Two Brønsted-Evans- Polanyi relationships are identified and an electron delocalization model is suggested to understand the general trend of reactivity for sp² carbon materials.

Some Understanding of Fischer-Tropsch Synthesis from Density Functional Theory Calculations

The combination of density functional theory (DFT) calculations and kinetic analyses is a very useful approach to study surface reactions in heterogeneous catalysis. The present paper reviews some recent work applying this approach to Fischer-Tropsch (FT) synthesis. Emphasis is placed on the following fundamental issues in FT synthesis: (i) reactive sites for both hydrogenation and C-C coupling reactions; (ii) reaction mechanisms including carbene mechanism, CO-insertion mechanism and hydroxyl-carbene mechanism; (iii) selectivity with a focus on CH(4) selectivity, alpha-olefin selectivity and chain growth probability; and (iv) activity.

NO Oxidation on Platinum Group Metals Oxides:First Principles Calculations Combined with Microkinetic Analysis

By combining density functional theory calculation and microkinetic analysis, NO oxidation on the platinum group metal oxides (PtO(2), IrO(2), OsO(2)) is investigated, aiming at shedding light on the activities of metal oxides and exploring the activity variations of metal oxides compared to their corresponding metals. A microkinetic model, taking into account the possible low diffusion of surface species on metal oxide surfaces, is proposed for NO oxidation. The resultant turnover frequencies of NO oxidation show that under the typical experimental condition, T = 600 K, p(O2) = 0.1 atm, p(NO) = 3 x 10(-4) atm, p(NO2) = 1.7 x 10(-4) atm; (i) IrO(2)(110) exhibits higher activity than PtO(2)(110) and OsO(2)(110), and (ii) compared to the corresponding metallic Pt, Ir, and Os, the activity of PtO(2) to catalyze NO oxidation is lower, but interestingly IrO(2) and OsO(2) exhibit higher activities. The reasons for the activity differences between the metals and oxides are addressed. Moreover, other possible reaction pathways of NO oxidation on PtO(2)(110), involving O(2) molecule (NO + O(2) -> OONO) and lattice bridge-O(2c), are also found to give low activities. The origin of the Pt catalyst deactivation is also discussed.

GRADIENT CORRECTIONS IN DENSITY-FUNCTIONAL THEORY CALCULATIONS FOR SURFACES - CO ON PD(110)

Ab initio total energy calculations have been performed for CO chemisorption on Pd(110). Local density approximation (LDA) calculations yield chemisorption energies which are significantly higher than experimental values but inclusion of the generalised gradient approximation (GGA) gives better agreement. In general, sites with higher coordination of the adsorbate to surface atoms lead to a larger degree of overbinding with LDA, and give larger corrections with GGA. The reason is discussed using a first-order perturbation approximation. It is concluded that this may be a general failure of LDA for chemisorption energy calculations. This conclusion may be extended to many surface calculations, such as potential energy surfaces for diffusion.

Strong second harmonic generation in SiC, ZnO, GaN two-dimensional hexagonal crystals from first-principles many-body calculations

The second harmonic generation (SHG) intensity spectrum of SiC, ZnO, GaN two-dimensional hexagonal crystals is calculated by using a real-time first-principles approach based on Green's function theory [Attaccalite et al., Phys. Rev. B: Condens. Matter Mater. Phys. 2013 88, 235113]. This approach allows one to go beyond the independent particle description used in standard first-principles nonlinear optics calculations by including quasiparticle corrections (by means of the GW approximation), crystal local field effects and excitonic effects. Our results show that the SHG spectra obtained using the latter approach differ significantly from their independent particle counterparts. In particular they show strong excitonic resonances at which the SHG intensity is about two times stronger than within the independent particle approximation. All the systems studied (whose stabilities have been predicted theoretically) are transparent and at the same time exhibit a remarkable SHG intensity in the range of frequencies at which Ti:sapphire and Nd:YAG lasers operate; thus they can be of interest for nanoscale nonlinear frequency conversion devices. Specifically the SHG intensity at 800 nm (1.55 eV) ranges from about 40-80 pm V(-1) in ZnO and GaN to 0.6 nm V(-1) in SiC. The latter value in particular is 1 order of magnitude larger than values in standard nonlinear crystals.

Some Attempts in the Rational Design of Heterogeneous Catalysts Using Density Functional Theory Calculations

Heterogeneous catalysis is of great importance both industrially and academically. Rational design of heterogeneous catalysts is highly desirable, and the computational screening and design method is one of the most promising approaches for rational design of heterogeneous catalysts. Herein, we review some attempts towards the rational catalyst design using density functional theory from our group. Some general relationships and theories on the activity and selectivity are covered, such as the Brønsted–Evans–Polanyi relation, volcano curves/surfaces, chemical potentials, optimal adsorption energy window and energy descriptor of selectivity. Furthermore, the relations of these relationships and theories to the rational design are discussed, and some examples of computational screening and design method are given.

Volumetric Carbon Content Calculations: A Geostatistical and Compositional Approach

This study applies spatial statistical techniques including cokriging to integrate airborne geophysical (radiometric) data with ground-based measurements of peat depth and soil organic carbon (SOC) to monitor change in peat cover for carbon stock calculations. The research is part of the EU funded Tellus Border project and is supported by the INTERREG IVA development programme of the European Regional Development Fund, which is managed by the Special EU Programmes Body (SEUPB). The premise is that saturated peat attenuates the radiometric signal from underlying soils and rocks. Contemporaneous ground-based measurements were collected to corroborate mapped estimates and develop a statistical model for volumetric carbon content (VCC) to 0.5 metres. Field measurements included ground penetrating radar, gamma ray spectrometry and a soil sampling methodology which measured bulk density and soil moisture to determine VCC. One aim of the study was to explore whether airborne radiometric survey data can be used to establish VCC across a region. To account for the footprint of airborne radiometric data, five cores were obtained at each soil sampling location: one at the centre of the ground radiometric equivalent sample location and one at each of the four corners 20 metres apart. This soil sampling strategy replicated the methodology deployed for the Tellus Border geochemistry survey. Two key issues will be discussed from this work. The first addresses the integration of different sampling supports for airborne and ground measured data and the second discusses the compositional nature of the VOC data.

Towards a spectroscopically accurate set of potentials for heavy hydride laser cooling candidates: effective core potential calculations of BaH

BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio $^{2}\Sigma^+$ potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction (MRCI+Q)is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron Effective Core-Potential (ECP) and even-tempered augmented polarized core-valence basis sets (aug-pCV$n$Z-PP, n= 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy De for the X$^2\Sigma^+$ state (extrapolated to the CBS limit) is 16895.12 cm$^{-1}$ (2.094 eV), which agrees within 0.1$\%$ of a revised experimental value of <16910.6 cm$^{-1}$, while the calculated re is within 0.03 pm of the experimental result.

Relativistic radiatively damped R -matrix calculations of the electron-impact excitation of helium-like iron and krypton

Electron-impact excitation data for He-like ions are of significant importance for diagnostic applications to both laboratory and astrophysical plasmas. Here we report on the first fully relativistic R -matrix calculations with radiation damping for the He-like ions Fe 24+ and Kr 34+ . Effective collision strengths for these two ions have been determined with and without damping over a wide temperature range for all transitions between the 49 levels through n = 5. We find that damping has a pronounced effect on the effective collision strengths for excitation to some of the low-lying levels, but its effect on excitation to the vast majority of levels is small. At the energy of a resonance peak, we also investigate the effect of radiation damping on the angular distribution of scattered electrons. Finally, we compare our results for Fe 24+ with an earlier intermediate coupling frame transformation R -matrix calculation with radiation damping by Whiteford et al ( J. Phys. B: At. Mol. Opt. Phys. 34 3179) and find good agreement, especially for excitation to the lower levels.