A DFT study of the chain growth probability in Fischer-Tropsch synthesis

The chain growth probability (alpha value) is one of the most significant parameters in Fischer-Tropsch (FT) synthesis. To gain insight into the chain growth probability, we systematically studied the hydrogenation and C-C coupling reactions with different chain lengths on the stepped Co(0001) surface using density functional theory calculations. Our findings elucidate the relationship between the barriers of these elementary reactions and the chain length. Moreover, we derived a general expression of the chain growth probability and investigated the behavior of the alpha value observed experimentally. The high methane yield results from the lower chain growth rate for C-1 + C-1 coupling compared with the other coupling reactions. After C-1, the deviation of product distribution in FT synthesis from the Anderson-Schulz-Flory distribution is due to the chain length-dependent paraffin/olefin ratio. (C) 2008 Elsevier Inc. All rights reserved.

Chain growth mechanism in Fischer-Tropsch synthesis: A DFT study of C-C coupling over Ru, Fe, Rh, and Re surfaces

A quantitative approach is used to understand the chain growth mechanism in FT synthesis on the Ru, Fe, Rh, and Re surfaces. The C-C coupling reactions are extensively calculated on the stepped metal surfaces. Combining the coupling barriers and reactant stabilities, we investigate the reaction rates of all possible C, + C-1 coupling pathways on the metal surfaces. It is found that (i) all the transition-state structures are similar on these surfaces, while some coupling barriers are very different; (ii) the dominant chain growth pathways on these surfaces are different: C + CH and CH + CH on Rh and Ru surfaces, C + CH3 on Fe surface, and C + CH on Re surface. The common features of the major coupling reactions together with those on the Co surface are also discussed.

A density functional theory study of the a-olefin selectivity in Fischer-Tropsch synthesis

We studied the alpha-olefin selectivity in Fischer-Tropsch (FT) synthesis using density functional theory (131717) calculations. We calculated the relevant elementary steps from C-2 to C-6 species. Our results showed that the barriers of hydrogenation and dehydrogenation reactions were constant with different chain lengths, and the chemisorption energies of alpha-olefins from DFT calculations also were very similar, except for C-2 species. A simple expression of the paraffin/olefin ratio was obtained based on a kinetic model. Combining the expression of the paraffin/olefin ratio and our calculation results, experimental findings are satisfactorily explained. We found that the physical origin of the chain length dependence of paraffin/olefin ratio is the chain length dependence of both the van der Waals interaction between adsorbed alpha-olefins and metal surfaces and the entropy difference between adsorbed and gaseous alpha-olefins, and that the greater chemisorption energy of ethylene is the main reason for the abnormal ethane/ethylene ratio. (c) 2008 Elsevier Inc. All rights reserved.

Reactivity of the 4d transition metals toward N hydrogenation and NH dissociation: A DFT-based HSAB analysis

The density functional theory (DFT) based hard-soft acid-base (HSAB) reactivity indices, including the electrophilicity index, have been successfully applied to many areas of molecular chemistry. In this work we test the applicability of such an approach to fundamental surface chemistry. We have considered, as prototypical surface reactions, both the hydrogenation of atomic nitrogen and the dissociative adsorption of the NH molecular radical. By use of a DFT methodology, the minimum energy reaction pathways, and corresponding reaction barriers, of the above reactions over Zr(001), Nb(110), Mo(110), Tc(001), Ru(001), Rh(111), and Pd(111) have been determined. By consideration of the chemical potential and chemical hardness of the surface metal atoms, and the principle of electronegativity equalization, it is found that the charge transferred to the NH radical during the process of dissociative adsorption correlates very well with that determined by Mulliken population analysis. Furthermore, it is found that the stability of the NH/surface transition state complex relates directly to this charge transfer and that the trend in transition state stability predicted by a HSAB; treatment correlates very strongly with that determined by DFT calculations. With regards to N hydrogenation, we find that during the course of the reaction, H loses cohesion to the surface, as it must migrate from a 3-fold hollow site to either a bridge or top site, to react with N. Partial density of states (PDOS) and Mulliken population analysis reveal that this loss of bonding is accompanied by charge transfer from H to the surface metal atoms. Moreover, by simple modeling, we show that the reaction barriers are directly proportional to this mandatory charge transfer. Indeed, it is found that the reaction barriers correlate very well with the electrophilicity index of the metal atoms.

