Protons colliding with crystalline ice: proton reflection and collision induced water desorption at low incidence energies

We present results of classical trajectory (CT) calculations on the sticking of protons to the basal plane (0001) face of crystalline ice, for normal incidence at a surface temperature (Ts) of 80 K. The calculations were performed for moderately low incidence energies (Ei) ranging from 0.05 to 4.0 eV. Surprisingly, significant reflection is predicted at low values of Ei (< or = 0.2 eV) due to repulsive electrostatic interactions between the incident proton and the surface water molecules with one of their H-atoms pointing upward toward the gas phase. The sticking probability increases with Ei and converges to unity for Ei > or = 0.8 eV. In the case of sticking, the proton is trapped in the ice forming a Zundel complex (H5O2+), with an average binding energy of 9.9 eV with a standard deviation of 0.5 eV, independent of the value of Ei. In nearly all sticking trajectories, the proton is implanted into the ice surface, with a penetration depth that increases with Ei. The strong interaction with the neighboring water molecules leads to a local rupture of the hydrogen bonding network, resulting in collision induced desorption of water (puffing), a process that occurs with significant probability even at the lowest Ei considered. The probability of water desorption increases with Ei. In nearly all trajectories in which water desorption occurs, a single three-coordinated water molecule is desorbed from the topmost monolayer.     

Experimental and theoretical investigation of HC5N adsorption on amorphous ice surface: simulation of the interstellar chemistry

HC 5N adsorbed on amorphous water ice at 10 K presents an interaction with the ice surface and induces the restructuring of the ice amorphous bulk. Warming up the sample induces the HC 5N desorption from the H 2O ice film, between 120 and 160 K, and the associated desorption energy is 90 kJ/mol. This value is in good agreement with that calculated E d (80 kJ/mol) and gives evidence that the amorphous ice surface is essentially dynamic. From theoretical calculations, it is shown that the HC 5N moiety presents a curvature and is no more linear and stabilized by two strong N...H bonds (2.09 and 2.29 A) and one H...O bond (1.84 A).     

Interstellar ice surface site modification induced by dicyanoacetylene adsorption

Dicyanoacetylene adsorbed on amorphous ice water at 10 K presents an interaction with the dangling H site and induces a s(4) adsorption site formation due to the restructuring of the ice bulk. Warming up the sample provokes the dicyanoacetylene desorption from the H(2)O ice film, which could be due to the beginning of the ice crystallization process. The desorption activation energy measured by temperature-programmed desorption (E(d) = 42 +/- 5 kJ x mol(-1)) is in good agreement with that calculated (E(d) = 46 kJ x mol(-1)) and gives evidence of a hydrogen-bonded adsorbed state on amorphous ice films.     

Determination of the sticking coefficient and scattering dynamics of water on ice using molecular beam techniques

The sticking coefficient for D(2)O impinging on crystalline D(2)O ice was determined for incident translational energies between 0.3 and 0.7 eV and for H(2)O on crystalline H(2)O ice at 0.3 eV. These experiments were done using directed molecular beams, allowing for precise control of the incident angle and energy. Experiments were also performed to measure the intensity and energy of the scattered molecules as a function of scattering angle. These results show that the sticking coefficient was near unity, slightly increasing with decreasing incident energy. However, even at the lowest incident energy, some D(2)O did not stick and was scattered from the ice surface. We observe under these conditions that the sticking probability asymptotically approaches but does not reach unity for water sticking on water ice. We also present evidence that the scattered fraction is consistent with a binary collision; the molecules are scattered promptly. These results are especially relevant for condensation processes occurring under nonequilibrium conditions, such as those found in astrophysical systems.     

The CO formation reaction pathway in steam methane reforming by rhodium

Three different pathways toward CO formation from adsorbed CH and O are compared by quantum-chemical density functional theory (DFT) calculations for planar and stepped Rh surfaces. The conventional pathway competes with the pathway involving a formyl (CHO) species. This holds for both types of surfaces. The barrier for carbon-oxygen bond formation for the planar surface (180 kJ/mol) is substantially higher than that for the stepped surface (90 kJ/mol). The reaction path through intermediate formyl formation competes with direct formation of CO from recombination via adsorbed C and O atoms. Calculations are used as a basis for the analysis of the overall kinetics of the methane steam reforming reaction as a function of the particle size and the metal.     

