Computational investigation of NH3 adsorption and dehydrogenation on a W-modified Fe(111) surface.

Hydrogen gas will play an important role in the future since it could be a replacement for gasoline, heating oil, natural gas, and other fuels. In previous reports ammonia (NH3), which has a high hydrogen content, provides a promising mode for the transferring and storing of hydrogen for its on-site generation. Therefore, the dehydrogenation of NH3 on a metal surface has been studied widely in the last few decades. In our study, we employed monolayer tungsten metal to modify the Fe(111) surface, denoted as W@Fe(111), and calculated the adsorption and dehydrogenation behaviors of NH3 on W@Fe(111) surface via first-principles calculations based on density functional theory (DFT). The three adsorption sites of the surface, top (T), 3-fold-shallow (S), and 3-fold-deep (D) were considered. The most stable structure of the NHx (x = 0-3) species on the surface of W@Fe(111) have been predicted. The calculated activation energies for NHx (x = 1-3) dehydrogenations are 19.29 kcal mol(-1) (for H2N-H bond activation), 29.17 kcal mol(-1) (for HN-H bond activation) and 27.94 kcal mol(-1) (for N-H bond activation), and the entire process is exothermic by 33.05 kcal mol(-1). To gain detailed knowledge of the catalytic processes of the NH3 molecule on the W@Fe(111) surface, the physical insights between the adsorbate/substrate interaction and interface morphology were subjected to a detailed electronic analysis.

Reaction mechanisms and kinetics of the iminovinylidene radical with NO: ab initio study.

The nitric oxide (NO) is a notorious compound for polluting environment. Recent year, removing nitric oxide from the atmosphere becomes a focus of the investigation. In our work, we study the iminovinylidene (HNCC) radical reacted with NO molecule. The mechanism and kinetic for reaction of the HNCC radical with the NO molecule is investigated via considering the possible channels of the N and O atoms of NO attacking the N and C atoms of the HNCC based on the high level ab initio molecular orbital calculations in conjunction with variational TST and RRKM calculations. The species involved have been optimized at the B3LYP/6-311++G(3df,2p) level and their single-point energies are refined by the CCSD(T)/aug-cc-PVQZ//B3LYP/6-311++G(3df,2p) method. The calculated potential energy surfaces indicated that energetically the most favorable channel for the HNCC + NO reaction was predicted to be the formation of HNC+CNO (P8) product via the addition reaction of the C atom of HNCC radical and the N atom of NO with the head to head orientation. To rationalize the scenario of the calculated results, we also employ the Fukui functions and HSAB theory to seek for a possible explanation. In addition, the reaction rate constants were calculated using VariFlex code, and the results show that the total rate coefficient, ktotal, at Ar pressure 760 Torr can be represented with an equation: ktotal = 6.433 × 10(-11) T (0.100) exp(0.275 kcal mol(-1)/RT) at T = 298-3000 K, in units of cm(3) molecule(-1) s(-1).

Quasiclassical trajectory calculations for Li(2(2)P(J)) + H2 → LiH(X(1)Σ+) + H: influence by vibrational excitation and translational energy.

Ab initio potential energy surfaces and the corresponding analytical energy functions of the ground 1A' and excited 2A' states for the Li(2(2)P) plus H(2) reaction are constructed. Quasiclassical trajectory calculations on the fitted energy functions are performed to characterize the reactions of Li(2(2)P) with H(2)(v = 0, j = 1) and H(2)(v = 1, j = 1) as well as the reaction when the vibrational energy is replaced by collision energy. For simplicity, the transition probability is assumed to be unity when the trajectories go through the crossing seam region and change to the lower surface. The calculated rotational distributions of LiH(v = 0) for both H(2)(v = 0, j = 1) and H(2)(v = 1, j = 1) reactions are single-peaked with the maximum population at j' = 7, consistent with the previous observation. The vibrational excitation of H(2)(v = 1) may enhance the reaction cross section of LiH(v' = 0) by about 200 times, as compared to a result of 93-107 reported in the experimental measurements. In contrast, the enhancement is 3.1, if the same amount of energy is deposited in the translational states. This endothermic reaction can be considered as an analog of late barrier. According to the trajectory analysis, the vibrational excitation enlarges the H-H distance in the entrance channel to facilitate the reaction, but the excess energy may not open up additional reaction configuration.

