Properties of Ir4 clusters in the gas phase and on oxide substrates

Results of a theoretical study on the properties of Ir4 clusters in the gas–phase and on oxide surfaces are presented. The work is based on density functional theory (DFT) within the generalized gradient approximation (GGA) and ultrasoft pseudopotentials. Properties of a small particle such as Ir4 cluster are entirely determined by its geometry. The already known result that the most stable form of Ir4 in the gas–phase is the square structure which is significantly more stable than the butterfly and tetrahedron is confirmed. This result is in contradiction with experiments which indicate that the oxide supported Ir4 adopts a tetrahedral configuration. It is shown in this thesis that the chemical environment has a strong influence on the relative stability of Ir4 clusters. On MgO(100) surface, the square isomer remains the most stable Ir4 structure, well separated in energy from the other two. Moreover, the tetrahedron is heavily distorted by the interaction with the surface oxygen. Presence of point defects (neutral and charged O vacancies) affects the energy ordering making tetrahedron and square very close in energy, but the structural distortion of the tetrahedron is even bigger and the predicted data do not correspond to experiments. On TiO2(110) the tetrahedron and square structures become degenerate and the butterfly becomes the least stable isomer. Moreover, structural distortions are very small, in agreement with experimental data. It is shown that the TiO2 surface influences the relative stability of the three isomers through a particularly strong electrostatic field. Interactions of Ir4 with H, C and O atoms as well as with CO molecules have been studied. Adsorption of a single C atom strongly influences the relative stability of the three isomers. Upon C adsorption, the butterfly becomes the most stable gas–phase isomer while on both surfaces the tetrahedron is the most probable structure. Adsorption of a single H or O atom does not produce the same effect. The interaction with CO molecules is also important given the experimental procedure used for producing supported Ir4 clusters. It is shown that on MgO(100), CO dissociation is as probable as the competing process CO desorption justifying the presence of carbon adatoms on Ir4 clusters which brings theoretical predictions in better agreement with experimental data.

Role of Adsorbed H, C, O, and CO on the Atomic Structure of Free and MgO(100)-Supported Ir-4 Clusters: An ab Initio Study

Alfonso Baldereschi, Zeljko Sljivancanin, Vladan Stevanovic

We applied density functional theory based on ultrasoft pseudopotentials to study the structural properties of Ir-4 clusters both in the gas phase and adsorbed on a MgO(100) surface. To resolve the discrepancy between experimental data which suggest a tetrahedral structure for Ir-4 and theoretical results which show a strong preference for the square isomer, we investigated the effect of several adsorbates on the equilibrium atomic structure of the clusters. Calculated binding energies of a single H, C, or O atom, as well as one CO or OH molecule, to three stable Ir-4 isomers indicate that C or CO adsorption significantly influences the relative stability of Ir-4 isomers. For MgO(100)-supported Ir-4, atomic carbon is able to change the isomer preference from the square to the tetrahedral geometry.

2010

Release kinetics of volatiles from smectite clay minerals

Pascal Clausen

