Sticking of CO to crystalline and amorphous ice surfaces

We present results of classical trajectory calculations on the sticking of hyperthermal CO to the basal plane (0001) face of crystalline ice Ih and to the surface of amorphous ice Ia. The calculations were performed for normal incidence at a surface temperature Ts = 90 K for ice Ia, and at Ts = 90 and 150 K for ice Ih. For both surfaces, the sticking probability can be fitted to a simple exponentially decaying function of the incidence energy, Ei: Ps = 1.0e(-Ei(kJ/mol)/90(kJ/mol)) at Ts = 90 K. The energy transfer from the impinging molecule to the crystalline and the amorphous surface is found to be quite efficient, in agreement with the results of molecular beam experiments on the scattering of the similar molecule, N2, from crystalline and amorphous ice. However, the energy transfer is less efficient for amorphous than for crystalline ice. Our calculations predict that the sticking probability decreases with Ts for CO scattering from crystalline ice, as the energy transfer from the impinging molecule to the warmer surfaces becomes less efficient. At high Ei (up to 193 kJ/mol), no surface penetration occurs in the case of crystalline ice. However, for CO colliding with the amorphous surface, a penetrating trajectory was observed to occur into a large water pore. The molecular dynamics calculations predict that the average potential energy of CO adsorbed to ice Ih is -10.1 +/- 0.2 and -8.4 +/- 0.2 kJ/mol for CO adsorbed to ice Ia. These values are in agreement with previous experimental and theoretical data. The distribution of the potential energy of CO adsorbed to ice Ia was found to be wider (with a standard deviation sigma of 2.4 kJ/mol) than that of CO interacting with ice Ih (sigma = 2.0 kJ/mol). In collisions with ice Ia, the CO molecules scatter at larger angles and over a wider distribution of angles than in collisions with ice Ih.     

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.     

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.     

Coverage-dependent sticking probability and desorption kinetics of water molecules on Rh(111)

The adsorption and desorption kinetics of water molecules on Rh(111) were investigated using temperature programed desorption (TPD). Water molecules on Rh(111) show coverage-dependent sticking probability; the initial sticking probability is estimated to be 0.46. In the desorption process, a dilute gaslike phase and two-dimensional islands of water coexist on the surface. Based on the model proposed by Kreuzer and Payne [Surf. Sci.200, L433 (1988)], the apparent fractional-order TPD spectra can be interpreted as first-order desorption from the coexistence of two phases on which the sticking probabilities are different. Based on this, the previous estimation of pre-exponential factors assuming half-order desorption [A. Beniya et al., J. Chem. Phys.125, 054717 (2006)] should be revised.     

Dewetting growth of crystalline water ice on a hydrogen saturated Rh(111) surface at 135 K

We investigated the water (D(2)O) adsorption at 135?K on a hydrogen pre-adsorbed Rh(111) surface using temperature programmed desorption and infrared reflection absorption spectroscopy (IRAS) in ultrahigh vacuum. With increasing the hydrogen coverage, the desorption temperature of water decreases. At the saturation coverage of hydrogen, dewetting growth of water ice was observed: large three-dimensional ice grains are formed. The activation energy of water desorption from the hydrogen-saturated Rh(111) surface is estimated to be 51 kJ/mol. The initial sticking probability of water decreases from 0.46 on the clean surface to 0.35 on the hydrogen-saturated surface. In IRAS measurements, D-down species were not observed on the hydrogen saturated surface. The present experimental results clearly show that a hydrophilic Rh(111) clean surface changes into a hydrophobic surface as a result of hydrogen adsorption.     

The role of molecular rotation in activated dissociative adsorption on metal surfaces

The role of molecular rotation in dissociative adsorption of H2 on the activated NiAl(110) metal surface is systematically investigated by means of classical dynamics calculations performed on ab initio six-dimensional potential energy surfaces. The calculations show that molecules rotate abruptly when they are close to the surface and that this rotation allows the molecules to adopt the orientation that is more convenient for dissociation (i.e., nearly parallel to the surface). Also, in reactive sectors of the NiAl(110) unit cell, there is an "angular threshold" below which molecules cannot dissociate. This angular threshold goes down as the incidence energy increases, which explains the rise of the dissociation probability and the fact that it reaches a value close to 1 at incidence energies of the order of 2 eV. The fact that switching on molecular rotation favors dissociation establishes a competition between dissociation and rotational excitation of reflected molecules above the dissociation threshold. Measurements on rotational excitation might thus bring indirect evidence on the dissociation dynamics. Sample calculations for nonactivated Pd(111) and activated Cu(110) metal surfaces suggest that some of these conclusions may be of general validity.     

