Pb Removal Efficiency by Calcium Carbonates: Biogenic versus Abiogenic Materials.

The sorption of heavy metals on mineral surfaces plays a key role in controlling the fate and bioavailability of harmful elements through dissolution-precipitation reactions. Here, we investigate the efficiency of Pb removal from highly contaminated waters by two calcium carbonate hard tissues, scallop shells (up to 99.9 mol %; -biocalcite) and cuttlefish bones (up to 90.0 mol %; bioaragonite), which template the precipitation of the highly insoluble mineral cerussite (PbCO3). The experiments show that both biomaterials are about five times more effective Pb scavengers (5 mmol of cerussite precipitated/g sample) than their inorganic counterparts (∼1 mmol/g). We relate this enhanced Pb scavenging capacity of biocarbonates to their composite organic-inorganic nature, which modulates their specific nano- and microstructural features and defines their larger surface areas, solubility, and reactivity compared to those of their inorganic counterparts. The oriented growth of cerussite progressively passivates the bioaragonite surface, reducing its long-term Pb scavenging capacity. In contrast, the randomly oriented growth of cerussite crystals on biocalcite prevents surface passivation and explains why biocalcite outperforms bioaragonite as a long-term Pb scavenger. The use of biocarbonates could be a key for designing more efficient decontamination strategies for heavy metal-polluted waters.

Thermodynamic model of MSWI flue gas cooling path: Effect of flue gas composition on heavy metal binding forms.

In the context of circular economy and heavy metal (HM) recovery from municipal solid waste incineration (MSWI) fly ash (FA), detailed knowledge of HM binding forms is required for achieving higher extraction rates. The FA mineralogy is still poorly understood due to its low grain size and low metal concentration. To investigate the HM binding forms, a sophisticated thermodynamic reactive transport model was developed to simulate ash-forming processes. The stability of different binding forms was investigated at different flue gas conditions (varying ratios of HCl, SO2, O2) by simulating the gas cooling path in closed system and dynamic open system, where the gas composition is changing upon cooling due to precipitation of solids. The simulations predict that at flue gas conditions of molar ratio S/Cl < 1, Cu and Zn precipitate as oxides (and Zn silicates) at approximately 650°C. At temperatures <300°C, Zn, Cu, Pb and Cd are predicted to precipitate as easily soluble chlorides. In flue gas with molar ratio S/Cl > 1, the HM precipitate as less soluble sulphates. The results indicate that the less soluble HM fraction in the electrostatic precipitator ash represent oxides and silicates that formed in the boiler section but were transported to the electrostatic precipitator. The model provides insight into the physical-chemical processes controlling the metal accumulation in the flue gas and FA during the cooling of the flue gas. The obtained data serve as valuable basis for improving metal recovery from MSWI FA.

Stability and Speciation of Hydrated Magnetite {111} Surfaces from Ab Initio Simulations with Relevance for Geochemical Redox Processes.

Magnetite is a common mixed Fe(II,III) iron oxide in mineral deposits and the product of (anaerobic) iron corrosion. In various Earth systems, magnetite surfaces participate in surface-mediated redox reactions. The reactivity and redox properties of the magnetite surface depend on the surface speciation, which varies with environmental conditions. In this study, Kohn-Sham density functional theory (DFT + U method) was used to examine the stability and speciation of the prevalent magnetite crystal face {111} in a wide range of pH and Eh conditions. The simulations reveal that the oxidation state and speciation of the surface depend strongly on imposed redox conditions and, in general, may differ from those of the bulk state. Corresponding predominant phase diagrams for the surface speciation and structure were calculated from first principles. Furthermore, classical molecular dynamics simulations were conducted investigating the mobility of water near the magnetite surface. The obtained knowledge of the surface structure and oxidation state of iron is essential for modeling retention of redox-sensitive nuclides.

Preceramic Polymers for Additive Manufacturing of Silicate Ceramics.

