A fast XPS study of propene decomposition over clean and sulphated Pt{111}

The thermal decomposition of propene over clean and sulphate precovered Pt{111} has been followed by Fast XPS. The saturation propene coverage over the clean surface is 0.21 mL at 90 K. Propene is stable up to 200 K, above which molecular desorption and dehydrogenation result in the formation of a stable propylidyne intermediate adlayer at 300 K. Propylidyne decomposes above 400 K eventually forming graphitic carbon above 800 K. Preadsorbed surface sulphate promotes room temperature propene combustion associated with the decomposition of a thermally unstable alkyl--sulphate complex. Propylidyne also forms as on clean Pt{111}, but is less reactive, its decomposition above 450 K triggering partial oxidation with residual surface oxygen to liberate gas phase CO.

Effect of methane concentration in hydrogen plasma on hydrogen impurity incorporation in thick large-grained polycrystalline diamond films

We investigate the impact of methane concentration in hydrogen plasma on the growth of large-grained polycrystalline diamond (PCD) films and its hydrogen impurity incorporation. The diamond samples were produced using high CH4 concentration in H2 plasma and high power up to 4350 W and high pressure (either 105 or 110 Torr) in a microwave plasma chemical vapor deposition (MPCVD) system. The thickness of the free-standing diamond films varies from 165 µm to 430 µm. Scanning electron microscopy (SEM), micro-Raman spectroscopy and Fourier-transform infrared (FTIR) spectroscopy were used to characterize the morphology, crystalline and optical quality of the diamond samples, and bonded hydrogen impurity in the diamond films, respectively. Under the conditions employed here, when methane concentration in the gas phase increases from 3.75% to 7.5%, the growth rate of the PCD films rises from around 3.0 µm/h up to 8.5 µm/h, and the optical active bonded hydrogen impurity content also increases more than one times, especially the two CVD diamond specific H related infrared absorption peaks at 2818 and 2828 cm−1 rise strongly; while the crystalline and optical quality of the MCD films decreases significantly, namely structural defects and non-diamond carbon phase content also increases a lot with increasing of methane concentration. Based on the results, the relationship between methane concentration and diamond growth rate and hydrogen impurity incorporation including the form of bonded infrared active hydrogen impurity in CVD diamonds was analyzed and discussed. The effect of substrate temperature on diamond growth was also briefly discussed. The experimental findings indicate that bonded hydrogen impurity in CVD diamond films mainly comes from methane rather than hydrogen in the gas source, and thus can provide experimental evidence for the theoretical study of the standard methyl species dominated growth mechanism of CVD diamonds grown with methane/hydrogen mixtures.

A new regime for high rate growth of nanocrystalline diamond films using high power and CH4/H2/N2/O2 plasma

In this work, we report high growth rate of nanocrystalline diamond (NCD) films on silicon wafers of 2 inches in diameter using a new growth regime, which employs high power and CH4/H2/N2/O2 plasma using a 5 kW MPCVD system. This is distinct from the commonly used hydrogen-poor Ar/CH4 chemistries for NCD growth. Upon rising microwave power from 2000 W to 3200 W, the growth rate of the NCD films increases from 0.3 to 3.4 μm/h, namely one order of magnitude enhancement on the growth rate was achieved at high microwave power. The morphology, grain size, microstructure, orientation or texture, and crystalline quality of the NCD samples were characterized by scanning electron microscopy (SEM), atomic force microscopy (AFM), X-ray diffraction, and micro-Raman spectroscopy. The combined effect of nitrogen addition, microwave power, and temperature on NCD growth is discussed from the point view of gas phase chemistry and surface reactions. © 2011 Elsevier B.V. All rights reserved.

Catalyst studies on the ring opening of tetrahydrofuran-dimethanol to 1,2,6-hexanetriol