Room-temperature kinetics of short-chain alkanethiol film growth on Ag(111) from the vapor phase

We have used time-of-flight (TOF) direct recoiling spectroscopy (DRS) to follow propanethiol adsorption at 300 K from the vapor phase on an Ag(111) surface, for exposures ranging from 10(-1) to 10(5) L. Results show that the adsorption proceeds with changes in the sticking coefficient, consistent with at least three phases. At low exposures, the alkanethiol molecules adsorb with high probability at defect sites, followed by a slower growth mode that essentially covers the whole surface. A third change in the sticking coefficient is associated with the final saturation stage, corresponding to a thicker layer related to molecules in a more upright orientation. The adsorption kinetics for hexanethiol is similar to that of propanethiol but taking place at higher rates, stressing the importance of the hydrocarbon chain length in the growth process. ISS-TOF measurements during thermal desorption show that most of the C, H, and S go away together, suggesting that the molecules adsorb and desorb from flat regions without S-C bond cleavage. Fitting the desorption maximum at 450 K with a first-order desorption curve gives a desorption energy of 1.43 eV. A small final S content that is correlated with the initial Ag(111) surface roughness is observed after desorption.     

Adsorption dynamics of water on Pt{110}-(1 x 2)

The dynamics of H(2)O adsorption on Pt{110}-(1 x 2) is studied using supersonic molecular beam and temperature programed desorption techniques. The sticking probabilities are measured using the King and Wells method at a surface temperature of 165 K. The absolute initial sticking probability s(0) of H(2)O is 0.54+/-0.03 for an incident kinetic energy of 27 kJmol. However, an unusual molecular beam flux dependence on s(0) is also found. At low water coverage (theta<1), the sticking probability is independent of coverage due either to diffusion in an extrinsic precursor state formed above bilayer islands or to incorporation into the islands. We define theta=1 as the water coverage when the dissociative sticking probability of D(2) on a surface predosed with water has dropped to zero. The slow falling H(2)O sticking probability at theta>1 results from compression of the bilayer and the formation of multilayers. Temperature programed desorption of water shows fractional order kinetics consistent with hydrogen-bonded islands on the surface. A remarkable dependence of the initial sticking probability on the translational (1-27 kJ/mol) and internal energies of water is observed: s(0) is found to be essentially a step function of translational energy, increasing fivefold at a threshold energy of 5 kJ/mol. The threshold migrates to higher energies with increasing nozzle temperature (300-700 K). We conclude that both rotational state and rotational alignment of the water molecules in the seeded supersonic expansion are implicated in dictating the adsorption process.     

Molecular-dynamics study of photodissociation of water in crystalline and amorphous ices

We present the results of classical dynamics calculations performed to study the photodissociation of water in crystalline and amorphous ice surfaces at a surface temperature of 10 K. A modified form of a recently developed potential model for the photodissociation of a water molecule in ice [S. Andersson et al., Chem. Phys. Lett. 408, 415 (2005)] is used. Dissociation in the top six monolayers is considered. Desorption of H(2)O has a low probability (less than 0.5% yield per absorbed photon) for both types of ice. The final outcome strongly depends on the original position of the photodissociated molecule. For molecules in the first bilayer of crystalline ice and the corresponding layers in amorphous ice, desorption of H atoms dominates. In the second bilayer H atom desorption, trapping of the H and OH fragments in the ice, and recombination of H and OH are of roughly equal importance. Deeper into the ice H atom desorption becomes less important and trapping and recombination dominate. Motion of the photofragments is somewhat more restricted in amorphous ice. The distribution of distances traveled by H atoms in the ice peaks at 6-7 Angstroms with a tail going to about 60 Angstroms for both types of ice. The mobility of OH radicals is low within the ice with most probable distances traveled of 2 and 1 Angstrom for crystalline and amorphous ices, respectively. OH is, however, quite mobile on top of the surface, where it has been found to travel more than 80 Angstroms. Simulated absorption spectra of crystalline ice, amorphous ice, and liquid water are found to be in very good agreement with the experiments. The outcomes of photodissociation in crystalline and amorphous ices are overall similar, but with some intriguing differences in detail. The probability of H atoms desorbing is 40% higher from amorphous than from crystalline ice and the kinetic-energy distribution of the H atoms is on average 30% hotter for amorphous ice. In contrast, the probability of desorption of OH radicals from crystalline ice is much higher than that from amorphous ice.     