Using effusive molecular beams and microcanonical unimolecular rate theory to characterize CH4 dissociation on Pt(111)

The dissociative sticking coefficient for CH4 on Pt(111) has been measured as a function of both gas temperature (Tg) and surface temperature (Ts) using effusive molecular beam and angle-integrated ambient gas dosing methods. The experimental results are used to optimize the three parameters of a microcanonical unimolecular rate theory (MURT) model of the reactive system. The MURT calculations allow us to extract transition state properties from the data as well as to compare our data directly to other molecular beam and thermal equilibrium sticking measurements. We find a threshold energy for dissociation of E0 = 52.5 +/- 3.5 kJ mol(-1). Furthermore, the MURT with an optimized parameter set provides for a predictive understanding of the kinetics of this C-H bond activation reaction, that is, it allows us to predict the dissociative sticking coefficient of CH4 on Pt(111) for any combination of Ts and Tg even if the two are not equal to one another, indeed, the distribution of molecular energy need not even be thermal. Comparison of our results to those from recent thermal equilibrium catalysis studies on CH4 reforming over Pt nanoclusters ( approximately 2 nm diam) dispersed on oxide substrates indicates that the reactivity of Pt(111) exceeds that of the Pt nanocatalysts by several orders of magnitude.     

Influence of submonolayer sodium adsorption on the photoemission of the Cu(111)/water ice surface

Photoemission from an ice film deposited on Cu(111) as a function of thickness has been observed in the presence and absence of sodium atoms at the surface-vacuum interface. For either adsorbate alone and photon energies below 4 eV, two-photon photoemission from the Cu(111) substrate dominates. The Cu(111) photoelectron spectrum is perturbed by low coverages of Na, and its intensity is strongly attenuated by a few monolayers of ice. For a low density amorphous ice film, strong charging effects are observed. For ice films annealed to yield either the dense amorphous or crystalline phase, this effect is absent. Deposition of only 0.02 monolayer of Na leads to a dramatic decrease in the threshold for photoemission to 2.3+/-0.2 eV. Thus, photoelectrons are generated by visible radiation in a one-photon process with a cross section that exceeds 10(-18) cm(2). The initial state for the photoemission is identified as a metastable surface trapped electron, which decays thermally with an activation energy of 10+/-2 kJ mol(-1). Quantum calculations are described which support this model and show that the Na atom is accommodated in the first layer of the ice surface.     

Adsorption of acetaldehyde on ice as seen from computer simulation and infrared spectroscopy measurements

Detailed investigation of the adsorption of acetaldehyde on I(h) ice is performed under tropospheric conditions by means of grand canonical Monte Carlo computer simulations and compared to infrared spectroscopy measurements. The experimental and simulation results are in a clear accordance with each other. The simulations indicate that the adsorption process follows Langmuir behavior in the entire pressure range of the vapor phase of acetaldehyde. Further, it was found that the adsorption layer is strictly monomolecular, and the adsorbed acetaldehyde molecules are bound to the ice surface by only one hydrogen bond, typically formed with the dangling H atoms at the ice surface, in agreement with the experimental results. Besides this hydrogen bonding, at high surface coverages dipolar attraction between neighboring acetaldehyde molecules also contributes considerably to the energy gain of the adsorption. The acetaldehyde molecules adopt strongly tilted orientations relative to the ice surface, the tilt angle being scattered between 50° and 90° (i.e., perpendicular orientation). The range of the preferred tilt angles narrows, and the preference for perpendicular orientation becomes stronger upon saturation of the adsorption layer. The CH(3) group of the acetaldehyde molecules points as straight away from the ice surface within the constraint imposed by the tilt angle adopted by the molecule as possible. The heat of adsorption at infinitely low coverage is found to be -36 ± 2 kJ/mol from the infrared spectroscopy measurement, which is in excellent agreement with the computer simulation value of -34.1 kJ/mol.     