The utilization of preceramic polymers (PCPs) to produce both oxide and non-oxide ceramics has caught significant interest, owing to their exceptional characteristics. Diverse types of polymer-derived ceramics (PDCs) synthesized by using various PCPs have demonstrated remarkable characteristics such as exceptional thermal stability, resistance to corrosion and oxidation at elevated temperatures, biocompatibility, and notable dielectric properties, among others. The application of additive manufacturing techniques to produce PDCs opens up new opportunities for manufacturing complex and unconventional ceramic structures with complex designs that might be challenging or impossible to achieve using traditional manufacturing methods. This is particularly advantageous in industries like aerospace, automotive, and electronics. In this review, various categories of preceramic polymers employed in the synthesis of polymer-derived ceramics are discussed, with a particular focus on the utilization of polysiloxane and polysilsesquioxanes to generate silicate ceramics. Further, diverse additive manufacturing techniques adopted for the fabrication of polymer-derived silicate ceramics are described.

Oxidation State and Structure of Fe in Nontronite: From Oxidizing to Reducing Conditions.

The redox reaction between natural Fe-containing clay minerals and its sorbates is a fundamental process controlling the cycles of many elements such as carbon, nutrients, redox-sensitive metals, and metalloids (e.g., Co, Mn, As, Se), and inorganic as well as organic pollutants in Earth's critical zone. While the structure of natural clay minerals under oxic conditions is well-known, less is known about their behavior under anoxic and reducing conditions, thereby impeding a full understanding of the mechanisms of clay-driven reduction and oxidation (redox) reactions especially under reducing conditions. Here we investigate the structure of a ferruginous natural clay smectite, nontronite, under different redox conditions, and compare several methods for the determination of iron redox states. Iron in nontronite was gradually reduced chemically with the citrate-bicarbonate-dithionite (CBD) method. 57Fe Mössbauer spectrometry, X-ray photoelectron spectroscopy (XPS), X-ray absorption near edge structure (XANES) spectroscopy including its pre-edge, extended X-ray absorption fine structure (EXAFS) spectroscopy, and mediated electrochemical oxidation and reduction (MEO/MER) provided consistent Fe(II)/Fe(III) ratios. By combining X-ray diffraction (XRD) and transmission electron microscopy (TEM), we show that the long-range structure of nontronite at the highest obtained reduction degree of 44% Fe(II) is not different from that of fully oxidized nontronite except for a slight basal plane dissolution on the external surfaces. The short-range order probed by EXAFS spectroscopy suggests, however, an increasing structural disorder and Fe clustering with increasing reduction of structural Fe.

Diffusion and Gas Flow Dynamics in Partially Saturated Smectites.

Clays and clay rocks are considered good natural and engineered barriers for deep geological disposal of nuclear waste worldwide. Metal corrosion and organic waste degradation in underground repositories generate significant amounts of gas that should be able to migrate through the multibarrier system to avoid potential pressure buildup, which could be compromising the integrity of the barriers and host rocks. The gas is expected to accumulate in larger pores and eventually form an interconnected network. Under such conditions, the migration of gas molecules takes place both in pore water films and gas-filled macropores. Therefore, mass fluxes depend on the distribution of gas molecules between the water-rich and gas-rich phases and their mobility in both compartments. Classical molecular dynamics (MD) simulations were employed to investigate the mobilities of He, H2, CO2, Ar, and CH4 in a Na-montmorillonite mesopore as a function of the degree of saturation, as well as evaluate the hydrodynamic behavior of the pore fluid in partially saturated clays. The diffusivity of the gas molecules was determined by observing the asymptotic behavior of the mean square displacement in the gas-rich phase and at the gas-water interface. The partition coefficient and Gibbs free energy were analyzed to investigate the transfer of gas molecules between the gas-rich and water-rich phases by observing the molecular trajectories as they cross the vapor-liquid interface. The results revealed that the diffusion coefficient in the gas phase increased with increasing gas-filled pore width and converged asymptotically toward the diffusion coefficient in the bulk state. It could be shown that the diffusion coefficient of gas molecules dissolved in the water films remained constant as long as the interacting water surface was in the bulk-liquid-like phase. This behavior changes in very thin water films. It was observed that the partitioning coefficient of gas molecules at the solid-liquid interface is nearly the same as that in the bulk-liquid-like phase. Partitioning is observed to be strongly dependent on the temperature and gas molecular weights. In the second part of the study, nonequilibrium molecular dynamics (NEMD) simulations were performed to investigate the mobility of gases in pressure-driven decoupled gas-phase dynamics (DGPD) and coupled gas and water phase dynamics (CGWPD) in a partially saturated Na-montmorillonite slit mesopore. The dynamic viscosity of the gas phase was calculated from NEMD simulations and indicated that the viscosity of the gas phase was almost the same in both methods (DGPD and CGWPD). The average slip length for gas molecules at the gas-water interface was also calculated, revealing that the slip-free boundary condition assumed in continuum models is generally invalid for microfluidics and that a slip boundary condition exists at the microscale for specific surface interactions. Finally, a Bosanquet-type equation was developed to predict the diffusion coefficient and dynamic viscosity of gas as a function of the average pore width, gas mean-free path, geometric factor, and thickness of the adsorbed water film.