The temperature dependence of the adsorption of NO on Pt{211}: A RAIRS and DFT investigation

RAIRS experiments have been performed to investigate the adsorption of NO on Pt{211}. Results show that adsorption is complex and strongly temperature dependent. At 307 K, three bands are seen at saturation with frequencies of 1801, 1609, and 1576 cm(-1). However, at 120 K only two bands, at 1688 and 1620 cm(-1), are observed. To help with the assignment of these vibrational bands, DFT calculations were also performed. The calculations show that a bridged NO species, bonded to the step edge, is the most stable species on the surface and gives rise to the band observed at 1610-1620 cm(-1). The calculations also suggest that the temperature dependence of NO adsorption on Ptf{211} can be assigned to NO dissociation which occurs at room temperature but not at 120 K. In particular, the RAIRS band observed at 1801 cm(-1), which is observed on adsorption at 307 K but not at 120 K, is tentatively assigned to the formation of an O-NO complex. This species forms when a NO molecule bonds on top of an O atom, which results from the dissociation of NO on the Pt{211} surface at room temperature.

Resonance Raman and DFT studies of tetra-tert-butyl porphine: Assignment of strongly enhanced distortion modes in a ruffled porphyrin

The free-base form of tetra-tert-butyl porphine (TtBP), which has extremely bulky meso substituents, is severely distorted from planarity, with a ruffling angle of 65.5degrees. The resonance Raman spectrum of TtBP (lambda(ex) = 457.9 nm) and its d(2), d(8), and d(10) isotopomers have been recorded, and while the spectra show high-frequency bands similar to those observed for planar meso-substituted porphyrins, there are several additional intense bands in the low-frequency region. Density functional calculations at the B3-LYP/6-31G(d) level were carried out for all four isotopomers, and calculated frequencies were scaled using a single factor of 0.98. The single factor scaling approach was validated on free base porphine where the RMS error was found to be 14.9 cm(-1). All the assigned bands in the high-frequency (> 1000 cm(-1)) region of TtBP were found to be due to vibrations similar in character to the in-plane skeletal modes of conventional planar porphyrins. In the low-frequency region, two of the bands, assigned as nu(8) (ca. 330 cm(-1)) and nu(16) (ca. 540 cm(-1)), are also found in planar porphyrins such as tetra-phenyl porphine (TPP) and tetra-iso-propyl porphine (IPP). Of the remaining three very strong bands, the lowest frequency band was assigned as gamma(12) (pyr swivel, obsd 415 cm(-1), calcd 407 cm(-1) in do). The next band, observed at 589 cm-1 in the do compound (calcd 583 cm(-1)), was assigned as a mode whose composition is a mixture of modes that were previously labeled gamma(13) (gamma(CmCaHmCa)) andy gamma(11) (pyr fold(asym)) in NiOEP. The final strong band, observed at 744 cm(-1) (calcd 746 cm(-1)), was assigned to a mode whose composition is again a mixture of gamma(11) and gamma(13), although here it is gamma(11) rather than gamma(13) which predominates. These bands have characters and positions similar to those of three of the four porphyrin ring-based, weak bands that have previously been observed for NiTPP. In addition there are several weaker bands in the TtBP spectra that are also

Analysis procedure for calculation of electron energy distribution functions from incoherent Thomson scattering spectra

Incoherent Thomson scattering (ITS) provides a nonintrusive diagnostic for the determination of one-dimensional (1D) electron velocity distribution in plasmas. When the ITS spectrum is Gaussian its interpretation as a three-dimensional (3D) Maxwellian velocity distribution is straightforward. For more complex ITS line shapes derivation of the corresponding 3D velocity distribution and electron energy probability distribution function is more difficult. This article reviews current techniques and proposes an approach to making the transformation between a 1D velocity distribution and the corresponding 3D energy distribution. Previous approaches have either transformed the ITS spectra directly from a 1D distribution to a 3D or fitted two Gaussians assuming a Maxwellian or bi-Maxwellian distribution. Here, the measured ITS spectrum transformed into a 1D velocity distribution and the probability of finding a particle with speed within 0 and given value v is calculated. The differentiation of this probability function is shown to be the normalized electron velocity distribution function. (C) 2003 American Institute of Physics.