Pressure amorphized ices--an atomistic perspective

We offer our viewpoint on the nature of amorphous ices produced by pressurization of crystalline ice Ih and the inter-relationship between them from an atomistic perspective. We argue that the transformation of high density amorphous (HDA) ice from crystalline ice is due to a mechanical process arising from the instability of the ice Ih structure. The densification of HDA upon thermal annealing under pressure is a relaxation process. The conversion of the densified amorphous ice to a lower density form (LDA) upon the release of pressure can be attributed to a similar process. It is speculated that amorphous ices are metastable frustrated structures due to the large activation barriers associated with proton reorientation in the formation of the underlying stable crystalline ice polymorphs.     

Sticking coefficient and processing of water vapor on organic-coated nanoaerosols

Organic-based aerosols may play an unexpectedly important role in the atmospheric processing of water vapor. In this paper, we report results of a molecular dynamics (MD) study on the sticking coefficient of water vapor on a coated particle (water droplet coated with fatty acid of radius approximately 4 nm). The sticking coefficient was found to be almost a constant (11-16%) for incident speeds around the most probable speed as opposed to 100% for water vapor incident on a pure water droplet. The sticking coefficient was found to increase with the size of the water cluster (water-Nmer) impinging the surface of the particle and was seen to approach 1 for impinging water clusters consisting of 10 molecules or more. We also computed the average energy transferred per collision for monomers impinging the surface of the coated particle and found the value to be 4.4581 kJ/mol. Most important perhaps was the fact that despite the lower sticking coefficient, the equilibrium vapor pressure of water over these inverted micelles was considerably lower due to the surface tension effects of the fatty acid layer. As such, these coated particles act as effective substrates for water vapor condensation and may play a role as cloud condensation nuclei.     

Interaction of atomic and molecular deuterium with a nonporous amorphous water ice surface between 8 and 30 K

Molecular and atomic interactions of hydrogen on dust grains covered with ice at low temperatures are key mechanisms for star formation and chemistry in dark interstellar clouds. We have experimentally studied the interaction of atomic and molecular deuterium on nonporous amorphous water ice surfaces between 8 and 30 K, in conditions compatible with an extrapolation to an astrophysical context. The adsorption energy of D(2) presents a wide distribution, as already observed on porous water ice surfaces. At low coverage, the sticking coefficient of D(2) increases linearly with the number of deuterium molecules already adsorbed on the surface. Recombination of atomic D occurs via a prompt reaction that releases molecules into the gas phase. Part of the newly formed molecules are in vibrationally excited states (v=1-7). The atomic recombination efficiency increases with the presence of D(2) molecules already adsorbed on the water ice, probably because these increase the sticking coefficient of the atoms, as in the case of incident D(2). We have measured the atomic recombination efficiency in the presence of already absorbed D(2), as it is expected to occur in the interstellar medium. The recombination efficiency decreases rapidly with increasing temperature and is zero at 13 K. This allows us to estimate an upper limit to the value of the atom adsorption energy E(a) approximately 29 meV, in agreement with previous calculations.     

On the dynamics of H2 adsorption on the Pt(111) surface

The dynamics and kinetics of the dissociation of hydrogen over the hexagonal close packed platinum (Pt(111)) surface are investigated using Car–Parrinello molecular dynamics and static density functional theory calculations of the potential energy surfaces. The calculations model the reference energy‐resolved molecular beam experiments, considering the degrees of freedom of the catalytic surface. Two‐dimensional potential energy surfaces above the main sites on Pt(111) are determined. Combined with Car–Parrinello trajectories, they confirm the dissociative adsorption of H₂ as the only adsorption pathway on this surface at H₂ incindence energies above 5 kJ/mol. A direct determination of energy‐resolved sticking coefficients from molecular dynamics is also performed, showing an excellent agreement with the experimental data at incidence energies in the 5–30 kJ/mol range. Application of dispersion corrections does not lead to an improvement in the prediction of the H₂ sticking coefficient. The adsorption reaction rate obtained from the calculated sticking coefficients is consistent with experimentally derived literature values.     