Classification of perovskite structural types with dynamical octahedral tilting.

Perovskites ABX3 with delocalized positions of the X atoms represent a distinct class of dynamically distorted structures with peculiar structural relations and physical properties. The delocalization originates from atoms crossing shallow barriers of the potential energy surface. Quantum mechanically, they can be treated similar to light atoms in diffusive states. Many of these perovskite structures are widely used functional materials thanks to their particular physical properties, such as superconductivity, ferroelectricity and photo-activity. A number of these properties are related to static or dynamic motion of octahedral units. Yet, a full understanding of the relationships between perovskite crystal structure, chemical bonding and physical properties is currently missing. Several studies indicate the existence of dynamic disorder generated by anharmonic motion of octahedral units, e.g. in halide perovskite structures. To simplify structural analysis of such systems we derive a set of space groups for simple perovskites ABX3 with dynamical octahedral tilting. The derived space groups extend the well established space group tables for static tiltings by Glazer [Acta Cryst. B (1972). 28, 3384-3392], Aleksandrov [Ferroelectrics (1976). 24, 801-805] and Howard & Stokes [Acta Cryst. B (1998). 54, 782-789]. Ubiquity of dynamical tilting is demonstrated by an analysis of the structural data for perovskites reported in recent scientific publications and the signature of dynamic tilting in the corresponding structures is discussed, which can be summarized as follows: (a) volume increase upon a lowering of temperature, (b) apparent distortion of octahedra (where Jahn-Teller distortions can be ruled out), (c) mismatch between observed instantaneous symmetry and average symmetry, (d) deviation of the experimental space group from the theoretically predicted structures for static tilting, (e) inconsistency of lattice parameters with those suggested by the theory of static tilts, and (f) large displacement parameters for atoms at the X and B sites. Finally, the possible influence of dynamic disorder on the physical properties of halide perovskites is discussed.

Pbx(OH)y cluster formation and anomalous thermal behaviour in STI framework-type zeolites.

For the first time, the structural investigation of a Pb-exchanged zeolite (Pb13.4(OH)10Al17.4Si54.6O144 ∙38H2O) with STI framework type, revealed a highly unusual and intriguing sudden volume increase under continuous heating. Understanding the fundamental mechanisms leading to such an unusual behaviour is essential for technological applications and interpretation of chemical bonding in zeolites. The dehydration was tracked in situ from 25 to 450 °C by single crystal X-ray diffraction, infrared and X-ray absorption spectroscopy. Further interpretation of the experimental observations was supported by ab initio molecular dynamics simulations. Initially, Pb-STI unit-cell volume contracts (ΔV = - 3.5%) from 25 to 100 °C. This agrees with the trend observed in STI zeolites. Surprisingly, at 125 °C, the framework expanded (ΔV = + 2%), adopting a configuration, which resembles that of the room temperature structure. Upon heating, the structure loses H2O but no de-hydroxylation occurred. The key mechanism leading to the sudden volume increase was found to be the formation of Pbx(OH)y clusters, which prevent the shrinking of the channels, rupture of the tetrahedral bonds and occlusion of the pores. This zeolite has therefore an increased thermal stability with respect to other STI metal-exchanged zeolites, with important consequences on its applications.