The metal catalyzed hydrogenolysis of the biomass-derived THF-dimethanol to 1,2,6-hexanetriol using heterogeneous catalysts was investigated. Bimetallic Rh-Re catalysts (4 wt% Rh and a Re/Rh (mol. ratio of 0.5) on a silica support gave the best performance and 1,2,6-hexanetriol was obtained in 84% selectivity at 31% conversion (120 C, 80 bar, 4 h); the selectivity reaches a maximum of 92% at 80 C. The product distribution at prolonged reaction times or higher temperatures or both shows the formation of diols and mono-alcohols, indicating that the 1,2,6-hexanetriol is prone to subsequent hydrodeoxygenation reactions. Different silica supports were investigated and optimal results were obtained with an amorphous silica featuring an intermediate surface area and an average mesopore size of about 6 nm. TPR and XPS surface analysis support the presence of mixed Rh and Re particles. The redox Reδ+/ReTotal surface ratio correlates with the conversion in a volcano type dependency. Both gas phase as well as Rh200Re1OH cluster DFT calculations support an acid-metal bifunctional mechanism and explain the products distribution. © 2013 Elsevier B.V. All rights reserved.

High pressure studies on group VI metal hexacarbonyl molecular solids

Group VI metal hexacarbonyls, M(CO)6 (M = Cr, Mo and W), are of extreme importance as catalysts in industry and also of fundamental interest due to the established charge transfer mechanism between the carbon monoxide and the metal. They condense to molecular solids at ambient conditions retaining the octahedral (Oh) symmetry of gas phase and have been extensively investigated by previous workers to understand their fundamental chemical bonding and possible industrial applications. However little is known about their behavior at high pressures which is the focus of this dissertation. Metal hexacarbonyls were subjected to high pressures in Diamond-Anvil cells to understand the pressure effect on chemical bonding using Raman scattering in situ. The high-pressure results on each of the three metal hexacarbonyls are presented and are followed by a critical analysis of the entire family. The Raman study was conducted at pressures up to 45 GPa and X-ray up to 58 GPa. This is followed by a discussion on infra red spectra in conjunction with Raman and X-ray analysis to provide a rationale for polymerization. Finally the probable synthesis of extremely reactive species under high-pressures and as identified via Raman is discussed. The high-pressure Raman scattering, up to 30 GPa, demonstrated the absence of Π-backbonding. The disappearance of parental Raman spectra for (M = Cr, Mo and W) at 29.6, 23.3 and 22.2 GPa respectively was attributed to the total collapse of the Oh symmetry. This collapse under high-pressure lead to metal-mediated polymeric phase characterized by Raman active δ(OCO) feature, originating from intermolecular vibrational coupling in the parent sample. Further increase in pressures up to 45 GPa, did not affect this feature. The pressure quenched Raman spectra, revealed various chemical groups non-characteristic of the parent sample and adsorption of CO in addition to the characteristic δ(OCO) feature. The thus recorded Raman, complemented with the far and mid-infrared pressure quenched spectra, reveal the formation of novel metal-mediated polymers. The X-ray diffraction on W(CO)6 up to 58 GPa revealed the generation of amorphous polymeric pattern which was retained back to ambient conditions.

Emergence and Fate of Siloxanes in Waste Streams: Release Mechanisms, Partitioning and Persistence in Three Environmental Compartments

Siloxanes are widely used in personal care and industrial products due to their low surface tension, thermal stability, antimicrobial and hydrophobic properties, among other characteristics. Volatile methyl siloxanes (VMS) have been detected both in landfill gas and biogas from anaerobic digesters at wastewater treatment plants. As a result, they are released to gas phase during waste decomposition and wastewater treatment. During transformation processes of digester or landfill gas to energy, siloxanes are converted to silicon oxides, leaving abrasive deposits on engine components. These deposits cause increased maintenance costs and in some cases complete engine overhauls become necessary. The objectives of this study were to compare the VMS types and levels present in biogas generated in the anaerobic digesters and landfills and evaluate the energetics of siloxane transformations under anaerobic conditions. Siloxane emissions, resulting from disposal of silicone-based materials, are expected to increase by 29% within the next 10 years. Estimated concentrations and the risk factors of exposure to siloxanes were evaluated based on the initial concentrations, partitioning characteristics and persistence. It was determined that D4 has the highest risk factor associated to bioaccumulation in liquid and solid phase, whereas D5 was highest in gas phase. Additionally, as siloxanes are combusted, the particle size range causes them to be potentially hazardous to human health. When inhaled, they may affix onto the alveoli of the lungs and may lead to development of silicosis. Siloxane-based COD-loading was evaluated and determined to be an insignificant factor concerning COD limits in wastewater. Removal of siloxane compounds is recommended prior to land application of biosolids or combustion of biogas. A comparison of estimated costs was made between maintenance practices for removal of siloxane deposits and installation/operation of fixed-bed carbon absorption systems. In the majority of cases, the installation of fixed-bed adsorption systems would not be a feasible option for the sole purpose of siloxane removal. However they may be utilized to remove additional compounds simultaneously.

Generalized atoms in molecules approach (GAIM approach)

This thesis involves two parts. The first is a new-proposed theoretical approach called generalized atoms in molecules (GAIM). The second is a computational study on the deamination reaction of adenine with OH⁻/nH₂O (n=0, 1, 2, 3) and 3H₂O. The GAIM approach aims to solve the energy of each atom variationally in the first step and then to build the energy of a molecule from each atom. Thus the energy of a diatomic molecule (A-B) is formulated as a sum of its atomic energies, EA and EB. Each of these atomic energies is expressed as, EA = Hᴬ + Vₑₑᴬᴬ + 1/2Vₑₑᴬ<>ᴮ EB = Hᴮ + Vₑₑᴮᴮ + 1/2Vₑₑᴬ<>ᴮ where; Hᴬ and Hᴮ are the kinetic and nuclear attraction energy of electrons of atoms A and B, respectively; Vₑₑᴬᴬ and Vₑₑᴮᴮ are the interaction energy between the electrons on atoms A and B, respectively; and Vₑₑᴬ<>ᴮ is the interaction energy between the electrons of atom A with the electrons of atom B. The energy of the molecule is then minimized subject to the following constraint, |ρA(r)dr + |ρB(r)dr = N where ρA(r) and ρB(r) are the electron densities of atoms A and B, respectively, and N is the number of electrons. The initial testing of the performance of GAIM was done through calculating dissociation curves for H₂, LiH, Li₂, BH, HF, HCl, N₂, F₂, and Cl₂. The numerical results show that GAIM performs very well with H₂, LiH, Li₂, BH, HF, and HCl. GAIM shows convergence problems with N₂, F₂, and Cl₂ due to difficulties in reordering the degenerate atomic orbitals Pₓ, Py, and Pz in N, F, and Cl atoms. Further work for the development of GAIM is required. Deamination of adenine results in one of several forms of premutagenic lesions occurring in DNA. In this thesis, mechanisms for the deamination reaction of adenine with OH⁻/nH₂O, (n = 0, 1, 2, 3) and 3H₂O were investigated. HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), and B3LYP/6-31+G(d) levels of theory were employed to optimize all the geometries. Energies were calculated at the G3MP2B3 and CBS-QB3 levels of theory. The effect of solvent (water) was computed using the polarizable continuum model (PCM). Intrinsic reaction coordinate (IRC) calculations were performed for all transition states. Five pathways were investigated for the deamination reaction of adenine with OH⁻/nH₂O and 3H₂O. The first four pathways (A-D) begin with by deprotonation at the amino group of adenine by OH⁻, while pathway E is initiated by tautomerization of adenine. For all pathways, the next two steps involve the formation of a tetrahedral intermediate followed by dissociation to yield products via a 1,3-hydrogen shift. Deamination with a single OH⁻ has a high activation barrier (190 kJ mol⁻¹ using G3MP2B3 level) for the rate-determining step. Addition of one water molecule reduces this barrier by 68 kJ mol⁻¹ calculated at G3MP2B3 level. Adding more water molecules decreases the overall activation energy of the reaction, but the effect becomes smaller with each additional water molecule. The most plausible mechanism is pathway E, the deamination reaction of adenine with 3H₂O, which has an overall G3MP2B3 activation energy of 139 and 137 kJ mol⁻¹ in the gas phase and PCM, respectively. This barrier is lower than that for the deamination with OH⁻/3H₂O by 6 and 2 kJ mol⁻¹ in the gas phase and PCM, respectively.