Interaction between water molecules and zinc sulfide nanoparticles studied by temperature-programmed desorption and molecular dynamics simulations

We have investigated the bonding of water molecules to the surfaces of ZnS nanoparticles (approximately 2-3 nm sphalerite) using temperature-programmed desorption (TPD). The activation energy for water desorption was derived as a function of the surface coverage through kinetic modeling of the experimental TPD curves. The binding energy of water equals the activation energy of desorption if it is assumed that the activation energy for adsorption is nearly zero. Molecular dynamics (MD) simulations of water adsorption on 3 and 5 nm sphalerite nanoparticles provided insights into the adsorption process and water binding at the atomic level. Water binds with the ZnS nanoparticle surface mainly via formation of Zn-O bonds. As compared with bulk ZnS crystals, ZnS nanoparticles can adsorb more water molecules per unit surface area due to the greatly increased curvature, which increases the distance between adjacent adsorbed molecules. Results from both TPD and MD show that the water binding energy increases with decreasing the water surface coverage. We attribute the increase in binding energy with decreasing surface water coverage to the increasing degree of surface under-coordination as removal of water molecules proceeds. MD also suggests that the water binding energy increases with decreasing particle size due to the further distance and hence lower interaction between adsorbed water molecules on highly curved smaller particle surfaces. Results also show that the binding energy, and thus the strength of interaction of water, is highest in isolated nanoparticles, lower in nanoparticle aggregates, and lowest in bulk crystals. Given that water binding is driven by surface energy reduction, we attribute the decreased binding energy for aggregated as compared to isolated particles to the decrease in surface energy that occurs as the result of inter-particle interactions.     

Theoretical calculations of CH4 and H2 associative desorption from Ni(111): could subsurface hydrogen play an important role?

The results of theoretical calculations of associative desorption of CH(4) and H(2) from the Ni(111) surface are presented. Both minimum-energy paths and classical dynamics trajectories were generated using density-functional theory to estimate the energy and atomic forces. In particular, the recombination of a subsurface H atom with adsorbed CH(3) (methyl) or H at the surface was studied. The calculations do not show any evidence for enhanced CH(4) formation as the H atom emerges from the subsurface site. In fact, there is no minimum-energy path for such a concerted process on the energy surface. Dynamical trajectories started at the transition state for the H-atom hop from subsurface to surface site also did not lead to direct formation of a methane molecule but rather led to the formation of a thermally excited H atom and CH(3) group bound to the surface. The formation (as well as rupture) of the H-H and C-H bonds only occurs on the exposed side of a surface Ni atom. The transition states are quite similar for the two molecules, except that in the case of the C-H bond, the underlying Ni atom rises out of the surface plane by 0.25 A. Classical dynamics trajectories started at the transition state for desorption of CH(4) show that 15% of the barrier energy, 0.8 eV, is taken up by Ni atom vibrations, while about 60% goes into translation and 20% into vibration of a desorbing CH(4) molecule. The most important vibrational modes, accounting for 90% of the vibrational energy, are the four high-frequency CH(4) stretches. By time reversibility of the classical trajectories, this means that translational energy is most effective for dissociative adsorption at low-energy characteristic of thermal excitations but energy in stretching modes is also important. Quantum-mechanical tunneling in CH(4) dissociative adsorption and associative desorption is estimated to be important below 200 K and is, therefore, not expected to play an important role under typical conditions. An unexpected mechanism for the rotation of the adsorbed methyl group was discovered and illustrated a strong three-center C-H-Ni contribution to the methyl-surface bonding.     

Off-normal incidence dissociative sticking of H2 on Cu(100) studied using six-dimensional quantum calculations