Density functional theory study of CO catalytic oxidation on Co_2B_2/TiO_2（110） surface

Titanium dioxide with CoB amorphous alloys nanoparticles deposited on the surface is known to exhibit higher catalytic activity than the CoB amorphous. A study of the structure of such system is necessary to understand this effect. A quantum chemical study of Co2B2 on the TiO2 (110) surface was studied using periodic slab model within the framework of density functional theory (DFT). The results of geometry optimization indicated that the most stable model of adsorption was Co2B2 cluster adsorbed on the hollow site of TiO2.The adsorption energy calculated for Co2B2 on the hollow site was 439.3 kJ/mol.The adsorption of CO and O2 was further studied and the results indicated that CO and O2 are preferred to adsorb on the Co2 site. Co-adsorption of CO and O2 shows that Co2B2/TiO2 is a good catalyst for the oxidation of CO to carbon dioxide in the presence of oxygen.     

Dynamically biased RRKM model of activated gas-surface reactivity: vibrational efficacy and rotation as a spectator in the dissociative chemisorption of CH4 on Pt(111)

The reactivity of CH(4) impinging on a Pt(111) surface was examined using a precursor-mediated microcanonical trapping model of dissociative chemisorption wherein the effects of rotational and vibrational energy could be explored. Dissociative sticking coefficients for a diverse range of non-equilibrium effusive beam, supersonic beam, and eigenstate-resolved experiments were simulated and an average relative discrepancy between theory and experiment of better than 50% was achieved by treating molecular rotations and translation parallel to the surface as spectator degrees of freedom, and introducing a dynamically-biased vibrational efficacy. The model parameters are {E(0) = 57.9 kJ mol(-1), s = 2, η(v) = 0.40} where E(0) is the apparent threshold energy for reaction, s is the number of surface oscillators participating in energy exchange within each gas-surface collision complex formed, and η(v) is the mean vibrational efficacy for reaction relative to normal translational energy which figures in the assembly of the active exchangeable energy which is available to surmount the activation barrier to dissociative chemisorption. GGA-DFT electronic structure calculations provided vibrational frequencies for the transition state for dissociative chemisorption. The asymmetry of the rotational state populations in supersonic and effusive molecular beam experiments allowed kinetic analysis to establish that taking rotation as a spectator degree of freedom is a good approximation. Surface phonons, rather than the incident molecules, are calculated to play the dominant role in supplying the energy required to overcome the activation barrier for dissociative chemisorption under the thermal equilibrium conditions relevant to high pressure catalysis. Over the temperature range 300 K ≤T≤ 1000 K, the thermal dissociative sticking coefficient is predicted to be well described by S(T) = S(0) exp(-E(a)/RT) where S(0) = 0.62 and E(a) = 62.6 kJ mol(-1).     

Adsorption and activation of CO over flat and stepped Co surfaces: a first principles analysis

The adsorption and activation of CO over flat Co{0001}, corrugated Co{1120}, and stepped Co{1012} and Co{1124} surfaces have been analyzed using periodic density functional theory calculations. CO strongly chemisorbs on all these surfaces but does not show a strong dependence on the surface structure. The calculated structure of adsorbed CO on Co{0001} at 1/3 monolayer (ML) of coverage was found to be in good agreement with the experiment. The barrier for CO dissociation over Co{0001} was found to decrease with decreasing CO coverage, taking on a value of 232 kJ/mol at 1/4 ML and 218 kJ/mol at 1/9 ML. The presence of the "zigzag" channel on Co{1120} enhances the reactivity slightly by reducing the barrier for CO dissociation to 195 kJ/mol. In contrast, the stepped Co{1012} and Co{1124} surfaces are much more active than the flat and corrugated surfaces. Both stepped surfaces provide direct channels for CO dissociation that do not have barriers with respect to gas-phase CO. In general the activation barriers lower as the reaction energies become more exothermic. Reconstruction of the step edges that occur in the product state, however, prevents a linear correlation between the reaction energy and the activation energy.     