Evaporation maps for ternary non-ideal liquid mixtures

This thesis deals with the evaporation of non-ideal liquid mixtures using a multicomponent mass transfer approach. It develops the concept of evaporation maps as a convenient way of representing the dynamic composition changes of ternary mixtures during an evaporation process. Evaporation maps represent the residual composition of evaporating ternary non-ideal mixtures over the full range of composition, and are analogous to the commonly-used residue curve maps of simple distillation processes. The evaporation process initially considered in this work involves gas-phase limited evaporation from a liquid or wetted-solid surface, over which a gas flows at known conditions. Evaporation may occur into a pure inert gas, or into one pre-loaded with a known fraction of one of the ternary components. To explore multicomponent masstransfer effects, a model is developed that uses an exact solution to the Maxwell-Stefan equations for mass transfer in the gas film, with a lumped approach applied to the liquid phase. Solutions to the evaporation model take the form of trajectories in temperaturecomposition space, which are then projected onto a ternary diagram to form the map. Novel algorithms are developed for computation of pseudo-azeotropes in the evaporating mixture, and for calculation of the multicomponent wet-bulb temperature at a given liquid composition. A numerical continuation method is used to track the bifurcations which occur in the evaporation maps, where the composition of one component of the pre-loaded gas is the bifurcation parameter. The bifurcation diagrams can in principle be used to determine the required gas composition to produce a specific terminal composition in the liquid. A simple homotopy method is developed to track the locations of the various possible pseudo-azeotropes in the mixture. The stability of pseudo-azeotropes in the gas-phase limited case is examined using a linearized analysis of the governing equations. Algorithms for the calculation of separation boundaries in the evaporation maps are developed using an optimization-based method, as well as a method employing eigenvectors derived from the linearized analysis. The flexure of the wet-bulb temperature surface is explored, and it is shown how evaporation trajectories cross ridges and valleys, so that ridges and valleys of the surface do not coincide with separation boundaries. Finally, the assumption of gas-phase limited mass transfer is relaxed, by employing a model that includes diffusion in the liquid phase. A finite-volume method is used to solve the system of partial differential equations that results. The evaporation trajectories for the distributed model reduce to those of the lumped (gas-phase limited) model as the diffusivity in the liquid increases; under the same gas-phase conditions the permissible terminal compositions of the distributed and lumped models are the same.