Role of lone-pair electron localization in temperature-induced phase transitions in mimetite.

The crystal structure of mimetite Pb5(AsO4)3Cl, a phosphate with apatite structure-type has been investigated in situ at 123, 173, 273, 288, 353 and 393 K. A careful inspection of the diffraction pattern and subsequent structure refinements indicated that mimetite transforms from the monoclinic to the hexagonal polymorph with increasing temperature. At 123 K, a monoclinic superstructure, mimetite-2M, with cell parameters a = 20.4487 (9),  b = 7.4362 (2), c = 20.4513 (9) Å, ß = 119.953 (6)°, V = 2694.5 (2) Å3 and space group P21 was observed. From 173 to 353 K, the reflections of the supercell were evident only along one direction of the corresponding hexagonal apatite-cell and the structure transforms to the polymorph mimetite-M with space group P21/b and unit-cell parameters a = 10.2378 (3), b = 20.4573 (7), c = 7.4457 (2) Å, ß = 120.039 (5)°, V = 1349.96 (9) Å3. Only at higher temperature, i.e. 393 K, does mimetite adopt the hexagonal space group P63/m characteristic of apatite structure-types. The role of the electron lone pairs of Pb atoms in the phase transition was investigated through the analysis of the electron localization function (ELF) calculated based on the DFT-geometry optimized structures of the three polymorphs. The changes in spatial distribution of the 6s2 electron density during the phase transitions were explored by means of the Wannier Function Centres (WFCs) derived from ab initio molecular dynamics trajectories. In the high-temperature hexagonal structure the 6s2 electrons are spherically symmetric relative to the position of Pb atoms. At low temperature the maximum of 6s2 electron density is displaced relative to the position of Pb atom contributing to the polar interaction in the monoclinic polymorphs.

A lab-on-a-chip approach integrating in-situ characterization and reactive transport modelling diagnostics to unravel (Ba,Sr)SO4 oscillatory zoning.

The co-precipitation of sulphate minerals such as celestine and barite is widely studied because their formation is ubiquitous in natural and anthropogenic systems. Co-precipitation in porous media results in crystallization of solid solutions yielding characteristics such as oscillatory zoning that are rarely observed in bulk solution or in batch experiments. In the past, the precipitation of compositionally-zoned (Ba,Sr)SO4 crystals was observed post-mortem in macroscopic silica gel counter-diffusion experiments. Their formation was originally explained by the difference in the solubility products of the end-members combined with diffusion-limited transport of solutes to the mineral-fluid interface, while a later study favored the idea of kinetically controlled reactions. With recent advances combining in-operando microfluidic experiments and reactive transport modelling, it is now possible to verify hypotheses on the driving forces of transport-coupled geochemical processes. We developed a "lab on a chip" experiment that enabled the systematic study of the nucleation and growth of oscillatory-zoned (Ba,Sr)SO4 crystals in a microfluidic reactor. The compositions of the solid solutions were determined by in-situ Raman spectroscopy. Our investigation shows (1) that the composition of the nucleating phases can be approximated using classical nucleation theory, (2) that the oscillatory zoning is not solely controlled by the limited diffusional transport of solutes, and (3) that nucleation kinetics plays a major role in the switch between different stoichiometric compositions. The zoning phenomena is governed by the complex interplay between the diffusion of reactants and the crystallization kinetics as well as other factors, e.g. surface tension and lattice mismatch.

Deciphering the molecular mechanism of water boiling at heterogeneous interfaces.

