Competitive Adsorption of Trace Gases on Ice at Tropospheric Temperatures: A Grand Canonical Monte Carlo Simulation Study.

The coadsorption of two atmospheric trace gases on ice is characterized by using, for the first time, grand canonical Monte Carlo (GCMC) simulations performed in conditions similar to those of the corresponding experiments. Adsorption isotherms are simulated at tropospheric temperatures by considering two different gas mixtures of 1-butanol and acetic acid molecules, and selectivity of the ice surface with respect to these species is interpreted at the molecular scale as resulting from a competition process between these molecules for being adsorbed at the ice surface. It is thus shown that the trapping of acetic acid molecules on ice is always favored with respect to that of 1-butanol at low pressures, corresponding to low coverage of the surface, whereas the adsorption of the acid species is significantly modified by the presence of the alcohol molecules in the saturated portion of the adsorption isotherm, in accordance with the experimental observations. The present GCMC simulations thus confirm that competitive adsorption effects have to be taken into consideration in real situations when gas mixtures present in the troposphere interact with the surface of ice particles.

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 two and five 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 analyzed at 228 K show that, at low coverage of the ice surface, the polar head of the adsorbed molecules tends to optimize its hydrogen bonding with the surrounding water, whereas the aliphatic chain lies 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 that 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 under tropospheric conditions, opens the way for its use in modeling studies also relevant to another (e.g., astrophysical) context.

DFT Study of the Formation of Atmospheric Aerosol Precursors from the Interaction between Sulfuric Acid and Benzenedicarboxylic Acid Molecules.

Dicarboxylic acids are ubiquitous products of the photooxidation of volatile organic compounds which are believed to play a significant role in the formation of secondary organic aerosols in the atmosphere. In this paper, we report high-level quantum investigations of the clustering properties of sulfuric acid and benzenedicarboxylic acid molecules. Up to four molecules have been considered in the calculations, and the behavior of the three isomers of the organic diacid species have been compared. The most stable geometries have been characterized together with the corresponding thermodynamic data. From an atmospheric point of view, the results of the DFT calculations show that the organic diacid molecules may significantly enhance the nucleation of small atmospheric clusters, at least from an energetic point of view. In this respect, the phthalic acid isomer seems more efficient than the two other isomers of the benzenedicarboxylic acid, in particular because the internal distance between the two carboxyl groups in the organic diacids appears to play an important role in the stabilization of the H-bond network inside the corresponding heterocluster formed with sulfuric acid molecules.

Influence of Onion-like Carbonaceous Particles on the Aggregation Process of Hydrocarbons.

Molecular dynamics simulations are performed to characterize the nucleation behavior of organic compounds in the gas phase. Six basic molecular species are considered-ethylene, propylene, toluene, styrene, ethylbenzene, and para-xylene-in interaction with onion-like carbon nanostructures that model soot nanoparticles (NPs) at room temperature. We identify a shell-to-island aggregation process during the physisorption of aromatic molecules on the soot surface: The molecules tend to first cover the NP in a shell, on top of which additional adsorbates form island-shaped aggregates. We present results for the binding energy, suggesting that the NPs lead to the formation of more stable molecular aggregates in comparison with the pure gas phase. Our findings describe a plausible microscopic mechanism for the active role of soot in the formation and growth of organic particulate matter.

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 (Ih) 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 Ih ice under identical thermodynamic conditions, due to the larger surface area exposed to the adsorption. Indeed, the monomolecular adsorption capacity of the LDA and Ih ice surfaces is found to be 10.5 ± 0.7 µmol/m2 and 9.4 µmol/m2, 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 Ih ice at 200 K.

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

The adsorption of methylamine at the surface of amorphous ice is studied at various temperatures, ranging from 20 to 200 K, by grand canonical Monte Carlo simulations under conditions that are characteristic to the interstellar medium (ISM). The results are also compared with those obtained earlier on crystalline ( Ih) ice. We found that methylamine has a strong ability of being adsorbed on amorphous ice, involving also multilayer adsorption. The decrease of the temperature leads to a substantial increase of this adsorption ability; thus, considerable adsorption is seen at 20-50 K even at bulk gas phase concentrations that are comparable with that of the ISM. Further, methylamine molecules can also be dissolved in the bulk amorphous ice phase. Both the adsorption capacity of amorphous ice and the strength of the adsorption on it are found to be clearly larger than those corresponding to crystalline ( Ih) ice, due to the molecular scale roughness of the amorphous ice surface as well as to the lack of clear orientational preferences of the water molecules at this surface. Thus, the surface density of the saturated adsorption monolayer is estimated to be 12.6 ± 0.4 µmol/m2, 20% larger than the value of 10.35 µmol/m2, obtained earlier for Ih ice, and at low enough surface coverages the adsorbed methylamine molecules are found to easily form up to three hydrogen bonds with the surface water molecules. The estimated heat of adsorption at infinitely low surface coverage is calculated to be -69 ± 5 kJ/mol, being rather close to the estimated heat of solvation in the bulk amorphous ice phase of -74 ± 7 kJ/mol, indicating that there are at least a few positions at the surface where the adsorbed methylamine molecules experience a bulk-like local environment.

