Adsorption of Methylamine on Amorphous Ice under Interstellar Conditions. A Grand Canonical Monte Carlo Simulation Study

Adsorption of C2-C5 alcohols on ice. A grand canonical Monte Carlo simulation study

In this paper, we report Grand Canonical Monte Carlo simulations performed to characterize the adsorption of four linear alcohol molecules, comprising between 2 and 5 carbon atoms (namely, ethanol, n-propanol, n-butanol, and n-pentanol) on crystalline ice in a temperature range typical of the Earth's troposphere.The adsorption details analysed at 228 K show that, at low coverage of the ice surface, the polar head of the adsorbed molecules tend to optimize its hydrogen bonding with the surrounding water, whereas the aliphatic chain lie more or less parallel to the ice surface. With increasing coverage, the lateral interactions between the adsorbed alcohol molecules lead to the reorientation of the aliphatic chains which tend to become perpendicular to the surface, the adsorbed molecules pointing thus their terminal methyl group up to the gas phase. When compared to the experimental data, the simulated and measured isotherms show a very good agreement, although a small temperature shift between simulations and experiments could be inferred from simulations at various temperatures. In addition, this agreement appears to be better for ethanol and n-propanol than for n-butanol and n-pentanol, especially at the highest pressures investigated, pointing to a possible slight underestimation of the lateral interactions between the largest alcohol molecules by the interaction potential model used. Nevertheless, the global accuracy of the approach used, as tested in tropospheric conditions, opens the way for its use in modeling studies also relevant to another (e.g., astrophysical) context.

Computer simulation investigation of the adsorption of acetamide on low density amorphous ice. An astrochemical perspective

The adsorption of acetamide on low density amorphous (LDA) ice is investigated by grand canonical Monte Carlo computer simulations at the temperatures 50, 100, and 200 K, characteristic of certain domains of the interstellar medium (ISM). We found that the relative importance of the acetamide-acetamide H-bonds with respect to the acetamide-water ones increases with decreasing temperature. Thus, with decreasing temperature, the existence of the stable monolayer, characterizing the adsorption at 200 K, is gradually replaced by the occurrence of marked multilayer adsorption, preceding even the saturation of the first layer at 50 K. While isolated acetamide molecules prefer to lay parallel to the ice surface to maximize their H-bonding with the surface water molecules, this orientational preference undergoes a marked change upon saturation of the first layer due to increasing competition of the adsorbed molecules for H-bonds with water and to the possibility of their H-bond formation with each other. As a result, molecules stay preferentially perpendicular to the ice surface in the saturated monolayer. The chemical potential value corresponding to the point of condensation is found to decrease linearly with increasing temperature. We provide, in analogy with the Clausius-Clapeyron equation, a thermodynamic explanation of this behavior and estimate the molar entropy of condensed phase acetamide to be 34.0 J/mol K. For the surface concentration of the saturated monolayer, we obtain the value 9.1 ± 0.8 µmol/m², while the heat of adsorption at infinitely low surface coverage is estimated to be -67.8 ± 3.0 kJ/mol. Our results indicate that the interstellar formation of peptide chains through acetamide molecules, occurring at the surface of LDA ice, might well be a plausible process in the cold (i.e., below 50 K) domains of the ISM; however, it is a rather unlikely scenario in its higher temperature (i.e., 100-200 K) domains.

Description of the Interfacial Behavior of Benzonitrile at Icy Surfaces by Grand Canonical Monte Carlo Simulations

The adsorption of benzonitrile at the surface of crystalline (I_h) and low-density amorphous (LDA) ice has been investigated by grand canonical Monte Carlo simulations at temperatures ranging from 50 to 200 K. It is found that, in spite of its rather large dipole moment of 4.5 D, benzonitrile molecules can only form a highly unsaturated monolayer on LDA ice, reaching not more than 50% of the surface concentration of the saturated monolayer even at the lowest temperature considered, and they practically do not adsorb on I_h ice. In spite of the observed weak ability of the benzonitrile molecules for being adsorbed, the estimated heat of adsorption at an infinitely low surface concentration of -66.8 ± 2.2 kJ/mol is rather large. This value includes the contribution of roughly -30 to -35 kJ/mol of a benzene ring, about -10 kJ/mol of a large molecular dipole moment, and about -20 to -25 kJ/mol of a benzonitrile-water H-bond, as estimated from comparisons with the heat of adsorption values of similar molecules. The surprisingly weak ability of benzonitrile for adsorption is thus attributed to the unusually strong cohesion between the molecules, considerably exceeding their adhesion to ice, as reflected in the 70-80 kJ/mol difference of the lateral and ice contributions to the binding energy of surface benzonitrile molecules in the presence of condensed benzonitrile.

Single-Particle Dynamics at the Intrinsic Surface of Aqueous Alkali Halide Solutions

The distribution of ions in the proximity of the liquid-vapor interface of their aqueous solution has been the subject of an intense debate during the last decade. The effects of ionic polarizability have been one of its salient aspects. Much less has been said about the corresponding dynamical properties, which are substantially unexplored. Here, we investigate the single-particle dynamics at the liquid-vapor interface of several alkali halide solutions, using molecular dynamics simulations with polarizable and nonpolarizable force fields and intrinsic surface analysis. We analyze the diffusion coefficient, residence time, and velocity autocorrelation function of water and ions and investigate how these properties depend on the molecular layer where they reside. While anions are found in the first molecular layer for relatively long times, cations are only making quick excursions into it, thanks to thermal fluctuations. The in-layer residence time of ions and their molar fraction in the layer turned out to be linearly dependent on each other. We interpret this unexpected result using a simple two-state model. In addition, we found that, unlike water and other neat molecular liquids that show a different diffusion mechanism at the surface than in the bulk of their liquid phase, ions do not enjoy enhanced mobility in the surface layer of their aqueous solution. This result indicates that ions in the surface layer are shielded by their nearest water neighbors from being exposed to the vapor phase as much as possible. Such positions are available for the ions at the negatively curved troughs of the molecularly rugged liquid surface.

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 Formamide at the Surface of Amorphous and Crystalline Ices under Interstellar and Tropospheric Conditions. A Grand Canonical Monte Carlo Simulation Study

The adsorption of formamide is studied both at the surface of crystalline (I_h) ice at 200 K and at the surface of low density amorphous (LDA) ice in the temperature range of 50-200 K by grand canonical Monte Carlo (GCMC) simulation. These systems are characteristic of the upper troposphere and of the interstellar medium (ISM), respectively. Our results reveal that while no considerable amount of formamide is dissolved in the bulk ice phase in any case, the adsorption of formamide at the ice surface under these conditions is a very strongly preferred process, which has to be taken into account when studying the chemical reactivity in these environments. The adsorption is found to lead to the formation of multimolecular adsorption layer, the occurrence of which somewhat precedes the saturation of the first molecular layer. Due to the strong lateral interaction acting between the adsorbed formamide molecules, the adsorption isotherm does not follow the Langmuir shape. Adsorption is found to be slightly stronger on LDA than I_h ice under identical thermodynamic conditions, due to the larger surface area exposed to the adsorption. Indeed, the monomolecular adsorption capacity of the LDA and I_h ice surfaces is found to be 10.5 ± 0.7 μmol/m² and 9.4 μmol/m², respectively. The first layer formamide molecules are very strongly bound to the ice surface, forming typically four hydrogen bonds with each other and the surface water molecules. The heat of adsorption at infinitely low surface coverage is found to be -105.6 kJ/mol on I_h ice at 200 K.