Iron Nitride Nanoparticles for Enhanced Reductive Dechlorination of Trichloroethylene.

Nitriding has been used for decades to improve the corrosion resistance of iron and steel materials. Moreover, iron nitrides (FexN) have been shown to give an outstanding catalytic performance in a wide range of applications. We demonstrate that nitriding also substantially enhances the reactivity of zerovalent iron nanoparticles (nZVI) used for groundwater remediation, alongside reducing particle corrosion. Two different types of FexN nanoparticles were synthesized by passing gaseous NH3/N2 mixtures over pristine nZVI at elevated temperatures. The resulting particles were composed mostly of face-centered cubic (γ'-Fe4N) and hexagonal close-packed (ε-Fe2-3N) arrangements. Nitriding was found to increase the particles' water contact angle and surface availability of iron in reduced forms. The two types of FexN nanoparticles showed a 20- and 5-fold increase in the trichloroethylene (TCE) dechlorination rate, compared to pristine nZVI, and about a 3-fold reduction in the hydrogen evolution rate. This was related to a low energy barrier of 27.0 kJ mol-1 for the first dechlorination step of TCE on the γ'-Fe4N(001) surface, as revealed by density functional theory calculations with an implicit solvation model. TCE dechlorination experiments with aged particles showed that the γ'-Fe4N nanoparticles retained high reactivity even after three months of aging. This combined theoretical-experimental study shows that FexN nanoparticles represent a new and potentially important tool for TCE dechlorination.

Adsorption and exchange reactions of iodine molecules at the alumina surface: modelling alumina-iodine reaction mechanisms.

Harnessing aluminum oxidation energy requires navigating the particle's passivation shell composed of alumina. The shell is a barrier to aluminum oxidation but can also exothermically react with halogenated species and therefore contribute to the overall energy generated during aluminum particle combustion. Fluorination reactions with alumina have been studied because fluorine is abundant in binder formulations that commonly surround aluminum particles in an energetic mixture. However, iodine has emerged as an alternative halogenated-based binder or oxidizer because iodine gas provides ancillary benefits such as chemical neutralization of biological agents or sterilization of contaminated environments. This study used density functional theory (DFT) calculations to evaluate potential reaction pathways for aluminum-iodine combustion. Relative to fluorinated fragments such as HF and F-, the adsorption energies associated with HI and I- are nearly triple the exchange reaction energy available from fluorination reactions with alumina (-189 and -278 kJ mol-1 for HI and I-, respectively). However, exchange reactions between iodinated species and the alumina surface are energetically unfavorable. These results explain that through adsorption, alumina surface exothermic reactions with iodine are more energetic than with fluorine fragments. Experiments performed with differential scanning calorimetry (DSC) confirm the higher magnitude of energy generated for iodination compared with fluorination reactions with alumina. Additionally, strong adsorption energies can promote synthesis of new shell chemistries. Adsorption in solution will promote alumina dissolution and iodine precipitation reactions to produce hydroxyl complexes and iodinated species synthesized on the surface of the particle, thereby replacing alumina with alternative passivation shell chemistry.

Reaction mechanism for fluorination reactions with hydroxylated alumina sites: Pathways promoting aluminum combustion.

Density functional theory calculations were used to reveal the mechanism for the fluorination reaction of active Lewis acid sites on alumina structures, which is important in understanding the pyrophoric processes involving Al particles. In this reaction, hydroxyl groups of active sites are replaced by fluorine anions. Alumina structures were represented by three aluminum aqua hydroxo clusters (labeled AlOOH), in which the Al atom had different coordination spheres, particularly four, five, or six. The F-bearing molecules HF, CH3F, and CF4 were taken as reactants for the fluorination reactions. The overall reaction was represented by four reaction steps as follows: (i) formation of the reaction complex, (ii) activation of the transition state (TS), (iii) deactivation of the TS with a formation of the product complex, and (iv) its decomplexation to individual products. The active reaction center of the TS structure is four-membered, in which two bonds break heterolytically and two form. The lowest reaction barriers were observed for the HF molecule, while the two other molecules had significantly higher reaction barriers. Similarly, the largest overall reaction energies (in absolute value) were found for HF, especially for the five- and six-coordinated Al centers. While the positive charge on the Al center remained almost constant throughout the reaction steps, large charge changes were observed for carbon bearing molecules with a formation of the carbenium cations in the TS step. Realizing the important role of HF in promoting exothermic reactions will enable new molecular design strategies for transforming energy release properties of aluminum powder fuels.

On glyphosate-kaolinite surface interactions. A molecular dynamic study.

Glyphosate is an important and widely used herbicide, its environmental behaviour being of scientific and public interest. Computational models of clay minerals and their interactions with small organic molecules are valuable in studying adsorption processes at an atomistic resolution. We analysed the adsorption process of glyphosate on kaolinite, a clay mineral with a high abundance in several types of soils (e.g., of subtropical or tropical origin), in terms of the adsorption strength. The molecular interactions are characterized by monitoring the occurrence of hydrogen bonds, the orientation of the molecular dipole relative to the interface and the interaction energy. Two different ionic forms of glyphosate were considered: neutral and anionic (-1). It was shown that the main mechanism of the binding of both glyphosate forms to the aluminol surface of kaolinite is through multiple hydrogen bonds. The standard free energy of adsorption of neutral glyphosate from water solution to the basal octahedral surface of kaolinite was computed at -5 kJ mol-1, whereas for the anionic form this quantity amounted to -14 kJ mol-1. Our finding showed that kaolinite has an important contribution to overall adsorption capacity of soils for glyphosate, specifically in its anionic form. HIGHLIGHTS: The adsorption free energy of glyphosate on a kaolinite surface is quantifiedInteractions are computed by quantum mechanics and by classical force fieldMolecular interactions are characterized in terms of hydrogen bonds and orientationsThe effect of polarization of the medium on the calculations is analysed.

A theoretical study of adsorption on iron sulfides towards nanoparticle modeling.

Surface modification of zero-valent iron (nZVI) nanoparticles, which are frequently used in the removal of chlorinated hydrocarbons from contaminated groundwater, can increase their surface stability without significant loss of reactivity. Sulfidation is a process during which thin iron sulfide phases are formed on nZVI particles. In this work, the adsorption capability of two iron sulfide minerals (mackinawite and pyrite) and ZVI with respect to two small polar molecules (H2O and H2S) and trichloroethylene (TCE) was modeled by using the quantum mechanics (QM) approach. High-level QM methods used on cluster models helped in benchmarking and validation of density functional theory methods used on periodic slab models of the (001) surfaces of iron sulfides and the (111) surface of ZVI. This careful computational treatment was necessary for achieving reliable results because modeled iron containing compounds represent computationally demanding systems. The results showed that adsorption was strongly affected by surface topology, accessibility of surface sites, and the shape of adsorbed molecular species. The mackinawite surface is practically hydrophobic having weak interactions with polar molecules (about -5/-6 kcal mol-1), in contrast to the surfaces of pyrite and ZVI (adsorption energies are about three times larger). On the other hand, the adsorption of weakly polar planar TCE molecule is relatively strong and similar for all three surfaces (in the range of -11 to -17 kcal mol-1). Moreover, it was shown that the dominant component of the adsorption energy of TCE had originated from dispersion interactions, which were less important for small polar molecules.

Molecular modelling of sorption processes of a range of diverse small organic molecules in Leonardite humic acid.

Soil organic matter (SOM) is abundant in the environment and plays an important role in several biogeochemical processes, including microbial activity, soil aggregation, plant growth and carbon storage. One of its key functions is the retention and release of various chemical compounds, primarily governed by the sorption process, which strongly affects the environmental fate of nutrients and pollutants. Sorption largely depends on the composition of SOM, as well as its structure, dynamics and the thermodynamic conditions. Although several approaches are available, experimental characterization of sorption mechanisms is not easy. Computational models for predicting sorption coefficients often require a wealth of experimental data for training and are only applicable to compounds and conditions related to the training dataset. Here, we use molecular dynamics (MD) simulations to study the sorption of a range of small organic compounds. As a model SOM system we use the standard Leonardite humic acid (LHA) sample, which physicochemical properties have recently been characterized computationally in detail. This model allowed us to estimate sorption propensities of the systems at two different hydration levels (water activities close to 0 and 1), showing a remarkable correlation with experimental data. Importantly, this molecular modelling approach based on perturbation free-energy calculations is rigorously derived from statistical thermodynamics and requires no experimental sorption data for training. It is therefore in principle applicable to any SOM model or thermodynamic condition. Moreover, the power of MD simulations to provide high-resolution insight into atomistic and molecular interactions was employed to explore how sorbate molecules associate with the LHA matrix and which contacts they form. The heteroatoms of both sorbate and sorbent play an important role and water molecules are identified as further key players in facilitating the sorption process. HIGHLIGHTS: Modelling of the sorption processes in soil organic matter at atomistic level.Rigorous, physics-based approach applicable to a range of SOM systems and conditions.Remarkable level of matching with experimental data with additional insight into the molecular mechanism.Interactions between the sorbate and local environment strongly affects the sorption process.

Interaction between Cellobiose Dehydrogenase and Lytic Polysaccharide Monooxygenase.

Lytic polysaccharide monooxygenases (LPMOs) are ubiquitous oxidoreductases, facilitating the degradation of polymeric carbohydrates in biomass. Cellobiose dehydrogenase (CDH) is a biologically relevant electron donor in this process, with the electrons resulting from cellobiose oxidation being shuttled from the CDH dehydrogenase domain to its cytochrome domain and then to the LPMO catalytic site. In this work, we investigate the interaction of four Neurospora crassa LPMOs and five CDH cytochrome domains from different species using computational methods. We used HADDOCK to perform protein-protein docking experiments on all 20 combinations and subsequently to select four complexes for extensive molecular dynamics simulations. The potential of mean force is computed for a rotation of the cytochrome domain relative to LPMO. We find that the LPMO loops are largely responsible for the preferred orientations of the cytochrome domains. This leads us to postulate a hybrid version of NcLPMO9F, with exchanged loops and predicted altered cytochrome binding preferences for this variant. Our work provides insight into the possible mechanisms of electron transfer between the two protein systems, in agreement with and complementary to previously published experimental data.

Polarization Effects in Simulations of Kaolinite-Water Interfaces.

Computational models of clay minerals and their interactions with the surrounding medium are highly valuable to study adsorption processes at an atomistic resolution, which may be relevant in different areas such as chromatography, environmental chemistry, and so forth. In this work, we analyzed the effect of the treatment of long-range interactions on the polarization of kaolinite-water interfacial systems in terms of structural, electric and dynamic properties, and hydrogen bonds. When using conventional three-dimensional (3D) Ewald summation, water molecules were more structured on the alumina interface of the kaolinite compared to simulations, in which the periodicity perpendicular to the plane was effectively removed. By applying an external electric field to the latter simulations, we were able to reproduce results obtained with the 3D Ewald summation. We proved that the induced polarization promotes adsorption processes of polar and charged species from aqueous solutions using a glyphosate molecule (widely used herbicide) as a particular example. The polarization phenomena raised from the simulation setup should be accounted for carefully in simulations of adsorption processes, which involve periodic images of polar interfaces. An inappropriate treatment of long-range electrostatic interactions can easily lead to artifacts and/or erroneous results.

Molecular Dynamics Simulations of the Standard Leonardite Humic Acid: Microscopic Analysis of the Structure and Dynamics.

Humic substances (HS) are abundant in the environment and play an important role in a number of biogeochemical processes including microbial activity, soil aggregation, plant growth, the retention and release of nutrients, the environmental fate of pollutants, and carbon storage. They are flexible, relatively small molecules forming supramolecular structures through weak interactions. Despite the great importance of understanding their behavior at the atomic level, computational modeling, a premier high-resolution technique providing great level of detail, has been surprisingly little-employed to study humic substances. Here, we use the recently developed Vienna Soil Organic-Matter Modeler to create representative models of a real HS sample, the standard Leonardite humic acid. Molecular dynamics simulations were used to probe the structure and dynamics of the system at a range of hydration levels. The studied systems were characterized in terms of their physicochemical properties, including density, dielectric properties, hydrogen bonding, etc. Moreover, the strength of sorption was estimated for three small organic compounds: benzaldehyde, propan-2-ol, and acetone. Strikingly, the HS models were validated against experimental data showing a remarkable agreement with calculated properties. Finally, we make the equilibrated models of the standard Leonardite humic acid, together with corresponding force-field parameters, available at the Vienna Soil Organic-Matter Modeler.

Adsorption of the herbicide 4-chloro-2-methylphenoxyacetic acid (MCPA) by goethite.

Interaction between the goethite surface and 4-chloro-2-methylphenoxyacetic acid (MCPA) herbicide was studied using density functional theory (DFT) calculations combined with molecular dynamics (MD). The important step made here lies in the use of a periodic DFT method enabling the study of a mineral surface of different protonation states, in strong contrast with previous molecular modeling studies limited to single protonation state corresponding to the point of zero charge. Different surface OH groups and MCPA proton states were used to mimic the strong effects of pH on the outer- and inner-sphere surface complexes that are theoretically possible, together with their binding energies, and their bond lengths. Modeling both a solvated and a protonated (110) goethite surface provided a major breakthrough in the acidic adsorption regime. An outer-sphere complex and a monodentate inner-sphere complex with the neutral MCPA molecule were found to be the most energetically stable adsorbate forms. MD modeling predicted that the latter forms via the sharing of the carbonyl oxygen between the MCPA carboxylate group and a singly coordinated surface hydroxyl group, releasing an H2O molecule. All the other complexes, including the bidentate inner-sphere complex, had higher relative energies and were therefore less likely. The two most likely DFT-optimized structures were used to constrain a surface complexation model applying the charge distribution multisite complexation (CD-MUSIC) approach. The adsorption constants for the complexes were successfully fitted to experimental batch equilibrium data.

Reductive Dechlorination of Chlorinated Ethenes at the Sulfidated Zero-Valent Iron Surface: A Mechanistic DFT Study.

Sulfidated nano- and microscale zero-valent iron (S-(n)ZVI) has shown enhanced selectivity and reactive lifetime in the degradation of chlorinated ethenes (CEs) compared to pristine (n)ZVI. However, varying effects of sulfidation on the dechlorination rates of structurally similar CEs have been reported, with the underlying mechanisms remaining poorly understood. In this study, we investigated the ß-dichloroelimination reactions of tetrachloroethene (PCE), trichloroethene (TCE), cis-1,2-dichloroethene (cis-DCE), and trans-1,2-dichloroethene (trans-DCE) at the S and Fe sites of several S-(n)ZVI surface models by using density functional theory. Dechlorination reactions were both kinetically and thermodynamically more favorable at Fe sites compared to S sites, indicating that maintaining the accessibility of reactive Fe sites is crucial for achieving high S-(n)ZVI reactivity with contaminants. At Fe sites adjacent to S atoms, the reactivity for CE dechlorination followed the order trans-DCE ≈ TCE > cis-DCE > PCE. PCE degradation was hindered at these sites due to the steric effects of S atoms. At the S sites, the energy barriers correlated with the CEs' energy of the lowest unoccupied molecular orbital in the order PCE < TCE < DCE isomers. Our findings reveal that the experimentally observed selectivity of S-(n)ZVI materials for individual CEs can be explained by an interplay of the varying reactivities of Fe and S sites in CE dechlorination reactions.

Intrinsic Effects of Sulfidation on the Reactivity of Zero-Valent Iron With Trichloroethene: A DFT Study.

Sulfidation represents a promising approach to enhance the selectivity and longevity of zero-valent iron (ZVI) in water treatment, particularly for nanoscale ZVI (nZVI). While previous mechanistic studies have primarily concentrated on the impact of sulfidation on the (n)ZVI hydrophobicity, the fundamental effects of sulfidation on the (n)ZVI reactivity with target contaminants remain poorly understood. Herein, we employed density functional theory to elucidate reaction mechanisms of trichloroethene (TCE) dechlorination at various (n)ZVI surface models, ranging from pristine Fe0 to regularly sulfidated Fe surfaces. Our findings indicate that sulfidation intrinsically hinders the TCE dechlorination by (n)ZVI, which aligns with prior observations of sulfur poisoning in transition metal catalysts. We further demonstrate that the positive effects of sulfidation emerge when the surface of (n)ZVI undergoes corrosion. Notably, S sites exhibit higher reactivity compared to the sites typically present on the surface of (n)ZVI oxidized in water. Additionally, S sites protect nearby Fe sites against oxidation and make them more selective for direct electron transfer. Overall, our results reveal that the reactivity of sulfidated (n)ZVI is governed by an interplay of intrinsic inhibitory effects and corrosion protection. A deeper understanding of these phenomena may provide new insights into the selectivity of sulfidated (n)ZVI for specific contaminants.

Iron nitride nanoparticles for rapid dechlorination of mixed chlorinated ethene contamination.

Sulfidation and, more recently, nitriding have been recognized as promising modifications to enhance the selectivity of nanoscale zero-valent iron (nZVI) particles for trichloroethene (TCE). Herein, we investigated the performance of iron nitride (FexN) nanoparticles in the removal of a broader range of chlorinated ethenes (CEs), including tetrachloroethene (PCE), cis-1,2-dichloroethene (cis-DCE), and their mixture with TCE, and compared it to the performance of sulfidated nZVI (S-nZVI) prepared from the same precursor nZVI. Two distinct types of iron nitride (FexN) nanoparticles, containing γ'-Fe4N and ε-Fe2-3N phases, exhibited substantially higher PCE and cis-DCE dechlorination rates compared to S-nZVI. A similar effect was observed with a CE mixture, which was completely dechlorinated by both types of FexN nanoparticles within 10 days, whereas S-nZVI was able to remove only about half of the amount, most of which being TCE. Density functional theory calculations further revealed that the cleavage of the first C-Cl bond was the rate-limiting step for all CEs dechlorinated on the γ'-Fe4N(001) surface, with the reaction barriers of PCE and cis-DCE being 29.9, and 40.8 kJ mol-1, respectively. FexN nanoparticles proved to be highly effective in the remediation of PCE, cis-DCE, and mixed CE contamination.

Adsorption of the hydrophobic organic pollutant hexachlorobenzene to phyllosilicate minerals.

Hexachlorobenzene (HCB), a representative of hydrophobic organic chemicals (HOC), belongs to the group of persistent organic pollutants (POPs) that can have harmful effects on humans and other biota. Sorption processes in soils and sediments largely determine the fate of HCB and the risks arising from the compound in the environment. In this context, especially HOC-organic matter interactions are intensively studied, whereas knowledge of HOC adsorption to mineral phases (e.g., clay minerals) is comparatively limited. In this work, we performed batch adsorption experiments of HCB on a set of twelve phyllosilicate mineral sorbents that comprised several smectites, kaolinite, hectorite, chlorite, vermiculite, and illite. The effect of charge and size of exchangeable cations on HCB adsorption was studied using the source clay montmorillonite STx-1b after treatment with nine types of alkali (M+: Li, K, Na, Rb, Cs) and alkaline earth metal cations (M2+: Mg, Ca, Sr, Ba). Molecular modeling simulations based on density functional theory (DFT) calculations to reveal the effect of different cations on the adsorption energy in a selected HCB-clay mineral system accompanied this study. Results for HCB adsorption to minerals showed a large variation of solid-liquid adsorption constants Kd over four orders of magnitude (log Kd 0.9-3.3). Experiments with cation-modified montmorillonite resulted in increasing HCB adsorption with decreasing hydrated radii of exchangeable cations (log Kd 1.3-3.8 for M+ and 1.3-1.4 for M2+). DFT calculations predicted (gas phase) adsorption energies (- 76 to - 24 kJ mol-1 for M+ and - 96 to - 71 kJ mol-1 for M2+) showing a good correlation with Kd values for M2+-modified montmorillonite, whereas a discrepancy was observed for M+-modified montmorillonite. Supported by further calculations, this indicated that the solvent effect plays a relevant role in the adsorption process. Our results provide insight into the influence of minerals on HOC adsorption using HCB as an example and support the relevance of minerals for the environmental fate of HOCs such as for long-term source/sink phenomena in soils and sediments.

Edge-Site-Free and Topological-Defect-Rich Carbon Cathode for High-Performance Lithium-Oxygen Batteries.

The rational design of a stable and catalytic carbon cathode is crucial for the development of rechargeable lithium-oxygen (LiO2 ) batteries. An edge-site-free and topological-defect-rich graphene-based material is proposed as a pure carbon cathode that drastically improves LiO2 battery performance, even in the absence of extra catalysts and mediators. The proposed graphene-based material is synthesized using the advanced template technique coupled with high-temperature annealing at 1800 °C. The material possesses an edge-site-free framework and mesoporosity, which is crucial to achieve excellent electrochemical stability and an ultra-large capacity (>6700 mAh g-1 ). Moreover, both experimental and theoretical structural characterization demonstrates the presence of a significant number of topological defects, which are non-hexagonal carbon rings in the graphene framework. In situ isotopic electrochemical mass spectrometry and theoretical calculations reveal the unique catalysis of topological defects in the formation of amorphous Li2 O2 , which may be decomposed at low potential (∼ 3.6 V versus Li/Li+ ) and leads to improved cycle performance. Furthermore, a flexible electrode sheet that excludes organic binders exhibits an extremely long lifetime of up to 307 cycles (>1535 h), in the absence of solid or soluble catalysts. These findings may be used to design robust carbon cathodes for LiO2 batteries.

The stability of the acetic acid dimer in microhydrated environments and in aqueous solution.

The thermodynamic stability of the acetic acid dimer conformers in microhydrated environments and in aqueous solution was studied by means of molecular dynamics simulations using the density functional based tight binding (DFTB) method. To confirm the reliability of this method for the system studied, density functional theory (DFT) and second order Møller-Plesset perturbation theory (MP2) calculations were performed for comparison. Classical optimized potentials for liquid simulations (OPLS) force field dynamics was used as well. One focus of this work was laid on the study of the capabilities of water molecules to break the hydrogen bonds of the acetic acid dimer. The barrier for insertion of one water molecule into the most stable cyclic dimer is found to lie between 3.25 and 4.8 kcal mol(-1) for the quantum mechanical methods, but only at 1.2 kcal mol(-1) for OPLS. Starting from different acetic acid dimer structures optimized in gas phase, DFTB dynamics simulations give a different picture of the stability in the microhydrated environment (4 to 12 water molecules) as compared to aqueous solution. In the former case all conformers are converted to the hydrated cyclic dimer, which remains stable over the entire simulation time of 1 ns. These results demonstrate that the considered microhydrated environment is not sufficient to dissociate the acetic acid dimer. In aqueous solution, however, the DFTB dynamics shows dissociation of all dimer structures (or processes leading thereto) starting after about 50 ps, demonstrating the capability of the water environment to break up the relatively strong hydrogen bridges. The OPLS dynamics in the aqueous environment shows--in contrast to the DFTB results--immediate dissociation, but a similar long-term behavior.

Assessment of ten DFT methods in predicting structures of sheet silicates: importance of dispersion corrections.

The performance of ten density functional theory (DFT) methods in a prediction of the structure of four clay minerals, in which non-bonding interactions dominate in the layer stacking (dispersive forces in talc and pyrophyllite, and hydrogen bonds in lizardite and kaolinite), is reported. In a set of DFT methods following functionals were included: standard local and semi-local (LDA, PW91, PBE, and RPBE), dispersion corrected (PW91-D2, PBE-D2, RPBE-D2, and vdW-TS), and functionals developed specifically for solids and solid surfaces (PBEsol and AM05). We have shown that the standard DFT functionals fail in the correct prediction of the structural parameters, for which non-bonding interactions are important. The remarkable improvement leading to very good agreement with experimental structures is achieved if the dispersion corrections are included in the DFT calculations. In such cases the relative error for the most sensitive lattice vector c dropped below 1%. Very good performance was also observed for both DFT functionals developed for solids. Especially, the results achieved with the PBEsol are qualitatively similar to those with DFT-D2.

Understanding of bonding and mechanical characteristics of cementitious mineral tobermorite from first principles.

This paper reports density functional theory study of the structural and mechanical properties of tobermorite mineral (9 Å phase) as one of the main component of cementitious materials in a concrete chemistry. Calculated bulk modulus and elastic constants reflect a relatively high resistance of the tobermorite structure with respect to external isostatic compression. Moreover, the elastic constants proved the anisotropic character of the tobermorite structure. The directions parallel to the axb plane are more resistant to the compression than the perpendicular direction. The largest contribution to this resistance comes from the "dreierketten" silicate chains. The bonding analysis linked macroscopic mechanical properties and the atomic structure of the tobermorite. It was found that polar covalent Si-O bonds are stiffer than iono-covalent Ca-O bonds. The SiO(4) tetrahedra are resistant with respect to the compression and the effect of external pressure is reflected by the large mutual tilting of these tetrahedra as it is shown by changes of the Si-O-Si bridging angles. Polyhedra with the seven-fold coordinated Ca(2+) cations undergo large structural changes. Especially, axial Ca-O bonds perpendicular to the axb plane are significantly shortened. Besides, it was shown that structural parameters, more or less in parallel orientation to the axb plane, are mainly responsible for the high resistance of the tobermorite structure to external pressure. The main mechanism of a dissipation of energy entered to the structure through the compression is proceeded by the tilting of the tetrahedra of the silicate chains and by large shortening of the axial Ca-O distances.

Molecular dynamics simulations of water molecule-bridges in polar domains of humic acids.

The stabilizing effect of water molecule bridges on polar regions in humic substances (HSs) has been investigated by means of molecular dynamics (MD) simulations. The purpose of these investigations was to show the effect of water molecular bridges (WAMB) for cross-linking distant locations of hydrophilic groups. For this purpose, a tetramer of undecanoid fatty acids connected to a network of water molecules has been constructed, which serve as a model for spatially fixed aliphatic chains in HSs terminated by a polar (carboxyl) group. The effect of environmental polarity has been investigated by using solvents of low and medium polarity in force-field MD. A nonpolar environment simulated by n-hexane was chosen to mimic the stability of WAMB in a hydrophilic hotspot surrounded by a nonpolar environment, while the more polar acetonitrile environment was chosen to simulate a more even distribution of polarity around the carboxylic groups and the water molecules. The dynamics simulations show that the rigidity of the oligomer chains is significantly enhanced as soon as the water cluster is large enough to comprise all four carboxyl groups. Increasing the temperature leads to evaporization processes which destabilize the rigidity of the tetramer-water cluster. Embedding it into the nonpolar environment introduces a pronounced cage effect which significantly impedes removal of water molecules from the cluster region. On the other hand, a polar environment facilitates their diffusion from the polar region. One important consequence of these simulations is that although the local water network is the stabilizing factor for the organic matter matrix, the degree of stabilization is additionally affected by the presence of nonpolar surroundings.

Ab initio study of structure and interconversion of native cellulose phases.

Dispersion-interaction corrected DFT simulations are performed to study the structure of two allomorphs of native cellulose I. Good agreement between theoretical and experimental data is achieved. Two H-bond patterns, previously identified experimentally, are examined for both allomorphs. The transition mechanism for the conversion between the phase I(α) and I(ß) is studied by means of constrained relaxations. New metastable intermediate phase occurring on the I(α) → I(ß) route is identified, and the corresponding structural data are reported.