Molecular Dynamics Simulations of the Interaction between Water Molecules and Aggregates of Acetic or Propionic Acid Molecules.

Water adsorption around small acetic and propionic acid aggregates has been studied by means of molecular dynamics simulation in the temperature range of 100-265 K as a function of the water content. Calculations have shown that acetic and propionic acid molecules behave similarly and that both the temperature and the water content have a strong influence on the behavior of the corresponding systems. Two situations have been evidenced for the acid-water aggregates, corresponding either to water adsorption on large acid grains at very low temperatures or to the formation of droplets consisting of acid molecules adsorbed at the surface of water aggregates at higher temperatures and high water content. At low water content and high temperature, only a partial mixing between water and acid molecules has been observed. The results of the present simulations emphasize the need for further experimental and simulation works to achieve a better characterization of the effects of both temperature and humidity on the behavior of organic aerosols in the troposphere.

The effect of anaesthetics on the properties of a lipid membrane in the biologically relevant phase: a computer simulation study.

Molecular dynamics simulations of the fully hydrated neat dipalmitoylphosphatidylcholine (DPPC) membrane as well as DPPC membranes containing four different general anaesthetic molecules, namely chloroform, halothane, diethyl ether and enflurane, have been simulated at two different pressures, i.e., at 1 bar and 1000 bar, at the temperature of 310 K. At this temperature the model used in this study is known to be in the biologically most relevant liquid crystalline (Lα) phase. To find out which properties of the membrane might possibly be related to the molecular mechanism of anaesthesia, we have been looking for properties that change in the same way in the presence of any general anaesthetic molecule, and change in the opposite way by the increase of pressure. This way, we have ruled out the density distribution of various groups along the membrane normal axis, orientation of the lipid heads and tails, self-association of the anaesthetics, as well as the local order of the lipid tails as possible molecular reasons of anaesthesia. On the other hand, we have found that the molecular surface area, and hence also the molecular volume of the membrane, is increased by the presence of any anaesthetic molecule, and decreased by the pressure, in accordance with the more than half a century old critical volume hypothesis. We have also found that anaesthetic molecules prefer two different positions along the membrane normal axis, namely the middle of the membrane and the outer edge of the hydrocarbon region, close to the polar headgroups. The increase of pressure is found to decrease the former, and increase the latter preference, and hence it might also be related to the pressure reversal of anaesthesia.

Methane clathrates in the solar system.

We review the reservoirs of methane clathrates that may exist in the different bodies of the Solar System. Methane was formed in the interstellar medium prior to having been embedded in the protosolar nebula gas phase. This molecule was subsequently trapped in clathrates that formed from crystalline water ice during the cooling of the disk and incorporated in this form into the building blocks of comets, icy bodies, and giant planets. Methane clathrates may play an important role in the evolution of planetary atmospheres. On Earth, the production of methane in clathrates is essentially biological, and these compounds are mostly found in permafrost regions or in the sediments of continental shelves. On Mars, methane would more likely derive from hydrothermal reactions with olivine-rich material. If they do exist, martian methane clathrates would be stable only at depth in the cryosphere and sporadically release some methane into the atmosphere via mechanisms that remain to be determined. In the case of Titan, most of its methane probably originates from the protosolar nebula, where it would have been trapped in the clathrates agglomerated by the satellite's building blocks. Methane clathrates are still believed to play an important role in the present state of Titan. Their presence is invoked in the satellite's subsurface as a means of replenishing its atmosphere with methane via outgassing episodes. The internal oceans of Enceladus and Europa also provide appropriate thermodynamic conditions that allow formation of methane clathrates. In turn, these clathrates might influence the composition of these liquid reservoirs. Finally, comets and Kuiper Belt Objects might have formed from the agglomeration of clathrates and pure ices in the nebula. The methane observed in comets would then result from the destabilization of clathrate layers in the nuclei concurrent with their approach to perihelion. Thermodynamic equilibrium calculations show that methane-rich clathrate layers may exist on Pluto as well. Key Words: Methane clathrate-Protosolar nebula-Terrestrial planets-Outer Solar System. Astrobiology 15, 308-326.

