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.

A density functional theory study of CO oxidation on Ru(0001) at low coverage

We have performed ab initio density functional theory calculations with the generalized gradient approximation to investigate CO oxidation on Ru(0001). Several reaction pathways and transition states are identified. A much higher reaction barrier compared to that on Pt(111) is determined, confirming that the Ru is very inactive for CO oxidation under UHV conditions. The origin of the reaction barrier was analyzed. It is found that in the transition state the chemisorbed O atom sits in an unfavorable bonding site and a significant competition for bonding with the same substrate atoms occurs between the CO and the chemisorbed O, resulting in the high barrier. Ab initio molecular dynamics calculations show that the activation of the chemisorbed O atom from the initial hcp hollow site (the most stable site) to the bridge site is the crucial step for the reaction. The CO oxidation on Ru(0001) via the Eley-Rideal mechanism has also been investigated. A comparison with previous theoretical work has been made. (C) 2000 American Institute of Physics. [S0021-9606(00)31223-5].

Site symmetry dependence of repulsive interactions between chemisorbed oxygen atoms on Pt{100}-(1x1)

Ab initio total energy calculations using density functional theory with the generalized gradient approximation have been performed for the chemisorption of oxygen atoms on a Pt{100}-(1 x 1) slab. Binding energies for the adsorption of oxygen on different high-symmetry sites are presented. The bridge site is the most stable at a coverage of 0.5 ML, followed by the fourfold hollow site. The atop site is the least stable. This finding is rationalized by analyzing the ''local structures'' formed upon oxygen chemisorption. The binding energies and heats of adsorption at different oxygen coverages show that pairwise repulsive interactions are considerably stronger between oxygen atoms occupying fourfold sites than those occupying bridge sites. Analysis of the partial charge densities associated with Bloch states demonstrates that the O-Pt bond is considerably more localized at the bridge site. These effects cause a sharp drop in the heats of adsorption for oxygen on hollow sites when the coverage is increased from 0.25 to 0.5 ML. Mixing between oxygen p orbitals and Pt d orbitals can be observed over the whole metal d-band energy range.

Fluxes of organic contaminants from the river catchment into, through and out of the Humber Estuary, UK

This paper summarises the work done on the distribution and reactivity of organic contaminants (simazine, atrazine, lindane, fluoranthene, pyrene, PCB 77, PCB 118) in the Humber Estuary and associated major rivers, as part of the LOIS programme. The preliminary flux calculations show that the most important contributors of selected organic contaminants were the rivers Trent (45% of simazine, 20% of atrazine), Aire (30% of simazine and 33% of atrazine), Don (36 and 37% of fluoranthene and pyrene) and Ouse (18% of fluoranthene and pyrene). For lindane and PCBs, the Aire and Ouse were the key sources. The water flow in all the rivers shows strong seasonal variations, as do the contaminant concentrations. As a result, the mean daily fluxes of these contaminants displayed a strong seasonality. Annual mean concentrations of simazine and atrazine decreased by more than 50% over the period 1994-1995 in most of the rivers, probably as a result of their restricted use in the UK. Mass balance calculations show that the Humber is a sink for atrazine, lindane, PCB 77 and PCB 118, although the degree of removal is generally much lower for atrazine and lindane than for PCB 77 and PCB 118. Mass balance results also show that the Humber can either be a source of fluoranthene and pyrene (in the suspended particulate phase), or a sink (in the dissolved phase), although overall the Humber acts as sink. The budget exercise represents an attempt to quantify the input and output of selected organic contaminants from catchment to ocean. However, due to limited data and assumptions involved in calculations, the estimates should be considered as an order of magnitude approximation. Further improvement both in resolution and accuracy is required.

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.

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.

Intermediate-coupling R -matrix calculations of electron-impact excitation of Fe 5+

For applications to laboratory and astrophysical plasmas, there is a great need for accurate electron-impact excitation data between individual levels in the lower charge-state ions of iron. Recently, we have reported on the first intermediate-coupling R -matrix calculation of electron-impact excitation in Fe 4+ , in which the close-coupling expansion of the target included levels from both ground and excited configurations (Ballance et al 2007 J. Phys. B: At. Mol. Opt. Phys. [/0953-4075/40/23/f01] 40 F327 , 2008 Europhys. News 39 14). In this paper, we present the results of two large intermediate-coupling Dirac R -matrix calculations of electron-impact excitation of Fe 5+ . The results from the two calculations, which differ only in the configuration–interaction expansions of the target, are compared. These comparisons provide some indication of the accuracy of the calculations and the resulting data should be useful in modelling plasmas containing iron.

Electron-impact excitation of Fe2: a comparison of intermediate coupling frame transformation, Breit-Pauli and Dirac R-matrix calculations

Modeling the spectral emission of low-charge iron group ions enables the diagnostic determination of the local physical conditions of many cool plasma environments such as those found in H II regions, planetary nebulae, active galactic nuclei etc. Electron-impact excitation drives the population of the emitting levels and, hence, their emissivities. By carrying-out Breit-Pauli and intermediate coupling frame transformation (ICFT) R-matrix calculations for the electron-impact excitation of Fe$^{2+}$ which both use the exact same atomic structure and the same close-coupling expansion, we demonstrate the validity of the application of the powerful ICFT method to low-charge iron group ions. This is in contradiction to the finding of Bautista et al. [Ap.J.Lett, 718, L189, (2010)] who carried-out ICFT and Dirac R-matrix calculations for the same ion. We discuss possible reasons.

Fine-structure photoionization cross sections of Fe II

We present the first calculation of fine-structure photoionization cross sections for the ground state of singly ionized Fe. These large-scale ab initio calculations, limited to the near-threshold region, were performed in the close-coupling approximation using a Dirac–Coulomb R -matrix method implemented within a modified version of the DARC package. Our calculated cross sections reproduce in detail the resonance structures observed in previous experimental determinations.

PAMOP: Petascale atomic, molecular and optical collision calculations

Petaflop architectures are currently being utilized efficiently to perform large scale computations in Atomic, Molecular and Optical Collisions. We solve the Schr\"odinger or Dirac equation for the appropriate collision problem using the R-matrix or R-matrix with pseudo-states approach. We briefly outline the parallel methodology used and implemented for the current suite of Breit-Pauli and DARC codes. In this report, various examples are shown from our theoretical results compared with experimental results obtained from Synchrotron Radiation facilities where the Cray architecture at HLRS is playing an integral part in our computational projects.

Dirac R -matrix calculations of electron-impact excitation of neon-like krypton

We have employed the Dirac R -matrix method to determine electron-impact excitation cross sections and effective collision strengths in Ne-like Kr 26+ . Both the configuration-interaction expansion of the target and the close-coupling expansion employed in the scattering calculation included 139 levels up through n = 5. Many of the cross sections are found to exhibit very strong resonances, yet the effects of radiation damping on the resonance contributions are relatively small. Using these collisional data along with multi-configuration Dirac–Fock radiative rates, we have performed collisional-radiative modeling calculations to determine line-intensity ratios for various radiative transitions that have been employed for diagnostics of other Ne-like ions.