Cyanoacetylene adsorption on amorphous and crystalline water ice films: investigation through matrix isolation and quantum study.

The structure and energy properties of the 1:1 complexes formed between cyanoacetylene and H2O (D2O) are investigated using FT-IR matrix isolation spectroscopy and ab initio calculations at the MP2/ 6-31G(d,p) level. Cyanoacetylene adsorption and desorption on amorphous ice film are monitored by FT-IR using the temperature-programmed desorption method. In an argon matrix, two types of 1:1 complexes are observed. The first one corresponds to the NH structure, which involves a hydrogen bond with the terminal nitrogen of cyanoacetylene. The second corresponds to the HO form, which involves a hydrogen bond from the cyanoacetylene to the oxygen of water. This last complex is the more stable (DeltaE = -8.1 kJ/mol.). As obtained in argon matrixes, two kinds of adsorption site are observed between HC3N and ice. The first one, stable between 25 and 45 K is characterized by a nu(OH) shift similar to the one observed in matrix for the NH complex. The second, stable at higher temperatures (between 45 and 110 K), corresponds to an interaction with the dangling oxygen site of ice and is similar to the HO complex observed in matrix. From theoretical calculations (DFT method combined with a plane wave basis set and ultrasoft pseudopotentials), it is shown that, for this adsorption site, the HC3N moiety is flattened on the ice surface and stabilized by a long-distance interaction ( approximately 3 A) between one dangling OH and the pi system of the C triple bond C triple bond. The HC3N desorption occurs between 110 and 140 K, and the associated desorption energy is 39 kJ/mol. This value is in good agreement with the first principle calculation based on density functional theory and ultrasoft pseudopotentials (34 kJ/mol). These calculations confirm the electrostatic nature of the interaction forces. A small amount of cyanoacetylene is incorporated into the bulk and desorbs at the onset of the ice crystallization near 145 K. In these two kinds of experiments, HC3N acts as both an electrophilic and a nucleophilic molecule.     

First-principles study of the adsorption of formaldehyde on the clean and atomic oxygen covered Cu(1 1 1) surface

The interaction of formaldehyde with the clean and atomic oxygen-covered Cu(1 1 1) surfaces has been studied by means of cluster model density functional calculations in which Cu₂₂(14,8) is used to represent the perfect Cu(1 1 1) surface. The calculations point towards a η¹-H₂CO---O orientation with the oxygen atom almost on top of a copper surface atom. The formaldehyde adsorption energy is of 22–25 kJ/mol and the internal geometry of adsorbed formaldehyde is almost identical to that of the molecule in the gas-phase. The C---O bond is almost parallel to the surface and the conformation with the molecular plane normal to the surface is slightly preferred to the conformation with the molecular plane nearly parallel to the surface. A Cu₂₂---O model where atomic oxygen is adsorbed on a fcc hollow site was used to study the co-adsorption and reaction of formaldehyde with atomic oxygen. Oxygen co-adsorption has a dramatic effect on the formaldehyde adsorption energy which is increased by 50%. The calculated energy barrier for the formation of the dioxymethylene intermediate species through the H₂CO+O→H₂CO₂ reaction is of 36 kJ/mol.     

Effusive molecular beam study of C2H6 dissociation on Pt(111)

The dissociative sticking coefficient for C2H6 on Pt(111) has been measured as a function of both gas temperature (Tg) and surface temperature (Ts) using effusive molecular beam and angle-integrated ambient gas dosing methods. A microcanonical unimolecular rate theory (MURT) model of the reactive system is used to extract transition state properties from the data as well as to compare our data directly with supersonic molecular beam and thermal equilibrium sticking measurements. We report for the first time the threshold energy for dissociation, E0 = 26.5 +/- 3 kJ mol(-1). This value is only weakly dependent on the other two parameters of the model. A strong surface temperature dependence in the initial sticking coefficient is observed; however, the relatively weak dependence on gas temperature indicates some combination of the following (i) not all molecular excitations are contributing equally to the enhancement of sticking, (ii) that strong entropic effects in the dissociative transition state are leading to unusually high vibrational frequencies in the transition state, and (iii) energy transfer from gas-phase rovibrational modes to the surface is surprisingly efficient. In other words, it appears that vibrational mode-specific behavior and/or molecular rotations may play stronger roles in the dissociative adsorption of C2H6 than they do for CH4. The MURT with an optimized parameter set provides for a predictive understanding of the kinetics of this C-H bond activation reaction, that is, it allows us to predict the dissociative sticking coefficient of C2H6 on Pt(111) for any combination of Ts and Tg even if the two are not equal to one another.     

