Interaction of D2 with H2O amorphous ice studied by temperature-programed desorption experiments

Theoretical Distribution of the Ammonia Binding Energy at Interstellar Icy Grains: A New Computational Framework

The binding energies (BE) of molecules on the interstellar grains are crucial in the chemical evolution of the interstellar medium (ISM). Both temperature-programmed desorption (TPD) laboratory experiments and quantum chemistry computations have often provided, so far, only single values of the BE for each molecule. This is a severe limitation, as the ices enveloping the grain mantles are structurally amorphous, giving rise to a manifold of possible adsorption sites, each with different BEs. However, the amorphous ice nature prevents the knowledge of structural details, hindering the development of a common accepted atomistic icy model. In this work, we propose a computational framework that closely mimics the formation of the interstellar grain mantle through a water by water accretion. On that grain, an unbiased random (but well reproducible) positioning of the studied molecule is then carried out. Here we present the test case of NH₃, a ubiquitous species in the molecular ISM. We provide the BE distribution computed by a hierarchy approach, using the semiempirical xTB-GFN2 as a low-level method to describe the whole icy cluster in combination with the B97D3 DFT functional as a high-level method on the local zone of the NH₃ interaction. The final ZPE-corrected BE is computed at the ONIOM(DLPNO-CCSD(T)//B97D3:xTB-GFN2) level, ensuring the best cost/accuracy ratio. The main peak of the predicted NH₃ BE distribution is in agreement with experimental TPD and computed data in the literature. A second broad peak at very low BE values is also present, which has never been detected before. It may provide the solution to a longstanding puzzle about the presence of gaseous NH₃ also observed in cold ISM objects.

Surface Temperature Dependence of Hydrogen Ortho-Para Conversion on Amorphous Solid Water

The surface temperature dependence of the ortho-to-para conversion of H_{2} on amorphous solid water is first reported. A combination of photostimulated desorption and resonance-enhanced multiphoton ionization techniques allowed us to sensitively probe the conversion on the surface of amorphous solid water at temperatures of 9.2-16 K. Within a narrow temperature window of 8 K, the conversion time steeply varied from ∼4.1×10^{3} to ∼6.4×10^{2}  s. The observed temperature dependence is discussed in the context of previously suggested models and the energy dissipation process. The two-phonon process most likely dominates the conversion rate at low temperatures.

Hindered rotational physisorption states of H2 on Ag(111) surfaces

We have investigated the physisorption states of H2 on Ag(111) surfaces. To clarify the accurate adsorption properties of H2 on Ag(111), we performed first-principles calculations based on spin-polarized density functional theory (DFT) with the semiempirical DFT-D2 method and the newly-developed exchange functional with the non-local correlation functional vdW-DF2 (rev-vdW-DF2). We constructed exhaustive potential energy surfaces, and revealed that non-negligible out-of-plane potential anisotropy with a perpendicular orientation preference exists even for H2 physisorption on planar Ag(111), as predicted by previous results of resonance-enhanced multiphoton ionization spectroscopy and temperature-programmed desorption experiments. Therefore, the molecular rotational ground states of ortho-H2 split into two energy levels in the anisotropic potential. The obtained adsorption energy and the number of bound states, including the zero-point energies and the rotational energy shift, agree with diffractive and rotationally mediated selective adsorption scattering resonance measurements. The origin of the potential anisotropy on Ag(111) is a combination of the London dispersion interaction and the virtual transition of the metal electron to the unoccupied molecular state.

Physisorption and desorption of H2, HD and D2 on amorphous solid water ice. Effect on mixing isotopologue on statistical population of adsorption sites

We study the adsorption and desorption of three isotopologues of molecular hydrogen mixed on 10 ML of porous amorphous water ice (ASW) deposited at 10 K.

Application of a diffusion-desorption rate equation model in astrochemistry

Desorption and diffusion are two of the most important processes on interstellar grain surfaces; knowledge of them is critical for the understanding of chemical reaction networks in the interstellar medium (ISM). However, a lack of information on desorption and diffusion is preventing further progress in astrochemistry. To obtain desorption energy distributions of molecules from the surfaces of ISM-related materials, one usually carries out adsorption-desorption temperature programmed desorption (TPD) experiments, and uses rate equation models to extract desorption energy distributions. However, the often-used rate equation models fail to adequately take into account diffusion processes and thus are only valid in situations where adsorption is strongly localized. As adsorption-desorption experiments show that adsorbate molecules tend to occupy deep adsorption sites before occupying shallow ones, a diffusion process must be involved. Thus, it is necessary to include a diffusion term in the model that takes into account the morphology of the surface as obtained from analyses of TPD experiments. We take the experimental data of CO desorption from the MgO(100) surface and of D2 desorption from amorphous solid water ice as examples to show how a diffusion-desorption rate equation model explains the redistribution of adsorbate molecules among different adsorption sites. We extract distributions of desorption energies and diffusion energy barriers from TPD profiles. These examples are contrasted with a system where adsorption is strongly localized--HD from an amorphous silicate surface. Suggestions for experimental investigations are provided.

Solid state chemistry of nitrogen oxides – Part II: surface consumption of NO2

Efficient surface destruction mechanisms (NO₂ + H/O/N), leading to solid H₂O, NH₂OH, and N₂O, can explain the non-detection of NO₂ in space.

Solid state chemistry of nitrogen oxides – Part I: surface consumption of NO

NO₂ is efficiently formed in the solid state via NO + O/O₂/O₃ reactions.

Small-angle neutron scattering study of micropore collapse in amorphous solid water

Amorphous solid water (ASW) is shown to undergo a micropore collapse from cylindrical pores (3D) to lamellae (2D) at >120 K using small-angle neutron scattering.

Efficient diffusive mechanisms of O atoms at very low temperatures on surfaces of astrophysical interest

At the low temperatures of interstellar dust grains, it is well established that surface chemistry proceeds via diffusive mechanisms of H atoms weakly bound (physisorbed) to the surface. Until recently, however, it was unknown whether atoms heavier than hydrogen could diffuse rapidly enough on interstellar grains to react with other accreted species. In addition, models still require simple reduction as well as oxidation reactions to occur on grains to explain the abundances of various molecules. In this paper we investigate O-atom diffusion and reactivity on a variety of astrophysically relevant surfaces (water ice of three different morphologies, silicate, and graphite) in the 6.5–25 K temperature range. Experimental values were used to derive a diffusion law that emphasizes that O atoms diffuse by quantum mechanical tunnelling at temperatures as low as 6.5 K. The rates of diffusion on each surface, based on modelling results, were calculated and an empirical law is given as a function of the surface temperature. The relative diffusion rates are k_H2Oice > k_sil > k_graph ≫ k_expected. The implications of efficient O-atom diffusion over astrophysically relevant time-scales are discussed. Our findings show that O atoms can scan any available reaction partners (e.g., either another H atom, if available, or a surface radical like O or OH) at a faster rate than that of accretion. Also, as dense clouds mature, H₂ becomes far more abundant than H and the O : H ratio grows, and the reactivity of O atoms on grains is such that O becomes one of the dominant reactive partners together with H.

Quantum simulations of the hydrogen molecule on ammonia clusters

Mixed ammonia-hydrogen molecule clusters [H2-(NH3)n] have been studied with the aim of exploring the quantitative importance of the H2 quantum motion in defining their structure and energetics. Minimum energy structures have been obtained employing genetic algorithm-based optimization methods in conjunction with accurate pair potentials for NH3-NH3 and H2-NH3. These include both a full 5D potential and a spherically averaged reduced surface mimicking the presence of a para-H2. All the putative global minima for n ≥ 7 are characterized by H2 being adsorbed onto a rhomboidal ammonia tetramer motif formed by two double donor and two double acceptor ammonia molecules. In a few cases, the choice of specific rhombus seems to be directed by the vicinity of an ammonia ad-molecule. Diffusion Monte Carlo simulations on a subset of the species obtained highlighted important quantum effects in defining the H2 surface distribution, often resulting in populating rhomboidal sites different from the global minimum one, and showing a compelling correlation between local geometrical features and the relative stability of surface H2. Clathrate-like species have also been studied and suggested to be metastable over a broad range of conditions if formed.

Remarkable nanoconfinement effects on chemical equilibrium manifested in nucleotide dimerization and H-D exchange reactions

Nanoconfinement entropic effects on chemical equilibrium involving a small number of molecules, which we term NCECE, are revealed by two widely diverse types of reactions. Employing statistical-mechanical principles, we show how the NCECE effect stabilizes nucleotide dimerization observed within self-assembled molecular cages. Furthermore, the effect provides the basis for dimerization even under an aqueous environment inside the nanocage. Likewise, the NCECE effect is pertinent to a longstanding issue in astrochemistry, namely the extra deuteration commonly observed for molecules reacting on interstellar dust grain surfaces. The origin of the NCECE effect is elucidated by means of the probability distributions of the reaction extent and related variations in the reactant-product mixing entropy. Theoretical modelling beyond our previous preliminary work highlights the role of the nanospace size in addition to that of the nanosystem size, namely the limited amount of molecules in the reaction mixture. Furthermore, the NCECE effect can depend also on the reaction mechanism, and on deviations from stoichiometry. The NCECE effect, leading to enhanced, greatly variable equilibrium "constants", constitutes a unique physical-chemical phenomenon, distinguished from the usual thermodynamical properties of macroscopically large systems. Being significant particularly for weakly exothermic reactions, the effects should stabilize products in other closed nanoscale structures, and thus can have notable implications for the growing nanotechnological utilization of chemical syntheses conducted within confined nanoreactors.

Molecular hydrogen desorption from amorphous surfaces at low temperature

We studied the desorption of hydrogen molecules from amorphous silicates of composition (Fe(x)Mg(1 - x))(2)SiO(4) (0 < x < 1) using thermal programmed desorption (TPD). Selected measurements of formation of molecular hydrogen on and desorption from a single crystal olivine sample were done for comparison. The experiments were conducted in conditions as close as technically possible to those found in selected interstellar medium environments, where the formation of molecular hydrogen takes place on dust grains. From molecular desorption data, we derive the energy distribution of binding sites using a direct inversion method. The application of this type of data to the study of elementary processes of migration of atoms and molecules on and ejection from disordered surfaces at low temperature is discussed.

D(2) desorption kinetics on amorphous solid water: from compact to porous ice films

The desorption kinetics of D(2) from amorphous solid water (ASW) films have been studied by the temperature-programmed desorption (TPD) technique in the 10-30 K temperature range. Compact (and nonporous) films were grown at 120 K over a copper substrate. Ultra-thin porous films were additionally grown at 10 K over the compact base. The TPD spectra from compact and from up to 20 monolayers (ML) porous films were compared. The simulation of the TPD experimental traces provides the corresponding D(2) binding-energy distributions. As compared to the compact case, the binding-energy distribution found for the 10 ML porous film clearly extends to higher energies. To study the transition from compact to porous ice, porous films of intermediate thicknesses (<10 ML), including ultra-thin films (<1 ML), were grown over the compact substrate. The thermal D(2) desorption peak was found to shift to higher temperatures as the porous ice network was progressively formed. This behavior can be explained by the formation of more energetic binding sites related to porous films. TPD spectra were also modelled by using a combination of the two energy distributions, one associated to a bare compact ice and the other associated to a 10 ML porous ice film. This analysis reveals a very fast evolution of the binding-energy distribution towards that of porous ice. Our results show that few ML of additional porous film are sufficient to produce a sample for which the D(2) adsorption can be described by the energy distribution found for the 10 ML porous film. These experiments then provide evidence that the binding energy of D(2) on ASW ice is primarily governed by the topological and morphological disorder of the surface at molecular scale.

Measurement of the adsorption energy difference between ortho- and para-D2 on an amorphous ice surface

Molecular hydrogen interaction on water ice surfaces is a major process taking place in interstellar dense clouds. By coupling laser detection and classical thermal desorption spectroscopy, it is possible to study the effect of rotation of D(2) on adsorption on amorphous solid water ice surfaces. The desorption profiles of ortho- and para-D(2) are different. This difference is due to a shift in the adsorption energy distribution of the two lowest rotational states. Molecules in J''=1 rotational state are on average more strongly bound to the ice surface than those in J''=0 rotational state. This energy difference is estimated to be 1.4+/-0.3 meV. This value is in agreement with previous calculation and interpretation. The nonspherical wave function J'' =1 has an interaction with the asymmetric part of the adsorption potential and contributes positively in the binding energy.

The formation of vibrationally excited HD from atomic recombination on cold graphite surfaces

HD molecules formed in v"=3 and v"=4 have been detected by laser spectroscopy when a cold (15 K) graphite surface is irradiated with H and D atoms. Population of the v"=3, J"=0-6 and v"=4, J"=0-6 levels has been detected and the average rotational temperatures of the nascent HD were determined. These results are compared with previous data collected for the formation of HD in v"=1 and 2 under similar conditions. This comparison indicates that the nascent HD flux increases with increasing vibrational quantum number for v"=1-4.