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.

Background: Knowledge of the stabilizing mechanisms of soil organic matter (SOM) is extremely important for numerous soil functions. For this, insight into the nature of organic matter through appropriate model concepts are crucial.Aims: For several years, a heated debate has emerged on the transformation and stabilization of SOM. In the present work, we try to contribute to this debate using molecular modeling and providing a comprehensive overview of the history of application of molecular modeling tools and developing structural concepts of SOM.Methods: Molecular modeling methods based on quantum and/or classical mechanics were used to model SOM and related properties including interactions with reactive surfaces of soil minerals.Results: Modeling of SOM aggregates revealed that hydrogen bonds and cation bridges are the main stabilizing factors in soil solution, whereas pH modifies the stability. The modeled supramolecular SOM aggregates exhibit physicochemical properties, similar to those of humic substances (HS) described in literature. The interactions of the HS models with surfaces in kaolinite nanopores led to a partial disintegration of the aggregates into individual molecules and/or smaller subaggregates.Conclusions: From the molecular modeling point of view, supramolecular microaggregate models that exhibit the properties of HS are stable in the soil solution. However, their binding to reactive mineral soil constituents can be also in the form of individual molecules or subaggregates. Thus, HS microaggregate stability is relative, depending on the interacting environment. This reconciles two points of view of HS: either as small molecules and/or supramolecular structures.

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 ⁻¹ , whereas for the anionic form this quantity amounted to −14 kJ mol ⁻¹ . 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 quantified Interactions are computed by quantum mechanics and by classical force field Molecular interactions are characterized in terms of hydrogen bonds and orientations The effect of polarization of the medium on the calculations is analysed

Soil organic matter (SOM) and various inorganic minerals represent key components of soils. During pedogenesis and due to biological activity these species interact, having a crucial impact on the formation of an aggregated soil structure with a hierarchical arrangement from nano to macro scale. In this process, the formation of organo–mineral microaggregates represents a dominant factor affecting soil functions and properties. This study focuses on the interactions between humic substances (HSs) and the mineral kaolinite as typical representatives of SOM and soil minerals. By performing classical molecular dynamics (MD) simulations on models of HSs and kaolinite, we demonstrate how two dominant but chemically different kaolinite surfaces affect the stability of HSs microaggregates. By analyzing volumetric, structural, and energetic properties of SOM–kaolinite models, we explain possible mechanisms of the formation of stable SOM–clay aggregates and show how a polarized environment affects the electrostatic interactions, stabilizing the microscopic structure of SOM–mineral aggregates. Our results showed that when stable aggregates of HSs are confined in kaolinite nanopores, their interactions with kaolinite surfaces disintegrate them into smaller subaggregates. These subaggregates are adsorbed more strongly on the polar aluminol surface of kaolinite compared to less the active hydrophobic siloxane surface.