Water boiling control evolution of natural geothermal systems is widely exploited in industrial processes due to the unique non-linear thermophysical behavior. Even though the properties of water both in the liquid and gas state have been extensively studied experimentally and by numerical simulations, there is still a fundamental knowledge gap in understanding the mechanism of the heterogeneous nucleate boiling controlling evaporation and condensation. In this study, the molecular mechanism of bubble nucleation at the hydrophilic and hydrophobic solid-water interface was determined by performing unbiased molecular dynamics simulations using the transition path sampling scheme. Analyzing the liquid to vapor transition path, the initiation of small void cavities (vapor bubbles nuclei) and their subsequent merging mechanism, leading to successively growing vacuum domains (vapor phase), has been elucidated. The molecular mechanism and the boiling nucleation sites' location are strongly dependent on the solid surface hydrophobicity and hydrophilicity. Then simulations reveal the impact of the surface functionality on the adsorbed thin water molecules film structuring and the location of high probability nucleation sites. Our findings provide molecular-scale insights into the computational aided design of new novel materials for more efficient heat removal and rationalizing the damage mechanisms.

Coupling of multiscale lattice Boltzmann discrete-element method for reactive particle fluid flows.

Reactive particulate systems are of prime importance in varieties of practical applications in process engineering. As an example this study considers extraction of phosphorous from waste water by calcium silicate hydrate particles in the P-RoC process. For such systems modeling has a large potential to help to optimize process conditions, e.g., particle-size distributions or volume flows. The goal of this study is to present a new generic modeling framework to capture relevant aspects of reactive particle fluid flows using combined lattice Boltzmann method and discrete-element method. The model developed is Euler-Lagrange scheme which consist of three-components viz., a fluid phase, a dissolved reactive substance, and suspended particles. The fluid flow and reactive mass transport are described in a continuum framework using volume-averaged Navier-Stokes and volume-averaged advection-diffusion-reaction equations, respectively, and solved using lattice Boltzmann methods. The volume-averaging procedure ensures correctness in coupling between fluid flow, reactive mass transport, and particle motion. The developed model is validated through series of well-defined benchmarks. The benchmarks include the validation of the model with experimental data for the settling of a single particle in a cavity filled with water. The benchmark to validate the multi-scale reactive transport involves comparing the results with a resolved numerical simulation. These benchmarks also prove that the proposed model is grid convergent which has previously not been established for such coupled models. Finally, we demonstrate the applicability of our model by simulating a suspension of multiple particles in fluid with a dissolved reactive substance. Comparison of this coupled model is made with a one-way coupled simulation where the influence of particles on the fluid flow and the reactive solution transport is not considered. This elucidates the need for the two-way coupled model.

Anharmonic effects on the dynamics of solid aluminium fromab initiosimulations.

Two approaches to simulations of phonon properties of solids beyond the harmonic approximation, the self-consistentab initiolattice dynamics (SCAILD) and decoupled anharmonic mode approximation (DAMA) are critically benchmarked against each other and molecular dynamics simulations using a density-functional-theory description of electronic states, and compared to experimental data for fcc aluminium. The temperature-dependence of phonon dispersion and the phonon density-of-states, heat capacity, and the mean atomic displacement for fcc aluminium are examined with these approaches at ambient pressure. A comparison of results obtained with the harmonic approximation to the ones predicted by SCAILD and DAMA reveal a negligible anharmonic contribution to phonon frequencies, a small, but significant influence on heat capacity, and a strong effect on atomic mean-square displacement. The phase space accessed with SCAILD and DAMA is reduced relative to molecular and harmonic lattice dynamics simulations. In particular the DAMA results are in good agreement with displacement amplitudes determined by the Debye-Waller factor in x-ray diffraction experiments.

Structural Evolution of Ag-LEV Zeolite upon Heating: An In Situ Single-Crystal X-ray Diffraction (SC-XRD) and X-ray Absorption Spectroscopy (XAS) Study.