Photodissociation of (ICN)(2) van der Waals dimer using velocity imaging technique.

Photodissociation of (ICN)(2) dimer from 265 to 270 nm are studied using time-of-flight mass spectrometry combined with velocity imaging technique. Both I(+) and I(2) (+) ions are found in the mass spectra. The I(2) (+) ions result from (1+1) resonant ionization of the neutral I(2) fragment that is produced in the photodissociation of dimer, but not from dissociative ionization of (ICN)(2); i.e., (ICN)(2) (+)+hnu-->I(2) (+)+2CN. The dissociation channels of I(2) (+) leading to I(+) are all found with parallel character. The total kinetic energy distributions and anisotropy parameters of the I(+) channels produced by (ICN)(2) are almost the same as those from a neutral I(2) sample, thereby confirming that the I(2) fragments are obtained in cold state. With the aid of ab initio calculations, a plausible dissociation mechanism is proposed.

Probing the ignored elimination channel of Br2 in the 248 nm photodissociation of 1,1-dibromoethylene by cavity ring-down absorption spectroscopy.

In the photodissociation of 1,1-C(2)H(2)Br(2) at 248 nm, the Br(2) elimination channel is probed by using cavity ring-down absorption spectroscopy (CRDS). In terms of spectral simulation, the vibrational population ratio of Br(2)(v = 1)/Br(2)(v = 0) is found to be 0.55+/-0.05, which indicates that the Br(2) fragment is vibrationally hot. The rotational population is thermally equilibrated with a Boltzmann temperature of 349+/-38 K. According to ab initio potential energy calculations, the obtained fragments are anticipated to result primarily from photodissociation of the ground electronic state that undergoes 1) H migration followed by three-center elimination, and 2) isomerization forming either trans- or cis-1,2-C(2)H(2)Br(2) from which Br(2) is eliminated. RRKM calculations predict that the Br(2) dissociation rates through the ground singlet state prevail over those through the triplet state. Measurements of temperature and Ar pressure dependence are examined to support the proposed pathway via internal conversion. The quantum yield for the Br(2) elimination reaction is determined to be 0.07+/-0.04. This result is smaller than that obtained in 1,2-C(2)H(2)Br(2), probably because the dissociation rates are slowed in the isomerization stage.

Rotational and vibrational state distributions of NaH in the reactions of Na(4 (2)S,3 (2)D, and 6 (2)S) with H2: Insertion versus harpoon-type mechanisms.