Molecular composition of organic aerosols in central Amazonia: an ultra-high-resolution mass spectrometry study

The Amazon Basin plays key role in atmospheric chemistry, biodiversity and climate change. In this study we applied nanoelectrospray (nanoESI) ultra-high-resolution mass spectrometry (UHRMS) for the analysis of the organic fraction of PM2.5 aerosol samples collected during dry and wet seasons at a site in central Amazonia receiving background air masses, biomass burning and urban pollution. Comprehensive mass spectral data evaluation methods (e.g. Kendrick mass defect, Van Krevelen diagrams, carbon oxidation state and aromaticity equivalent) were used to identify compound classes and mass distributions of the detected species. Nitrogen- and/or sulfur-containing organic species contributed up to 60 % of the total identified number of formulae. A large number of molecular formulae in organic aerosol (OA) were attributed to later-generation nitrogen- and sulfur-containing oxidation products, suggesting that OA composition is affected by biomass burning and other, potentially anthropogenic, sources. Isoprene-derived organosulfate (IEPOX-OS) was found to be the most dominant ion in most of the analysed samples and strongly followed the concentration trends of the gas-phase anthropogenic tracers confirming its mixed anthropogenic–biogenic origin. The presence of oxidised aromatic and nitro-aromatic compounds in the samples suggested a strong influence from biomass burning especially during the dry period. Aerosol samples from the dry period and under enhanced biomass burning conditions contained a large number of molecules with high carbon oxidation state and an increased number of aromatic compounds compared to that from the wet period. The results of this work demonstrate that the studied site is influenced not only by biogenic emissions from the forest but also by biomass burning and potentially other anthropogenic emissions from the neighbouring urban environments.

Reactions of atomic and molecular ions with acetone, 1,1,1- trifluoroacetone, and hexafluoroacetone: An investigation of the effects of molecular structure on the dynamics and kinetics of ion-molecule reactions

A selected ion flow tube study of the reactions of a series of gas-phase atomic cations (S⁺, Xe⁺, O⁺, Kr⁺, N⁺, Ar⁺ and Ne⁺) and molecular ions (SF _n⁺ (n = 1-5), CF_n⁺ (n = 1-3), CF₂Cl⁺, H₃O⁺, NO⁺, N ₂O⁺, CO₂⁺, CO⁺, and N₂⁺) spanning a large range of recombination energies (6.3-21.6 eV), with acetone, 1,1,1-trifluoroacetone, and hexafluoroacetone has been undertaken with the objective of exploring the nature of the reaction ion chemistry as the methyl groups in acetone are substituted for CF₃. The reaction rate coefficients and product ion branching ratios for all 66 reactions, measured at 298 K, are reported. The experimental reaction rate coefficients are compared to theoretically calculated collisional values. Several distinct reaction processes were observed among the large number of reactions studied, including charge transfer (non-dissociative and dissociative), abstraction, ion-molecule associations and, in the case of the reactions involving the reagent ion H₃O⁺, proton transfer. 

Ionic liquids as modulators of fragrance release in consumer goods