Ag-exchanged zeolites are known to have improved sorption and catalytic properties compared to the raw natural material. The application range of Ag-exchanged zeolites is linked to the stability of the exchanged form and its structural evolution at high temperatures. In this study, we investigated the thermal stability of a Ag-exchanged zeolite with an LEV framework type. The dehydration path was monitored in situ by single-crystal X-ray diffraction (SC-XRD) and X-ray absorption fine structure spectroscopy (XAFS). The experimental data were compared with those extrapolated from molecular dynamics (MD) trajectories. Our results showed that Ag-exchanged levyne (Ag-LEV) follows a different dehydration path compared to that of the natural levyne (Ca-LEV). Between 25 and 350 °C, the unit cell volume contraction was -4% with respect to that measured at room temperature. Upon dehydration, Ag-LEV transformed to the LEV B topology: such transformation is accompanied by the change from R 3̅ m to R 3m space group and by the onset of the rupture of one T-O-T connection at 250 °C. Differently from Ca-LEV, no additional change to LEV B' configuration was detected. XAFS analysis indicated that each Ag is approximately surrounded by four oxygen atoms between 2.15 and 2.40 Å. This local environment was maintained up to 650 °C. Weak Ag+1-Ag+1 interactions, detected in the whole investigated temperature range, are mainly ascribed to the presence of Ag-erionite (Ag-ERI) intergrown with Ag-LEV. No reduction to elemental Ag0 occurred upon heating.

Fundamental Research on Radiochemistry of Geological Nuclear Waste Disposal.

Currently, 5 · 1019 Bq of radioactive waste originating from the use of nuclear power for energy production, and medicine, industry and research, is maintained in Switzerland at intermediate storage facilities. Deep geological disposal of nuclear waste is considered as the most reliable and sustainable long-term solution worldwide. Alike the other European countries, the Swiss waste disposal concept embarks on the combination of engineered and geological barriers. The disposal cell is a complex geochemical system. The radionuclide mobility and consequently radiological impact depend not only on their chemical speciation but also on the background concentration of other stable nuclides and their behaviour in the natural environment. The safety assessment of the repository is thus a complex multidisciplinary problem requiring knowledge in chemical thermodynamics, structural chemistry, fluid dynamics, geo- and radiochemistry. Broad aspects of radionuclide thermodynamics and geochemistry are investigated in state-of-the-art radiochemical laboratories at the Paul Scherrer Institute. The research conducted over the last 30 years has resulted in a fundamental understanding of the radionuclides release, retention and transport mechanism in the repository system.

Iron Adsorption on Clays Inferred from Atomistic Simulations and X-ray Absorption Spectroscopy.

The atomistic level understanding of iron speciation and the probable oxidative behavior of iron (Feaq2+ → Fesurf3+) in clay minerals are fundamental for environmental geochemistry of redox reactions. Thermodynamic analyses of wet chemistry data suggest that iron adsorbs on the edge surfaces of clay minerals at distinct structural sites commonly referred as strong and weak sites (with high and low affinity, respectively). In this study, we applied ab initio molecular dynamics simulation to investigate the structure and the stability of the edge surfaces of trans- and cis-vacant montmorillonites. These structures were further used to evaluate the surface complexation energy and to calculate reference ab initio X-ray absorption spectra (XAS) for distinct inner-sphere complexes of iron. The combination of ab initio simulations and XAS allowed us to reveal the Fe-complexation mechanism and to quantify the Fe partitioning between the high and low affinity sites as a function of the oxidation state and loadings. Although iron is mostly present in the Fe3+ form, Fe2+ increasingly co-adsorbs at increasing loadings. Ab initio structure relaxations of several different clay structures with substituted Fe2+/Fe3+ in the bulk or at the surface site showed that the oxidative sorption of ferrous iron is an energetically favored process at several edge surfaces of the Fe-bearing montmorillonite.

Characterization of Structural Iron in Smectites - An Ab Initio Based X-ray Absorption Spectroscopy Study.