The focus of this thesis was on the mechanisms of release of volatile molecules from smectite clay minerals. The weight loss was monitored by thermogravimetry, differential calorimetry, and mass spectrometry under isothermal conditions and under ramp temperature conditions. Modeling of the clay-volatile release systems was performed with the help of finite element method calculations at the macroscopic level and ab initio calculations at the molecular level. A first comparison between finite element method simulations and experimental data for the evaporation of bulk volatile liquids under a convective gas flow showed the calculations to be in good agreement with experiments. Results of the simulations were used to develop a semi-analytical model explaining the dependence of the evaporation on the total pressure, the carrier gas flow, the temperature, and material constants. In particular, its temperature dependence could be approximated to good accuracy by an Arrhenius-type equation derived from the semi-analytical model. Differential calorimetry measurements of the heats of vaporization showed equilibrium conditions at the surface of the liquids to be satisfied. The same approach was extended to the release of model volatiles (water, ethanol, ethyl acetate and toluene) from smectite clays. At high coverage, the release was found to be close to that for the bulk liquids. Its decrease with time followed the behaviour observed in the respective curves of the gas/condensed phase partition coefficients, the equilibrium desorption isotherms. Equilibrium condition at the surface of the sample was evidenced by a comparison between the measured heats of vaporization and the equilibrium desorption isotherms. The differences observed in the measured equilibrium desorption isotherms of the volatile on the smectite clays could be rationalized by the use of ab initio calculations of the binding of the volatiles on the surface of a sodium smectite clay. At low coverage, the differences could be attributed to their differences in binding energies with the clay counter ions, which could be explained in terms of the chemical nature of the interacting species. At high coverage, the differences could be related to the properties of the bulk liquids, because high coverage corresponds to high activities of the volatiles on the clay samples. The differences observed in the measured rates of release of the volatiles from the smectite clays could be rationalized by the use of finite element calculations and the semi-analytical model developed for the bulk liquids, taking as input parameters the measured equilibrium desorption isotherms, that determines the volatile gas phase concentration at the surface of the sample. The slower rates of release of ethanol and ethyl acetate from the smectite clay, compared to water, could be explained from their differences in their equilibrium desorption isotherms from the clay, though diffusion effects also played a minor role. However, the slower rate of release of toluene, compared to water, could only be explained by its slow diffusion in the smectite clay. Diffusion effects were found to be enhanced with an increase in the size of the aggregate, such that, in addition to the binding strength of the counter ions, the swelling capacity of the clay, allowing the clay to clump, was also found to be an important factor for the controlled release of volatiles from the clay. Ion-exchange of the clay with lithium cations was found to be optimal in terms of ionic strength of the cation and swelling capacity, as compared to other metallic cations and small organic cations. This clay modification was tested on the release of malodorous compounds determined in cat urine by mass spectrometry and sniffing experiments. The smectite clay modified by ion-exchange with lithium cations showed improved controlled release of volatile organic compounds found in cat urine, compared to the sodium smectite, and improved clumping capacity compared to the calcium smectite. The calcium smectite clay used by Nestlé Purina could be modified by column exchange with lithium cations in high quantities. However the potential toxicity of the lithium is unknown. Additional information was obtained by ab initio calculations on the interaction of water with the smectite. The calculations of the binding of one to six water molecules on the surface of a sodium clay showed the water molecules bind to the cation and form hydrogen bonds between each other and with the clay oxygens, leading to the formation of different configurations for the same number of water molecules on the clay surface. With an increase in the number of water molecules on the surface of the clay, intermolecular interaction between the water molecules became gradually more important than sodium-water interactions. The binding energy was also found to be dependent on the location of the charge substitution in the clay (tetrahedral or octahedral) and to the formation of hydrogen bonds between the water molecule and a basal oxygen linked to a tetrahedral replacement. With an increase in the number of water molecules bound to one counter ion, the sodium counter ion was found to move away from the clay surface. The dynamic calculations of water molecules in the clay interlayer showed both the water molecules and the cations to be mobile. The average diffusion coefficient of the water molecules was reduced compared to that of bulk water, due to restrictions imposed by both the clay platelets and the counter ions. The structure of the water molecules showed to be similar to that of bulk water with a preferential orientation of their hydroxyl group almost perpendicular to the clay surface. The difference in the state of order of the water molecules around the sodium cation was interpreted as a symptom of hysteresis as observed in the measured equilibrium sorption isotherm.

EPFL2009

Metal Clusters on Surfaces

Michael Hugentobler

The morphology, stability and reactivity of nanometre-sized particles (Au, Pt) is investigated on three different supports. Gold nanoparticles with a diameter comprised between 4 and 6nm are stabilized in nanosized pits of well defined depth in highly oriented pyrolytic graphite (HOPG). These pits are produced by creation of artificial defects, followed by etching under a controlled oxygen atmosphere. At low Au coverage, clusters are found on the edges of the hexagonal pits maximizing the contact to dangling bonds on graphite multisteps. Larger coverage results in Au beads of well defined shape and with a constant bead density per unit length. Most remarkable is the stability of these nanostructures under ambient conditions. Temperatures as high as 650 K do not alter the morphology of the gold clusters. Above 700 K, the gold clusters catalyze the etching of graphite layers when heated under atmospheric conditions. The depth of the etching channel is defined by the depth of the pit and the channel width can be controlled by the cluster size. Hydrogen, oxygen and water are dissociated on the Au clusters at low temperatures (> 400 K). CO and CO2 is produced above 700 K which confirms the etching of the multilayer graphene sheets. The so-called water-gas shift reaction has been observed where oxygen and water dissociate on the Au clusters at low temperatures. Thermal desorption spectroscopy (TDS) measurements on the Pt=YSZ system have been performed and the desorption peaks of CO and O2 have been evidenced. Carbon monoxide desorbs from the surface at a temperature of 535 K and oxygen at 705 K, in good agreement with values reported in literature. Catalytic measurements have confirmed the known strong catalytic activity of Pt. The CO oxidation takes place at temperatures between 450 K and 700 K. The adsorption and desorption temperature of the two educts have been confirmed during the catalysis measurements. Size-selected Aun clusters (n = 5, 7) are deposited from the gas phase at room temperature with well-controlled kinetic energy on rutile TiO2 . Two different surface reconstructions have been used as support: TiO2(110) – (1 × 1) and TiO2(110) – (1 × 2). Systematic studies of morphology, stability and adsorption sites have been performed using scanning tunnelling microscopy (STM). The mean size of the clusters is maintained during deposition at a kinetic energy of 7.1eV per atom, indicating that there is no fragmentation. Au clusters are not stable neither on the TiO2(110) – (1 × 1) nor on the TiO2(110) – (1 × 2) reconstruction during annealing of the substrate to 800 K. A progressive transition into a pronounced 3-d formation is observed. The size distribution of the particles is small and the growth mode (Ostwald ripening) independent of the surface reconstruction. The higher corrugation of the TiO2(110) – (1 × 2) reconstruction results in a pronounced stabilization of the clusters on terraces, in contrast to the TiO2(110) – (1 × 1) reconstruction where 90% of the clusters are found on step edges at elevated temperatures. The catalytic activity of the Au clusters sets in during the 2d - 3d transition. These measurements have to be confirmed and a full understanding of the process should be found.