Breit-Pauli R-matrix calculation of fine-structure effective collision strengths for the electron impact excitation of Mg V

Context. Electron-impact excitation collision strengths are required for the analysis and interpretation of stellar observations.
Aims. This calculation aims to provide effective collision strengths for the Mg V ion for a larger number of transitions and for a greater temperature range than previously available, using collision strength data that include contributions from resonances.
Methods. A 19-state Breit-Pauli R-matrix calculation was performed. The target states are represented by configuration interaction wavefunctions and consist of the 19 lowest LS states, having configurations 2s22p4, 2s2p5, 2p6, 2s22p33s, and 2s22p33p. These target states give rise to 37 fine-structure levels and 666 possible transitions. The effective collision strengths were calculated by averaging the electron collision strengths over a Maxwellian distribution of electron velocities.
Results. The non-zero effective collision strengths for transitions between the fine-structure levels are given for electron temperatures in the range = 3.0 - 7.0. Data for transitions among the 5 fine-structure levels arising from the 2s22p4 ground state configurations, seen in the UV range, are discussed in the paper, along with transitions in the EUV range – transitions from the ground state 3P levels to 2s2p5?3P levels. The 2s22p4?1D–2s2p5?1P transition is also noted. Data for the remaining transitions are available at the CDS.

Electron-impact excitation of Cr II A theoretical calculation of collision and effective collision strengths for forbidden transitions

Context. Absorption or emission lines of Cr II are observed in a wide variety of astrophysical spectra and accurate atomic data are urgently needed to interpret these lines. Many of these data are impossible to measure experimentally and a full theoretical treatment is the only means by which these data can be obtained.

Aims. In this paper, we present collision strengths and effective collision strengths for electron-impact excitation of Cr II for forbidden transitions among the lowest-lying 74 fine-structure levels. Effective collision strengths have been computed for 18 individual electron temperatures of astrophysical importance, ranging from 2000-100 000 K.

Methods. The parallel suite of R-matrix packages, RMATRX II, which has recently been extended to allow for the inclusion of relativistic effects, has been used in the present work to compute the collision strengths and effective collision strengths for electron-impact excitation of Cr II. We concentrate in this publication on low-lying forbidden lines among the lowest 74 jj fine-structure levels with configurations 3d(5) and 3d(4)4s, although atomic data has been evaluated for all 39 060 transitions among the 280 jj levels of configurations 3d(5), 3d(4)4s and 3d(4)4p. This work constitutes the largest evaluation ever performed for this ion involving 1932 coupled channels.

Results. Collision and effective collision strengths are presented for all transitions among the lowest 74 J pi states of Cr II and comparisons made with the work of Bautista et al. (2009). While the effective collision strengths agree well for some transitions, significant discrepancies exist for others. We believe that the present atomic data represents the most accurate, most sophisticated and most complete data set for electron-impact excitation of Cr II and we would recommend them to astrophysicists and plasma physicists in their application work. We would expect that the effective collision strengths presented for the important low-lying forbidden lines are accurate to within 15%.

MEASUREMENT AND CALCULATION OF ABSOLUTE SINGLE AND MULTIPLE CHARGE EXCHANGE CROSS SECTIONS FOR Feq+ IONS IMPACTING H2O

Charge exchange (CE) plays a fundamental role in the collisions of solar- and stellar-wind ions with lunar and planetary exospheres, comets, and circumstellar clouds. Reported herein are absolute cross sections for single, double, triple, and quadruple CE of Feq+ (q = 5-13) ions with H2O at a collision energy of 7q keV. One measured value of the pentuple CE is also given for Fe9+ ions. An electron cyclotron resonance ion source is used to provide currents of the highly charged Fe ions. Absolute data are derived from knowledge of the target gas pressure, target path length, and incident and charge-exchanged ion currents. Experimental cross sections are compared with new results of the n-electron classical trajectory Monte Carlo approximation. The radiative and non-radiative cascades following electron transfers are approximated using scaled hydrogenic transition probabilities and scaled Auger rates. Also given are estimates of cross sections for single capture, and multiple capture followed by autoionization, as derived from the extended overbarrier model. These estimates are based on new theoretical calculations of the vertical ionization potentials of H2O up to H2O10+.