Fe-bearing clay minerals are abundant in argillaceous rocks as their redox-active structural iron may control the sorption mechanism of redox sensitive elements on the surface of clay minerals. The extent and efficiency of the redox reactions depend on the oxidation state (Fe2+/Fe3+ ratio) and structural distribution of the substituting cations in the TOT-layer of clay minerals. Even smectites with similar structure originating from different locations might have a distinct arrangement of isomorphic substitutions (e.g., individual iron or Fe-Fe pairs). In this study, the proportion of different iron distribution in Milos-, Wyoming-, and Texas-montmorillonite was determined by combining X-ray absorption spectroscopy (XAS) with ab initio calculations. The relaxed atomic structures of the smectite models with different arrangement of individual Fe atoms and Fe-Fe/Fe-Mg clusters served as the basis for the calculations of the XAS spectra. The combination of simulation results and measured Fe K-edge XAS spectra of Wyoming-, Milos- and Texas-montmorillonites suggested that iron is present as Fe3+ in the octahedral sheet. Fe3+ in Texas-montmorillonite has a tendency to form clusters, while no definitive statement about clustering or avoidance of Fe-Fe and Fe-Mg pairs can be made for Milos- and Wyoming-montmorillonite.

Combined XAFS Spectroscopy and Ab Initio Study on the Characterization of Iron Incorporation by Montmorillonite.

Iron occurs in clay minerals in both ferric and ferrous forms. Depending on its oxidation state and the environmental conditions, it can participate in redox reactions and influence the sorption processes at surfaces of clay minerals. Knowing the oxidation state and the preferential structural position of Fe2+ and Fe3+ is essential for the detailed understanding of the mechanism and kinetics of such processes. In this study, molecular dynamics (MD) calculations based on density functional theory (DFT+U) were applied to simulate the incorporated Fe in bulk montmorillonite and to explain the measured Fe K-edge X-ray absorption fine structure (XAFS) spectra. The analysis of the experimental data and simulation results suggested that iron in montmorillonite is preferentially incorporated as Fe3+ into the octahedral layer. The simulations showed that there is no preferential occupation of cis- or trans-sites by Fe2+ and Fe3+ in bulk montmorillonite. A very good agreement between the ab initio simulated and the measured XAFS spectra demonstrate the robustness of the employed simulation approach.

Molecular Mechanism of Moisture-Induced Transition in Amorphous Cellulose.

We investigate the influence of adsorbed water on amorphous cellulose structure and properties, within the full range of moisture content from the dry state to saturation, by molecular dynamics simulation. Increasing water content results in overall swelling, a substantial decrease in stiffness, and higher diffusivity of the water molecules. The obtained sorption curve as well as the range of swelling and weakening are confirmed by experiments. The measured properties undergo a noticeable change at about 10% of moisture content, which suggests that a transition occurs in the porous system, indicating that the sorption process is stepwise. Our analysis of water network formation reveals that the onset of percolation coincides with the moisture content at which a transition in the material properties is observed. An in-depth analysis of the molecular mechanism of hydrogen bonding, van der Waals interactions, and water network in the two regimes enhances the understanding of the adsorption process.

Mobility of Na and Cs on montmorillonite surface under partially saturated conditions.

Cs migration in soils at contaminated sites or in clay-rich backfill of waste disposal sites can take place under partially saturated conditions. To understand the molecular mechanism of Cs migration in partially saturated clays, Grand Canonical Monte Carlo simulations were applied to model adsorption of water films onto external surfaces of Cs and Na montmorillonites as function of partial water pressure. The surface complexation and diffusivity of Cs and Na at different partial water pressure was obtained by molecular dynamics simulations. The results suggest that ion mobility in adsorbed water films on external basal surfaces of clay is similar to that in the near-surface water of a saturated pore as far as the thickness of the adsorbed water film is more than two water layers. At lower partial water pressure (i.e., in thinner water films) the ion mobility dramatically decreases. In contrast, the average water mobility in thin water film is higher than in the water-saturated system due to enhanced mobility of water molecules close to vapor-film interface. The results of the simulations were applied to interpret recent laboratory measurements of tritiated water and Cs diffusivity in Callovo-Oxfordian Claystones under partially saturated conditions.