When organic esters or alcohols were dissolved in each of three novel ionic liquids (which have no effective vapour pressure), the vapour–liquid equilibria (as measured by infrared spectroscopy of the gas phase) revealed significant positive deviation from Raoult’s law for a wide range of perfume raw materials. The addition of water amplified the repulsive effect of the ionic liquid matrix, and this was exemplified by a series of ternary phase diagrams

A New Technique for Studying Vapour–liquid Equilibria of Multi-Component Systems

A Fourier transform infrared gas-phase method is described herein and capable of deriving the vapour pressure of each pure component of a poorly volatile mixture and determining the relative vapour phase composition for each system. The performance of the present method has been validated using two standards (naphthalene and ferrocene), and a Raoult’s plot surface of a ternary system is reported as proof-of-principle. This technique is ideal for studying solutions comprising two, three, or more organic compounds dissolved in ionic liquids as they have no measurable vapour pressures.

Catalizzatori polifunzionali per la sintesi di Metil Isobutil Chetone

In this work, we have examined the activity and selectivity of new catalysts for the single-stage production of methyl isobutyl ketone (MIBK, 4- methyl-2-pentanone) from acetone (both in liquid and gas phase), using a fixed bed reactor operated in the temperature range between 373 and 473 K. The main reaction pathways for the synthesis of MIBK from acetone are given in Fig.1. The first step is the self condensation of acetone to diacetone alcohol (DAA, 4-hydroxy-4-methyl-2-pentanone); the second step is the dehydration of DAA to mesityl oxide (MO, 4-methyl-3-penten-2-one); the final step is the selective hydrogenation of the carbon–carbon double bond of MO to form MIBK. The most commonly observed side reactions are over-condensations and unselective hydrogenations (also shown in Fig.1). Two types of catalysts were studied: i)Pd supported on MgO-SiO2 mixed oxides with ratio of Mg to Si, synthetized using Ohnishi’s method and ii)Pd supported on alumina doped with 5% or 10% of MgO. The different Mg-Si and Mg-Al catalysts were characterized by different techniques (XRD, BET, SEM, NH3-TPD and CO2-TPD) and tested under different conditions in the condensation of acetone to diacetone alcohol and its dehydration to mesityl oxide to enhance the activity. Palladium was chosen as metal component, and its hydrogenation activity was studied. A low hydrogenation activity negatively affects the acetone conversion and promotes the production of mesityl oxide. Hydrogenation conditions being too severe may favor the unwanted hydrogenation of acetone to 2-propanol and of MIBK to methyl isobutyl carbinol (MIBC, 4-methyl-2-pentanol) but this effect is less detrimental to the MIBK selectivity than an unsufficient hydrogenation activity.

The FIREDASS (Fire Detection and Suppression Simulation) Model

The FIREDASS (FIRE Detection And Suppression Simulation) project is concerned with the development of fine water mist systems as a possible replacement for the halon fire suppression system currently used in aircraft cargo holds. The project is funded by the European Commission, under the BRITE EURAM programme. The FIREDASS consortium is made up of a combination of Industrial, Academic, Research and Regulatory partners. As part of this programme of work, a computational model has been developed to help engineers optimise the design of the water mist suppression system. This computational model is based on Computational Fluid Dynamics (CFD) and is composed of the following components: fire model; mist model; two-phase radiation model; suppression model and detector/activation model. The fire model - developed by the University of Greenwich - uses prescribed release rates for heat and gaseous combustion products to represent the fire load. Typical release rates have been determined through experimentation conducted by SINTEF. The mist model - developed by the University of Greenwich - is a Lagrangian particle tracking procedure that is fully coupled to both the gas phase and the radiation field. The radiation model - developed by the National Technical University of Athens - is described using a six-flux radiation model. The suppression model - developed by SINTEF and the University of Greenwich - is based on an extinguishment crietrion that relies on oxygen concentration and temperature. The detector/ activation model - developed by Cerberus - allows the configuration of many different detector and mist configurations to be tested within the computational model. These sub-models have been integrated by the University of Greenwich into the FIREDASS software package. The model has been validated using data from the SINTEF/GEC test campaigns and it has been found that the computational model gives good agreement with these experimental results. The best agreement is obtained at the ceiling which is where the detectors and misting nozzles would be located in a real system. In this paper the model is briefly described and some results from the validation of the fire and mist model are presented.

Séquestration géologique du CO2 par carbonatation minérale dans les résidus miniers