Water and formic acid aggregates: a molecular dynamics study.

Water adsorption around a formic acid aggregate has been studied by means of molecular dynamics simulations in a large temperature range including tropospheric conditions. Systems of different water contents have been considered and a large number of simulations has allowed us to determine the behavior of the corresponding binary formic acid-water systems as a function of temperature and humidity. The results clearly evidence a threshold temperature below which the system consists of water molecules adsorbed on a large formic acid grain. Above this temperature, formation of liquid-like mixed aggregates is obtained. This threshold temperature depends on the water content and may influence the ability of formic acid grains to act as cloud condensation nuclei in the Troposphere.

On the abundances of noble and biologically relevant gases in Lake Vostok, Antarctica.

Motivated by the possibility of comparing theoretical predictions of Lake Vostok's composition with future in situ measurements, we investigated the composition of clathrates that are expected to form in this environment from the air supplied to the lake by melting ice. To establish the best possible correlation between the lake water composition and that of air clathrates formed in situ, we used a statistical thermodynamic model based on the description of the guest-clathrate interaction by a spherically averaged Kihara potential with a nominal set of potential parameters. We determined the fugacities of the different volatiles present in the lake by defining a "pseudo" pure substance dissolved in water owning the average properties of the mixture and by using the Redlich-Kwong equation of state to mimic its thermodynamic behavior. Irrespective of the clathrate structure considered in our model, we found that xenon and krypton are strongly impoverished in the lake water (a ratio in the 0.04-0.1 range for xenon and a ratio in the ≈ 0.15-0.3 range for krypton) compared to their atmospheric abundances. Argon and methane were also found to be depleted in the Lake Vostok water by factors in the 0.5-0.95 and 0.3-0.5 ranges, respectively, compared to their atmospheric abundances. On the other hand, the carbon dioxide abundance was found to be substantially enriched in the lake water compared to its atmospheric abundance (by a factor in the 1.6-5 range at 200 residence times). The comparison of our predictions of the CO2 and CH4 mole fractions in Lake Vostok with future in situ measurements will allow disentangling between the possible supply sources.

Molecular dynamics simulations of the water adsorption around malonic acid aerosol models.

Water nucleation around a malonic acid aggregate has been studied by means of molecular dynamics simulations in the temperature and pressure range relevant for atmospheric conditions. Systems of different water contents have been considered and a large number of simulations have allowed us to determine the phase diagram of the corresponding binary malonic acid-water systems. Two phases have been evidenced in the phase diagrams corresponding either to water adsorption on a large malonic acid grain at low temperatures, or to the formation of a liquid-like mixed aggregate of the two types of molecules, at higher temperatures. Finally, the comparison between the phase diagrams simulated for malonic acid-water and oxalic acid-water mixtures emphasizes the influence of the O : C ratio on the hydrophilic behavior of the aerosol, and thus on its ability to act as a cloud condensation nucleus, in accordance with recent experimental conclusions.

Anesthetic molecules embedded in a lipid membrane: a computer simulation study.

The effect of four general anesthetic molecules, i.e., chloroform, halothane, diethyl ether and enflurane, on the properties of a fully hydrated dipalmitoylphosphatidylcholine (DPPC) membrane is studied in detail by long molecular dynamics simulations. Furthermore, to address the problem of pressure reversal, the effect of pressure on the anesthetic containing membranes is also investigated. In order to ensure sufficient equilibration and adequate sampling, the simulations performed have been at least an order of magnitude longer than the studies reported previously in the literature on general anesthetics. The results obtained can help in resolving several long-standing contradictions concerning the effect of anesthetics, some of which were the consequence of too short simulation time used in several previous studies. More importantly, a number of seeming contradictions are found to originate from the fact that different anesthetic molecules affect the membrane structure differently in several respects. In particular, halothane, being able to weakly hydrogen bound to the ester group of the lipid tails, is found to behave in a markedly different way than the other three molecules considered. Besides, we also found that two changes, namely lateral expansion of the membrane and increasing local disorder in the lipid tails next to the anesthetic molecules, are clearly induced by all four anesthetic molecules tested here in the same way, and both of these effects are reverted by the increase in pressure.

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.