Six-dimensional quantum calculations of the sticking probability for H2 hitting a Cu(100) surface with off-normal incidence are presented. The multiconfiguration time-dependent Hartree approach is employed for an efficient wave-packet propagation. The sticking probability is calculated for different initial momenta parallel to the surface. In contrast with the picture described in the literature, the sticking probability was found to depend on the parallel momentum. The results are explained by the topology of the potential-energy surface, which shows significant corrugation with a moderate variation of the barrier height with the surface site.     

Molecular dynamics simulations of D2O ice photodesorption

Molecular dynamics (MD) calculations have been performed to study the ultraviolet (UV) photodissociation of D(2)O in an amorphous D(2)O ice surface at 10, 20, 60, and 90 K, in order to investigate the influence of isotope effects on the photodesorption processes. As for H(2)O, the main processes after UV photodissociation are trapping and desorption of either fragments or D(2)O molecules. Trapping mainly takes place in the deeper monolayers of the ice, whereas desorption occurs in the uppermost layers. There are three desorption processes: D atom, OD radical, and D(2)O molecule photodesorption. D(2)O desorption takes places either by direct desorption of a recombined D(2)O molecule, or when an energetic D atom produced by photodissociation kicks a surrounding D(2)O molecule out of the surface by transferring part of its momentum. Desorption probabilities are calculated for photoexcitation of D(2)O in the top four monolayers and are compared quantitatively with those for H(2)O obtained from previous MD simulations of UV photodissociation of amorphous water ice at different ice temperatures [Arasa et al., J. Chem. Phys. 132, 184510 (2010)]. The main conclusions are the same, but the average D atom photodesorption probability is smaller than that of the H atom (by about a factor of 0.9) because D has lower kinetic energy than H, whereas the average OD radical photodesorption probability is larger than that of OH (by about a factor of 2.5-2.9 depending on ice temperature) because OD has higher translational energy than OH for every ice temperature studied. The average D(2)O photodesorption probability is larger than that of H(2)O (by about a factor of 1.4-2.3 depending on ice temperature), and this is entirely due to a larger contribution of the D(2)O kick-out mechanism. This is an isotope effect: the kick-out mechanism is more efficient for D(2)O ice, because the D atom formed after D(2)O photodissociation has a larger momentum than photogenerated H atoms from H(2)O, and D transfers momentum more easily to D(2)O than H to H(2)O. The total (OD + D(2)O) yield has been compared with experiments and the total (OH + H(2)O) yield from previous simulations. We find better agreement when we compare experimental yields with calculated yields for D(2)O ice than when we compare with calculated yields for H(2)O ice.     

H2 dissociation on individual Pd atoms deposited on Cu(111)

We present a Molecular Dynamics (MD) study based on Density Functional Theory (DFT) calculations for H(2) interacting with a Pd-Cu(111) surface alloy for low Pd coverages, Θ(Pd). Our results show, in line with recent experimental data, that single isolated Pd atoms evaporated on Cu(111) significantly increase the reactivity of the otherwise inert pure Cu surface. On top of substitutional Pd atoms in the Pd-Cu(111) surface alloy, the activation energy barrier for H(2) dissociation is smaller than the lowest one found on Cu(111) by a factor of two: 0.25 eV vs. 0.46 eV. Also in agreement with experiments, our DFT-MD calculations show that a large fraction of the dissociating H atoms efficiently spillover from Pd (i.e. the active sites), thanks to their extra kinetic energy due to the ~0.50 eV chemisorption exothermicity. Still, our DFT-MD calculations predict a dissociative sticking probability for low energy H(2) molecules that is much smaller than the estimated value from scanning tunneling microscopy experiments. Thus, further theoretical and experimental investigations are required for a complete understanding of H(2) dissociation on low-Θ(Pd) Pd-Cu(111) surface alloys.     

NON-REACTIVE ORIENTATIONS OF MOLECULES AT SURFACES

The role of the orientation of a molecule in its interaction with a surface is examined for the specific case of NO interaction with Pt(111). For this system molecular chemisorption occurs, mediated by a strong chemisorption well. Experimental results concerning sticking, angular distributions of scattered molecules, steric effects in scattering, and rotational excitation will be presented. Classical trajectory calculations using a model potential energy hypersurface can reproduce most experimental findings. Analysis of the trajectories shows that there is a strong orientation dependence of rotational excitation and sticking. The O-end of the molecule turns out to be non-reactive. The N-end of the molecule is very reactive. Its behaviour can almost be described using statistical methods.     