La carbonatation minérale dans les résidus miniers est un moyen sûr et permanent de séquestrer le CO2 atmosphérique. C’est un processus naturel et passif qui ne nécessite aucun traitement particulier et donc avantageux d’un point de vue économique. Bien que la quantité de CO2 qu’il soit possible de séquestrer selon ce processus est faible à l’échelle globale, dans le cadre d’un marché du carbone, les entreprises minières pourraient obtenir des crédits et ainsi revaloriser leurs résidus. À l’heure actuelle, il y a peu d’informations pour quantifier le potentiel de séquestration du CO2 de façon naturelle et passive dans les piles de résidus miniers. Il est donc nécessaire d’étudier le phénomène pour comprendre comment évolue la réaction à travers le temps et estimer la quantité de CO2 qui peut être séquestrée naturellement dans les piles de résidus. Plusieurs travaux de recherche se sont intéressés aux résidus miniers de Thetford Mines (Québec, Canada), avec une approche principalement expérimentale en laboratoire. Ces travaux ont permis d’améliorer la compréhension du processus de carbonatation, mais ils nécessitent une validation à plus grande échelle sous des conditions atmosphériques réelles. L’objectif général de cette étude est de quantifier le processus de carbonatation minérale des résidus miniers sous des conditions naturelles, afin d’estimer la quantité de CO2 pouvant être piégée par ce processus. La méthodologie utilisée repose sur la construction de deux parcelles expérimentales de résidus miniers situées dans l’enceinte de la mine Black Lake (Thetford Mines). Les résidus miniers sont principalement constitués de grains et de fibres de chrysotile et lizardite mal triés, avec de petites quantités d’antigorite, de brucite et de magnétite. Des observations spatiales et temporelles ont été effectuées dans les parcelles concernant la composition et la pression des gaz, la température des résidus, la teneur en eau volumique, la composition minérale des résidus ainsi que la chimie de l’eau des précipitations et des lixiviats provenant des parcelles. Ces travaux ont permis d’observer un appauvrissement notable du CO2 dans les gaz des parcelles (< 50 ppm) ainsi que la précipitation d’hydromagnésite dans les résidus, ce qui suggère que la carbonatation minérale naturelle et passive est un processus potentiellement important dans les résidus miniers. Après 4 ans d’observations, le taux de séquestration du CO2 dans les parcelles expérimentales a été estimé entre 3,5 et 4 kg/m3/an. Ces observations ont permis de développer un modèle conceptuel de la carbonatation minérale naturelle et passive dans les parcelles expérimentales. Dans ce modèle conceptuel, le CO2 atmosphérique (~ 400 ppm) se dissout dans l’eau hygroscopique contenue dans les parcelles, où l’altération des silicates de magnésium forme des carbonates de magnésium. La saturation en eau dans les cellules est relativement stable dans le temps et varie entre 0,4 et 0,65, ce qui est plus élevé que les valeurs de saturation optimales proposées dans la littérature, réduisant ainsi le transport de CO2 dans la zone non saturée. Les concentrations de CO2 en phase gazeuse, ainsi que des mesures de la vitesse d’écoulement du gaz dans les cellules suggèrent que la réaction est plus active près de la surface et que la diffusion du CO2 est le mécanisme de transport dominant dans les résidus. Un modèle numérique a été utilisé pour simuler ces processus couplés et valider le modèle conceptuel avec les observations de terrain. Le modèle de transport réactif multiphase et multicomposant MIN3P a été utilisé pour réaliser des simulations en 1D qui comprennent l’infiltration d’eau à travers le milieu partiellement saturé, la diffusion du gaz, et le transport de masse réactif par advection et dispersion. Même si les écoulements et le contenu du lixivat simulés sont assez proches des observations de terrain, le taux de séquestration simulé est 22 fois plus faible que celui mesuré. Dans les simulations, les carbonates précipitent principalement dans la partie supérieure de la parcelle, près de la surface, alors qu’ils ont été observés dans toute la parcelle. Cette différence importante pourrait être expliquée par un apport insuffisant de CO2 dans la parcelle, qui serait le facteur limitant la carbonatation. En effet, l’advection des gaz n’a pas été considérée dans les simulations et seule la diffusion moléculaire a été simulée. En effet, la mobilité des gaz engendrée par les fluctuations de pression barométrique et l’infiltration de l’eau, ainsi que l’effet du vent doivent jouer un rôle conséquent pour alimenter les parcelles en CO2.