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.

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.

Photoionization cross section calculations for the halogen-like ions Kr + and Xe +

Photoionization cross section calculations on the halogen-like ions; Kr + and Xe + have been performed for a photon energy range from each ion threshold to 15 eV, using large-scale close-coupling calculations within the Dirac--Coulomb R -matrix approximation. The results from our theoretical work are compared with recent measurements made at the ASTRID merged-beam set-up at the University of Aarhus in Denmark and from the Fourier transform ion cyclotron resonance trap method at the SOLEIL synchrotron radiation facility in Saint-Aubin, France Bizau et al (2011 J. Phys. B: At. Mol. Opt. Phys. 44 055205) and the advanced light source M{ü}ller (2012 private communication), Aguliar et al (2012 J. Phys.: Conf. Ser . at press). For each of these complex ions our theoretical cross section results over the photon energy range investigated are seen to be in excellent agreement with experiment. Resonance energy positions and quantum defects of the prominent Rydberg resonances series identified in the spectra are compared with experiment for these complex halogen-like ions.

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.

The 2p interloper resonances in collisions

Ab initio electron scattering calculations using the R -matrix approach have been performed for within a three-state valence configuration-interaction model (VCI). The lowest three electronic target states ( , and the ) of this molecular nitrogen cation are included in the close-coupling method, with each state being represented by a valence CI approximation. From a detailed analysis of the resonance structure found in our work for the symmetries we find four prominent Rydberg series of the type , , , and a interloper resonance. This interloper molecular resonance associated with the B state of is seen to cause distortions of the resulting resonance spectra. A comparison of our total cross sections for the X - B transition shows excellent agreement with the available experimental data.

Collisional-radiative calculations of He line emission in low-temperature plasmas

We present spectral modeling results for neutral helium. Our underlying atomic data contains radiative transition rates that are generated from atomic structure calculations and electron-impact excitation rates, that are determined from both the standard R-matrix method and the R-matrix with pseudostates RMPS method. In this paper, we focus on transitions of particular importance to diagnostic line ratios. For example, our calculated rate coefficient for the electron-impact transition 1s3s 1S→1s3p 1P, which has a pronounced effect on the 728.1 nm diagnostic spectral line, is found to be in good agreement with previous experimental mea- surements. We also consider transitions from the 1s2 1S ground and 1s2s 3S terms to terms of the n=4 shell. They are found to be affected significantly by coupling of the bound states to the target continuum continuum coupling, which is included in our RMPS calculation, but not in our standard R-matrix calculation. We perform collisional-radiative calculations to determine spectral line intensity ratios for three ratios of particular interest, namely the 504.8 nm/471.3 nm, 492.2 nm/471.3 nm, and 492.2 nm/504.8 nm line ratios. Comparing our results determined from the RMPS excitation rates with those from the standard R-matrix excitation rates, we find that continuum coupling affects the rate coefficients significantly, leading to different values for all three line ratios. We also compare our modeling results with spectral measurements taken recently on the Auburn Helicon plasma device, finding that the ground and metastable populations are not in equilibrium, and that the experimental measurements are more consistent with the 1s2s 3S metastable term populations being short lived in the plasma.