Calcium adsorption on MgO(100): energetics, structure, and role of defects

The adsorption of Ca on the MgO(100) surface at 300 K has been studied using microcalorimetry, in combination with LEED, AES, ISS, work function, sticking probability measurements, and density functional theory (DFT) calculations. The MgO(100) thin films (approximately 4 nm thick) were grown epitaxially on a 1 microm thick Mo(100) single-crystal. The sticking probability of Ca on MgO(100) at 300 K is unity. On the basis of AES and ISS measurements, it was determined that Ca grows mainly as 3D particles on the MgO(100) surface with a density of approximately 1 x 10(12) islands/cm2. Ca adsorbs initially at defect sites with a very high heat of adsorption (approximately 410 kJ/mol). DFT calculations attribute this high initial heat to Ca binding to kink sites (376 kJ/mol), step sites (205 kJ/mol), and lower concentrations of stronger binding sites. The heat of adsorption decreases rapidly with coverage, reaching a minimum of 162 kJ/mol at approximately 0.3 ML, where Ca is mainly adding to small 3D Ca clusters. Afterward, it increases to the value of bulk Ca heat of sublimation (178 kJ/mol) at approximately 1.2 ML, attributed to the increase in stability with increasing Ca particle size. A 1.0 eV decrease of the work function with Ca coverage from 0 to 0.3 ML indicates that Ca adsorbed at defects is cationic, in agreement with calculations showing that Ca donates electron density to the MgO. Light ion sputtering of the MgO(100) surface generates point defects, but these do not change the heat of adsorption versus coverage, implying that they do not nucleate Ca particles. Oxygen vacancies are a likely candidate; DFT calculations show that F and F+ center vacancies bind Ca more weakly than terrace sites. More extensive sputtering creates extended defects (such as steps and kinks) that adsorb Ca with heats of adsorption up to approximately 400 kJ/mol, similar to that at the intrinsic defect sites.     

Microcanonical transition state theory for activated gas-surface reaction dynamics: application to H2/CU(111) with rotation as a spectator

A microcanonical unimolecular rate theory (MURT) model incorporating quantized surface vibrations and Rice-Ramsperger-Kassel-Marcus rate constants is applied to a benchmark system for gas-surface reaction dynamics, the activated dissociative chemisorption and associative desorption of hydrogen on Cu(111). Both molecular translation parallel to the surface and rotation are treated as spectator degrees of freedom. MURT analysis of diverse experiments indicates that one surface oscillator participates in the dissociative transition state and that the threshold energy for H2 dissociation on Cu(111) is E0 = 62 kJ/mol. The spectator approximation for rotation holds well at thermally accessible rotational energies (i.e., for Er less than approximately 40 kJ/mol). Over the temperature range from 300 to 1000 K, the calculated thermal dissociative sticking coefficient is ST = S0 exp(-Ea/kBT) where S0 = 1.57 and Ea = 62.9 kJ/mol. The sigmoid shape of rovibrational eigenstate-resolved dissociative sticking coefficients as a function of normal translational energy is shown to derive from an averaging of the microcanonical sticking coefficient, with threshold energy E0, over the thermal surface oscillator distribution of the gas-surface collision complexes. Given that H2/Cu(111) is one of the most dynamically biased of gas-surface reactive systems, the simple statistical MURT model simulates and broadly rationalizes the H2/Cu(111) reactive behavior with remarkable fidelity.