By using a pump-probe technique, the nascent rotational and vibrational state distributions of NaH are obtained in the Na(4 (2)S,3 (2)D, and 6 (2)S) plus H(2) reactions. The rotational distributions for the Na(4 (2)S,3 (2)D) reactions yield a bimodal feature with a major component peaking at J=20-22, similar to that obtained previously in the 4 (2)P reaction, whereas the Na(6 (2)S) reaction gives rise to a distinct distribution with a much lower rotational temperature. The vibrational populations (v=0-4) for these 4 (2)S, 3 (2)D, and 6 (2)S reactions are characterized by corresponding temperatures of 1692+/-120, 819+/-35, and 5329+/-350 K. Due to a significant contribution of configurational mixing between different states with the same symmetry, the collision species initiated from the 4 (2)S and 3 (2)D states are anticipated to track along the entrance surface in a near C(2v) symmetry, then undergo nonadiabatic transition to the inner limb of the reactive 2A(') surface. In contrast, the reaction pathway for the Na(6 (2)S) state with a significantly reduced ionization energy is anticipated to follow a harpoon-type mechanism via a (near) collinear configuration. The increased atomic size of Na may hinder the insertion approach.

Reaction dynamics of Ca(4s3d 1D2)+CH4-->CaH(X 2Sigma+)+CH3: reaction pathway and energy disposal for the CaH product.

The reaction pathway for Ca(4s3d 1D2)+CH4-->CaH(X 2Sigma+)+CH3 has been investigated by using a pump-probe technique in combination with potential-energy surface (PES) calculations. The nascent product distributions of CaH have been characterized with Boltzmann rotational temperatures of 1013+/-102 and 834+/-70 K for the v=0 and 1 levels, respectively, and a Boltzmann vibrational temperature of 1313+/-173 K. The rotational and vibrational energy partitions in CaH have been estimated to be 461+/-45 and 252+/-15 cm(-1), respectively. According to the PES calculations, the pathway favors an insertion mechanism. Ca(3 1D2) approaches CH4 along an attractive potential surface in a C2v (or Cs) symmetry and then the collision complex undergoes nonadiabatic transition to the reactive ground-state surface. An Arrhenius plot shows a potential-energy requirement of 2695+/-149 cm(-1), which accounts for the endothermicity of 2930 cm(-1) for the reaction scheme. The Ca-C bond distance in the transition state structure is short enough to allow for tight orbital overlap between CaH and CH3. The strong coupling between the moieties renders the energy transfer sufficient from CaH into the CH3 radical. As compared to the Ca(4 1P1) reaction, the dissociation lifetime of the intermediate complex with less excess energy is prolonged so as to cause much less vibrational energy disposal into CaH.

Reaction pathway for the nonadiabatic reaction of Ca(4s3d 1D)+H2-->CaH(X 2Sigma+)+H.

The reaction pathway and the nascent CaH product distribution in the reaction Ca(4s3d (1)D)+H(2)-->CaH(X (2)Sigma(+))+H are obtained using a pump-probe technique. The Ca atom is first prepared in the 3 (1)D state by a two-photon absorption, and then in brief time delay the laser-induced fluorescence of the reaction product CaH is monitored. The CaH(v=0,1) distributions appear to be single peaked, as characterized by Boltzmann rotational temperature of 807+/-38 K (v=0) and 684+/-77 K (v=1). The vibrational population ratio of CaH(v=0)/CaH(v=1) is determined to be 3.3+/-0.1, while the v=2 population is not detectable. The fractions of the available energy partitioning into rotation, vibration, and translation are estimated to be 0.36+/-0.05, 0.28+/-0.04, and 0.36+/-0.05, respectively. With the aid of the potential energy surfaces calculations, the current reaction should favor a near C(2v) collision configuration. The temperature dependence measurement yields a positive slope, indicative of the reaction occurrence without any potential barrier. The colliding species are anticipated to follow an attractive 1B(2) (or 2A') surface and then transit nonadiabatically to the reactive ground state surface.

Influence of vibrational excitation on the nonadiabatic reactions of metal atoms with H2.

The reactions of alkaline earth metal atoms, Mg(3s3p 1P1) and Ca(4s4p 1P1), with H2(v = 1, j) are studied using a pump-probe technique combined with stimulated Raman pumping and coherent anti-Stokes Raman spectroscopy. For the Ca(4 1P1) case, the energy deposited in the v = 1 level enlarges the H2 bond distance to help facilitate the reaction without opening an additional pathway. For the Mg(3 1P1) case, the vibrational excitation of H2 leads to enhancement of the low rotational component of the rotational distribution and the MgH(v = 0)/MgH(v = 1) ratio. These results can be predicted with quasi-classical trajectory calculations and interpreted with a kinematic collision model.

Reaction pathway and potential barrier for the CaH product in the reaction of Ca(4s4p1P1) + H2 --> CaH(X2Sigma+) + H.

The nascent CaH product in the reaction Ca(4s4p1P1) + H2 --> CaH(X2Sigma+) + H is obtained using a pump-probe technique. The CaH(v = 0,1) distributions, with a population ratio of CaH(v = 0)/CaH(v = 1) = 2.7+/-0.2, may be characterized by low Boltzmann rotational temperature. According to Arrhenius theory, the temperature dependence measurement yields a potential barrier of 3820+/-480 cm(-1) for the current reaction. As a result of the potential energy surfaces (PES) calculations, the reaction pathway favors a Ca insertion into the H2 bond along a (near) C2v geometric approach. As the H2 bond is elongated, the configurational mixing between the orbital components of the 4p and nearby low-lying 3d state with the same symmetry makes significant the nonadiabatic transition between the 5A' and 2A' surface in the repulsive limbs. Therefore, the collision species are anticipated to track along the 5A' surface, then undergo nonadiabatic transition to the inner limb of the 2A' surface, and finally cross to the reactive 1A' surface. The observed energy barrier probably accounts for the energy requirement to surmount the repulsive hill in the entrance. The findings of the nascent CaH product distributions may be reasonably interpreted from the nature of the intermediate structure and lifetime after the 2A'-1A' surface transition. The distinct product distributions between the Ca(4 1P1) and Mg(3 1P1) reactions with H2 may also be realized with the aid of the PES calculations.