Dynamics of analyte binding onto a metallophthalocyanine: NO/FePc

The gas-surface reaction dynamics of NO impinging on an iron(II) phthalocyanine (FePc) monolayer were investigated using King and Wells sticking measurements. The initial sticking probability was measured as a function of both incident molecular beam energy (0.09-0.4 eV) and surface temperature (100-300 K). NO adsorption onto FePc saturates at 3% of a monolayer for all incident beam energies and surface temperatures, suggesting that the final chemisorption site is confined to the Fe metal centers. At low surface temperature and low incident beam energy, the initial sticking probability is 40% and decreases linearly with increasing beam energy and surface temperature. The results are consistent with the NO molecule sticking onto the FePc molecules via physisorption to the aromatics followed by diffusion to the Fe metal center, or precursor-mediated chemisorption. The adsorption mechanism of NO onto FePc was confirmed by control studies of NO sticking onto metal-free H2Pc, inert Au111, and reactive Al111.     

Ab initio MD simulations of a prototype of methyl chloride hydrolysis with explicit consideration of three water molecules: a comparison of MD trajectories with the IRC path

Ab initio molecular dynamics simulations at the Hartree-Fock/6-31G level of theory are performed on methyl chloride hydrolysis with explicit consideration of one solute and two solvent water molecules at a temperature of 298 K. The reaction involves the formation of a reactant complex and the energy surface to the transition state is found to be simple. Two types of trajectories toward the product are observed. In the first type, the system reaches an intermediate complex (complex-P1) region after two nearly concerted proton transfers involving the attacking water molecule and the solvent water molecules. These trajectories resemble the intrinsic reaction coordinate trajectory. The thermal motion of the atoms leads the system to another intermediate complex (complex-P2) region. A second type of trajectory is found in which the system reaches the complex-P2 region directly after the proton transfers. In both of these forward trajectories, back proton transfers lead the system to a final complex-F region which resembles protonated methanol. Received: 3 July 1998 / Accepted: 2 September 1998 / Published online: 15 February 1999     

The interaction of hyperthermal argon atoms with CO-covered Ru(0001): scattering and collision-induced desorption

Hyperthermal Ar atoms were scattered under grazing incidence (θ(i) = 60°) from a CO-saturated Ru(0001) surface held at 180 K. Collision-induced desorption involving the ejection of fast CO (～1 eV) occurs. The angularly resolved in-plane CO desorption distribution has a peak along the surface normal. However, the angular distribution varies with the fractional coverage of the surface. As the total CO coverage decreases, the instantaneous desorption maximum shifts to larger outgoing angles. The results are consistent with a CO desorption process that involves lateral interaction with neighboring molecules. Furthermore, the data indicate that the incident Ar cannot readily penetrate the saturated CO overlayer. Time-of-flight measurements of scattered Ar exhibit two components-fast and slow. The slow component is most evident when scattering from the fully covered surface. The ratio and origin of these components vary with the CO coverage.     

Mechanistic study of proton transfer and HD exchange in ice films at low temperatures (100-140 K)

We have examined the elementary molecular processes responsible for proton transfer and HD exchange in thin ice films for the temperature range of 100-140 K. The ice films are made to have a structure of a bottom D(2)O layer and an upper H(2)O layer, with excess protons generated from HCl ionization trapped at the D(2)OH(2)O interface. The transport behavior of excess protons from the interfacial layer to the ice film surface and the progress of the HD exchange reaction in water molecules are examined with the techniques of low energy sputtering and Cs(+) reactive ion scattering. Three major processes are identified: the proton hopping relay, the hop-and-turn process, and molecular diffusion. The proton hopping relay can occur even at low temperatures (<120 K), and it transports a specific portion of embedded protons to the surface. The hop-and-turn mechanism, which involves the coupling of proton hopping and molecule reorientation, increases the proton transfer rate and causes the HD exchange of water molecules. The hop-and-turn mechanism is activated at temperatures above 125 K in the surface region. Diffusional mixing of H(2)O and D(2)O molecules additionally contributes to the HD exchange reaction at temperatures above 130 K. The hop-and-turn and molecular diffusion processes are activated at higher temperatures in the deeper region of ice films. The relative speeds of these processes are in the following order: hopping relay>hop and turn>molecule diffusion.