Sticking of Hyperthermal CO to the (0001) Face of Crystalline Ice

Adsorption of CO and N2 molecules at the surface of solid water. A grand canonical Monte Carlo study

The adsorption of carbon monoxide and nitrogen molecules at the surface of four forms of solid water is investigated by means of grand canonical Monte Carlo simulations. The trapping ability of crystalline Ih and low-density amorphous ices, along with clathrate hydrates of structures I and II, is compared at temperatures relevant for astrophysics. It is shown that when considering a gas phase that contains mixtures of carbon monoxide and nitrogen, the trapping of carbon monoxide is favored with respect to nitrogen at the surface of all solids, irrespective of the temperature. The results of the calculations also indicate that some amounts of molecules can be incorporated in the bulk of the water structures, and the molecular selectivity of the incorporation process is investigated. Again, it is shown that incorporation of carbon monoxide is favored with respect to nitrogen in most of the situations considered here. In addition, the conclusions of the present simulations emphasize the importance of the strength of the interactions between the guest molecules and the water network. They indicate that the accuracy of the corresponding interaction potentials is a key point, especially for simulating clathrate selectivity. This highlights the necessity of having interaction potential models that are transferable to different water environments.

Adsorption of H2 on amorphous solid water studied with molecular dynamics simulations

We investigated the behavior of H2, the main constituent of the gas phase in dense clouds, after collision with amorphous solid water (ASW) surfaces, one of the most abundant chemical species of interstellar ices. We developed a general framework to study the adsorption dynamics of light species on interstellar ices. We provide binding energies and their distribution, sticking probabilities for incident energies between 1 meV and 60 meV, and thermal sticking coefficients between 10 and 300 K for surface temperatures from 10 to 110 K. We found that the sticking probability depends strongly on the adsorbate kinetic energy and the surface temperature, but hardly on the angle of incidence. We observed finite sticking probabilities above the thermal desorption temperature. Adsorption and thermal desorption should be considered as separate events with separate time scales. Laboratory results for these species have shown a gap in the trends attributed to the differently utilized experimental techniques. Our results complement observations and extend them, increasing the range of gas temperatures under consideration. We plan to use our method to study a variety of adsorbates, including radicals and charged species.

Surface chemical reactions induced by well-controlled molecular beams: translational energy and molecular orientation control

I review our recent studies of chemical reactions on single-crystalline Cu and Si surfaces induced by hyperthermal oxygen molecular beams and by oriented molecular beams, respectively. Studies of oxide formation on Cu induced by hyperthermal molecular beams suggest that the translational energy of the incident molecules plays a significant role. The use of hyperthermal molecular beams enables us to open up new chemical reaction paths, and to develop new methods for the fabrication of thin films. Oriented molecular beams also demonstrate the possibility for controlling surface chemical reactions by varying the orientation of the incident molecules. The steric effects found on Si surfaces hint at new ways of achieving material fabrication on Si surfaces. Controlling the initial conditions of incoming molecules is a powerful tool for creating new materials on surfaces with well-controlled chemical reactions.

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.

Efficient penetration of the basal plane (0001) face of ice Ih by HF at Ts=150 K: Dependence on incidence energy, incidence angle, and rotational energy

Classical trajectory simulations are carried out to investigate the influence of incidence energy, incidence angle, and rotational energy on the penetration of the basal plane (0001) face of ice I(h) by HF at a surface temperature (T(s)) of 150 K. The interaction of HF with ice is modelled by pair interactions, with the pair potential fitted to ab initio (Hartree-Fock+MP2) calculations. The penetration of ice by HF occurs already at very low incidence energies, viz., E(i)>/=20 kJ mol(-1). This is much lower than the threshold incidence energy obtained for penetration of ice by HCl (E(i) approximately 96.5 kJ mol(-1)); the calculated average barrier to penetration of ice by HF is 16.0 kJ mol(-1) and is much lower than that previously reported for HCl. As was the case for HCl, penetration of ice by HF decreases with decreasing incidence energy and increasing incidence angle. Though in general, the penetration probability is independent of the molecule's initial rotational energy, penetration beyond the second bilayer (deep penetration) is suppressed by initial rotation. This suggests that, like was found for HCl, the steering operative in deep penetration is inhibited by initial rotation. Finally, because HF is a weak acid experimental observation of HF penetrated into ice may well be possible using infrared spectroscopy, and we suggest experiments along this line.

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.