Water adsorption around oxalic acid aggregates: a molecular dynamics simulation of water nucleation on organic aerosols.

The phase behaviour of binary oxalic acid-water mixtures has been investigated by means of computer simulation techniques. Such mixtures play an important role in atmospheric processes, since the hydrogen bonding ability of oxalic acid molecules allows them to form aerosol particles. Water can in turn be readily adsorbed on the surface of such aerosol particles, which results in the formation of small ice grains. These grains are thus considered to be acting as cloud condensation nuclei, giving rise to the formation of ice clouds.

The ice-vapor interface and the melting point of ice I(h) for the polarizable POL3 water model.

We use molecular dynamics simulations to determine the melting point of ice I(h) for the polarizable POL3 water force field (Dang, L. X. J. Chem. Phys.1992, 97, 2659). Simulations are performed on a slab of ice I(h) with two free surfaces at several different temperatures. The analysis of the time evolution of the total energy in the course of the simulations at the set of temperatures yields the melting point of the POL3 model to be T(m) = 180 ± 10 K. Moreover, the results of the simulations show that the degree of hydrogen-bond disorder occurring in the bulk of POL3 ice is larger (at the corresponding degree of undercooling) than in ice modeled by nonpolarizable water models. These results demonstrate that the POL3 water force field is rather a poor model for studying ice and ice-liquid or ice-vapor interfaces. While a number of polarizable water models have been developed over the past years, little is known about their performance in simulations of supercooled water and ice. This study thus highlights the need for testing of the existing polarizable water models over a broad range of temperatures, pressures, and phases, and developing a new polarizable water force field, reliable over larger areas of the phase diagram.

Adsorption of hydroxyacetone on pure ice surfaces.

The adsorption of hydroxyacetone molecules at the surface of ice is investigated by means of flow-tube reactor measurements in the temperature range: 213-253 K. The number of molecules adsorbed per surface unit is conventionally plotted as a function of the absolute gas concentration of hydroxyacetone and is compared to that previously obtained for acetone and ethanol. The enthalpy of adsorption and the monolayer capacity at the ice surface are determined. In addition, molecular dynamics simulations are performed to support the experimental results. However, it is shown that the available interaction potential between hydroxyacetone and ice is not accurate enough to allow a robust detailed analysis of the adsorption process. Finally, a rapid estimation of the hydroxyacetone partitioning between the gas phase and ice shows that in the densest ice clouds, up to 29% of hydroxyacetone could be adsorbed on pure ice surfaces at 203 K.

Water adsorption on oxidized single atomic vacancies present at the surface of small carbonaceous nanoparticles modeling soot.

Quantum calculations are used to study the interaction of water molecules with carbonaceous clusters containing one single carbon atom vacancy. This is a simple but realistic way to model the active surfaces found in soot emitted by aircrafts. Prior to water adsorption, the atomic vacancy is oxidised by an approaching oxygen molecule, which is also likely to occur behind planes. The results of the calculations show that this oxidation process results in the formation of one ketone-like site and one epoxide-like site around the atomic vacancy. These sites may act as nucleation centers for water molecules, which are, however, physisorbed on the oxidized surface, leading to very weak charge transfer with the surface. Although less attractive for water than, for instance, a carboxyl-like site, the ketone-like site can also participate in the hydrophilic behavior of soot primary particles. In contrast, the epoxide-like site formed around the vacancy shows a very low affinity for water molecules.

Water adsorption isotherms on porous onionlike carbonaceous particles. Simulations with the grand canonical Monte Carlo method.

The grand canonical Monte Carlo method is used to simulate the adsorption isotherms of water molecules on different types of model soot particles. These soot models are constructed by first removing atoms from onion-fullerene structures in order to create randomly distributed pores inside the soot, and then performing molecular dynamics simulations, based on the reactive adaptive intermolecular reactive empirical bond order (AIREBO) description of the interaction between carbon atoms, to optimize the resulting structures. The obtained results clearly show that the main driving force of water adsorption on soot is the possibility of the formation of new water-water hydrogen bonds with the already adsorbed water molecules. The shape of the calculated water adsorption isotherms at 298 K strongly depends on the possible confinement of the water molecules in pores of the carbonaceous structure. We found that there are two important factors influencing the adsorption ability of soot. The first of these factors, dominating at low pressures, is the ability of the soot of accommodating the first adsorbed water molecules at strongly hydrophilic sites. The second factor concerns the size and shape of the pores, which should be such that the hydrogen bonding network of the water molecules filling them should be optimal. This second factor determines the adsorption properties at higher pressures.